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ABSTRACT 

A computer code for predicting incident thermal radiation from defined 
plume gas properties in either axl symmetric or cylindrical coordinate systems. 
The radiation model is a statistical band model for exponential line strength 
distribution with Lorentz/Doppler line shapes for 5 gaseous species (H^O* CO^, 
CO, HCl and HF) and an approximate (non-scattering) treatment of carbon 
particles. The Curtis-Godson approximation is used for inhomogeneous gases, 
but a subroutine is available for using Young's Intuative Derivative method 
for HgO with Lorentz line shape and exp. -tailed-inverse line strength distri- 
bution. The geometry model provides integration over a hemisphere with up 
to 6 individually oriented identical axisymmetric plumes, a single 3-D plume, 
or a combined 3-D/axi symmetric plume for approximating clusters at high alti- 
tudes. Shading surfaces may be used in any of 7 shapes, and a conical limit 
may be defined for the plume to set individual line-of«sight limits. Inter- 
mediate coordinate systems may be specified to simplify input of plumes and 
shading surfaces. 
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1 INTRODUCTION 

The 6ASRAD code is the latast in a long series of evolutionary develop- 
ments by NASA/MSFC to provide improved prediction of radiation transfer from 
rocket exhaust plumes. Some of the coding is unchanged from that for the 
earliest version of the program* (Ref. 1) for single constituent axisymmetric 
plumes. This initial code was followed by a modification to include four 
radiating species (HjO, CO 2 , CO and carbon) (Ref. 2), a three-dimensional 
geometry for a single species (Ref. 3), and a ccmibined code for four species 
using either axisymmetric or three-dimensional plume property formats (Refs. 4 
and 5). This final version was expanded to Include HCl and HF with provisions 
for adding a sixth gaseous species if desired (Ref. 6). 

The current code uses a statistical band model with an exponential line 
strength distribution for combined Lorentz and Doppler line shapes. It uses 
the Curtis-Godson approximation for inhomogeneous gases, so it can be easily 
modified for approximating an exponentially- tailed-inverse line strength distri- 
bution. A subroutine is included incorporating Young's Intuitive Derivative 
method (Ref. 7) for HjO, but the subroutine available in earlier versions for 
modeling radiation using multiple line groups to represent several energy levels 
has been omitted. Other changes in the radiation portion of the program include 
the ability to input band model data rather than using it in data statements to 
allow comparisons and use of the most recent data. 

Significant changes have been made in the geometry procedures, but the code 
is still designed for heat transfer applications using lines-of-sight representing 


* Original coding by A. L. McGarrity at the MSFC Computation Laboratory. 
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a hemisphere over a point rather than detection applications using parallel 
lines-of-sight through the plume. Shading surfaces can be described using 
seven geometric shapes, and plumes can have conical boundaries specified which 
automatically adjust line-of-sight limits to the region inside the cone. 

Several identical, independently oriented, axi symmetric plumes can be specified 
if it is assumed they do not interfere, and a three-dimensional plume descrip- 
tion is available for describing flov/fields which are not axi symmetric. 

The following sections describe the general organization of the GASRAD code, 
the radiation modeling procedures, and the geometry modeling procedures. The 
Appendix Includes operating information such as input and output formats, sample 
problems, and descriptions of the program elements. The operating descriptions 
are intended specifically for the MSFC UNIVAC 1108. But the program is coded 
entirely in FORTRAN, so it may be adapted to other systems with minor changes. 
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2 GENERAL ORGANIZATION 

This section introduces the overall problem geometry which will be 
referred to in following sections. In addition, the general arrangements 
of the program design will be outlined to both orient the reader and define 
the goal of the overalay structure with the limitations inherent in the 
approach chosen. 

All geometry is related to a central (or reference) coordinate system 
(CCS) which may be conveniently located to suit the problem. For example, 

Space Shuttle predictions were made with the orbiter coordinate system, so 
surface points could be easily located from vehicle drawings. However, if a 
single engine on a test stand is analyzed, it is usually most convenient to 
choose the central system at the engine gimbal point or the center of the 
exit plane (if no gimballing is planned). Intermediate coordinate systems 
(ICS) may be located in the CCS by specifying a location and three axis rota- 
tions. Then the point-of-interest (POI), plumes, and shading surfaces can be 
located in either the CCS or a designated ICS by specifying a location and 
three axis rotations. 

A typical arrangement is illustrated in Fig. 2.1. In this illustration an 
intermediate system is used to facilitate location of the plume and the shading 
surface representing the nozzle. Use of a gimbal point ICS also provides a 
convenient means of locating surfaces on the nozzle and providing gimballing of 
the entire arrangement by rotating the ICS, The UVW coordinate system at the 
point-of-interest is used to define a hemispherical system for locating lines- 
of-sight. Positions are specified by elevation and azimuth angles 6 and <t> relative 
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Figure 2.1 Overall Geometry Arrangement 
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to the surface normal W and the U axis. If a conical plume boundary is 
specified, the integration limits on each line-of-sight (SZERO* and SMIN) 
are defined by the program within the limits SZERO and SMAX specified as 
input. 

In the general operation of the program, two large blocks of data 
representing the flowfield gas properties and the band model parameters 
must be stored for computation. Storing both data sets simultaneously 
would require a large memory, a limited flowfield description, or a rela- 
tively small spectral range, so the program was separated. Flowfield data 
are loaded, and as gas property interpolations are completed for each line- 
of-sight, it is written on an external storage device. When the flowfield 
interpolations are complete, the band model data are loaded in place of the 
flowfield, and gas properties are read in one line-of-sight at a time for 
the band model radiation calculation. This procedure does not use excessive 
CPU time, but increasing costs associated with use of external stores com- 
bined with decreasing internal storage costs may make other procedures 
more attractive for future applications. 

General features of the program organization and operation are illus- 
trated by the schematic in Fig. 2,2 which omits details within the flowfield 
interpolation and band model prediction blocks. Program control and integra- 
tion limits and intervals are input in a small MAIN subroutine which controls 
the top level of program overlay. MAIN passes control to FLOW which controls 
overlay of subroutines used in the flowfield Input and interpolation. In 
these subroutines, the remaining geometric data are read and the problem limits 
and geometry are output. Flowfield properties may be read from cards or a 
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binary data file. But if they are input via cards, a binary data file is 
prepared before the interpolation process begins. As the flowfield proper- 
ties are input, they are printed in either a long or short form to verify 
and record input. 

After gas property data along each line-of-sight has been transferred 
to the external data file, MAIN transfers control to RAD which controls over- 
lay of subroutines in the band model prediction. Band model data are input 
in card format to simplify additions or changes to the data, but in normal 
operation the card formatted data are stored on tape with the program and 
added to the runstream by control system corranand. The band model parameters 
are printed using either a short or long form to verify and record input. 

As radiation on each line-of-sight is calculated, details of the results and 
integration limits are printed, and the position {e,ii>) and radiance are trans- 
ferred to an external store for later use by VIEW. Upon completion of all 
1 ines-of-sight, the total flux for the integration is printed followed by 
the spectral distribution. Spectral data also may be transferred to an 
external file, if requested, for use by the spectral plot program. Control 
returns to MAIN which calls the VIEW subroutine. 

VIEW prepares a printer-plot representing the intensity in a polar format. 
The radius in the plot represents the elevation angle 6, and the angle repre- 
sents the azimuth angle <(). The digits 0-9 are used to represent truncated values 
of radiance/peak-radiance. If lines are blocked before entering the plume they 
are represented by B, while if they miss the plume, the position is blank. 

Some detail discrepancies occur in averaging the results into the rectangular 
grid required for printing, but the larger scale characteristics shown help to 
verify the point orientation and integration limits and illustrate the major 
radiation sources. 
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3 RADIATION MODEL 

Although the problem geometry is handled before the radiation calcu- 
lation in the computer code, the radiation model will be discussed first 
to establish the gas properties necessary and the options used for spatial 
integration which both affect the geometry procedures. This section will 
describe the band model formulations used, the available band model data, 
and considerations in selecting spectral and spatial intervals for integra- 
tion. 

3.1 BAND MODELS 

Band models have been developed as a simplification to monochromatic 
(line-by-line) calculation methods which become difficult and expensive in 
computing time as the number of spectral lines to be included becomes large. 

The objective of a band model is to provide an approximate mean value of 
transmission over a small spectral interval (usually 5 to 25 cm“^) using a 
few parameters to describe the behavior of lines in the interval. The following 
discussion will: briefly introduce concepts used in describing gaseous radiation, 
introduce the homogeneous gas band models used in the GASRAD program, describe 
inhomogeneous gas approximations, and summarize the radiation prediction models 
available in the GASRAD program. 

Gaseous Radiation Fundamentals 

A basic parameter in gaseous radiation is the monochromatic absorption 
coefficient, k(w,s) defined by 

dN(u),s)/ds = -k(m,s) p(s) N(o),s) (1) 

PREGEDIMO PAGE MOT FILMED 
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where N(w,s) is the spectral radiance at the wavenumber* w and position s, 
and p(s) is the concentration of the radiating species per unit path length. 

The functional dependence of k and p on s implies a dependence on gas property 
variations for a general case. For a homogeneous gas» these two parameters are 
not a function of s, but are determined by the gas temperature, pressure, and 
species concentration. Integration of Eq. 1 over a path length s defines the 
transmissivity 

t(w,s) = N(oo,s)/N(oo,o) (2) 

s 

= exp[-/ k(o),s)p(s)ds] 

which, for a homogeneous gas, is 

t(w,u) = exp[-k(w)u] (3) 

where u represents the concentration, ps, in units inverse to those of k. 

Spectral lines representing energy transitions in a molecule are broadened 
by collision and Doppler effects. The collision broadening is adequately 
described by the Lorentz line shape while the broadening due to thermal motion 
of the molecules is represented by the Doppler line shape. The line-strength 
or intensity describes the integrated absorption coefficient over the broadened 
line 

CO 

S = / k(o3)dw (4) 

0 

while the integrated absorptance is a function of both k and u 

oo oo 

W = / a{{jL),u)da3 = / [1 -T(w,u)]da) 

0 0 

[l-exp(-k(w)u)]dw (5) 

This is usually termed the equivalent-width since it represents the spectral 
width of an equivalent black (a=l ) line. The variation of W with u is the 
curve-of-growth for the line. 


* Wavenumber (cnrM is the inverse of wavelength (ym), so w = lOV^. 
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Representative line shapes and curves-of-growth are illustrated in Fig. 3.1 
for a single line. In the weak-line limit (Su/Y«l)the curves-of-growth 
approach a linear limit. As the absorption increases, the two line shapes 
produce different results. The Lorentz shape, which has significant absorption 
in the wings of the line, approaches a square-root limit, while the narrower 
Doppler line produces a much lower equivalent width. which increases as 
Un(Su/y)]^. 

Band Models for Homogeneous Gases 

Band models are generally characterized by assumptions made for the line strength 
(or intensity) distribution. They are commonly divided into three classes: 
the regular model, the statistical (or random) models, and the mixed models. 

The regular model considers absorption by identical, equally spaced lines. 

This approximates the spectra of many diatomic molecules such as HF or NO. 
Statistical models assume that lines are randomly distributed with a specified 
intensity distribution. With proper choices of intensity distribution functions, 
random models approximate features of polyatomic molecules, such as HgO and 
CO 2 . Mixed models interpolate between the limits of regularity and randomness 
of the other two models. For example, a mixed model may use a random super- 
position of regular models. 

Three parameters are used in band models to characterize the line structure 
and the resulting variation of mean transmissivity with path length in a small 
spectral interval. Representative line half-widths and spacing, y and d, are 
used to characterize the spectral structure varying from isolated lines (y « d) 
to strongly overlapping lines (y » d), and the mean line strength is represented 
by the absorption coefficient k" which is proportional to the line-strength- to- 
spacing ratio. S/d. 
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(w-(0i) /y 

(a) Line Shape (y = line half-width at half height) 



Su/Y 

(b) Curve-of-Growth 


Figure 3.1 Line Shape and Curve-of-Growth for Lorentz and Doppler Lines. 
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The GASRAD program uses statistical band models to best describe HjO 
and COz which are the gases of primary interest in heat transfer, but this 
model is not usually expected to provide the best representation of diatomic 
gases. Statistical band models have been developed based on four assumptions 
concerning the probability of occurrance, P(S), of a line strength, S: 

A) Equal strength, P(S)a'S{S-^, 

B) Exponential distribution, P{S)«x exp(-S), 

C) Inverse distribution, P(S)oc 1/S, 

D) Exponentially- tailed inverse, P(S)cc S"^exp(-S), 

where 5 in (A) refers to the Dirac delta function*. These distributions have 
been studied for both Lorentz and Doppler line shapes. For the Lorentz shape, 
the results can be represented by 

(-In t) = 3l T(x) (6) 

where = 2iTYL/d, x = T<u/3l» and f(x) is determined by the probability function 
used. Because of variations in constants, a similar generalized Doppler rela- 
tion must be stated in proportional form 

(-in t)oc 3d g(y) (7) 

where 3D“'n‘YD/'^» y “Bi/3d» and g(y) is determined by the intensity distribution. 

The parameters F and 1/d are usually tabulated as functions of wave- 
number (w) and temperature (T), while algebraic expressions are available for 
Y|_ and yq* The expression normally used for Y|_ (Ref. 10) is 

U1 ' [Ci Y1l(273/Tr'1 + I CJ Y1j(273/T)'”'>] (8) 

where Y-fi represents resonant self broadening of the ith species and Y-|j 
represents broadening by foreign gases and nonresonant self broadening. The 
pressure, temperature, and mole fraction are indicated by p, T, and c, and 
estimates for the coefficients (y) and exponents (n) are given in Table 3.1. 

An expression for Doppler broadening can be derived from the Boltzman distribu- 

*6(S-S)=0 when S/S and /6(S-S)dS-l when the integration limits include S. 
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tion for molecular velocities. With the constants evaluated, it is 

YDi = 0.3581 X 10“® w(T/M^)^ (9) 

where represents the molecular mass in unified atomic mass units (e.g. 
gm/gm-mole). 

For given values of IT, y» and d (determined by species concentrations, 
pressure, and temperature), the functions f(x) or g(y).can be viewed as repre- 
senting the relative variation of (-In t) with u. This variation is termed the 
"curve-of -growth" for the band model which is illustrated for the various statis- 
tical distributions in Fig. 3,1. These curves-of -growth are similar to those 
for a single line, and for statistical band models, -In t represent a mean 
equivalent-width (iT) multiplied by the number of lines in the interval (1/d), 
W/d. Although the curves are similar to those for a single line, there are 
differences in implications caused by modeling a group of lines. For band 
models, the linear limit can be approached by strongly overlapping lines or weak 
lines, so that the net result is Fu « 8 l- Conversely, the square root limit 
requires a relatively small line-width to spacing ratio, 3l «Fu. As a result, 
the linear limit is approached for high pressures (which produce large Y[_) or 
high line densities (such as in the 4.4pm band of CO 2 ). The Doppler line 
curves-of-growth also approach the linear Unit for small y, but approach two 
different limits for large y. 

The differences between the Lorentz curves-of-growth are not large, so 
Models B and D which have simple algebraic solutions are often used as approxi- 
mations of Models A and C. Which have more complex solutions. All of the 
Doppler line models require lengthy series solutions, so two algebraic approxi- 
mations have been used: one for solutions E(y) and F(y) , and one for solutions 

G(y) and H(y). In the Doppler case, the model chosen could make a significant 
difference for large y in cases of relatively hot low-pressure gases where 
Doppler broadening is significant. 
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The band models considered for use in the GASRAD program are represented 
by the following approximations from Ref. 10 with aj_ = yl/^ and 30 = yd/^- 


Exponential Line Strength 

Lorentz (W/d)^ = X|_ = Fu//l+¥u/(4aL) (10a) 

Doppler (W<1)d = Xd = 1-70 ap «K[l+(Fu/(1.70ao))M}’^ (10b) 

Exponentially-Tailed-Inverse Line Strength 

Lorentz (lT/d)j_ = X|_ = 2aL [(l+Fu/aL)M]. (11a) 

Doppler (W/d)o = Xq = O.937a0{Ml+(Fu/(O.937))®'®]}^'2 


Either of these equation sets are applied to the i radiating species in the 
gas. Then the Lorentz and Doppler components are combined to compute to trans- 
missivity using the form 

ZnT - -I k-j u-j (1 - y-j (12a) 

yi - [1 - (Xi/wr^ + Cl - {X0/ku)2]“"-i (12b) 

For a single species, these equations will reduce to In t -X^ if X[_ > 3 Xq 

or tn X -Xq if Xp > 3X|_, and will give the linear limit. In x - -Fu, if 

X^ = Fu or Xq Fu. Representative values of optical depth for HjO are 

illustrated in Fig. 3.3 showing the effect of pressure on the curve-of -growth 

and the relative importance of the Doppler broadened optical depth. 

Radiance along a line-of-sight t to a point at s=o is given by 

N((o,L) =-/ N°(oj,s)[d7(w,s)/ds]ds (13) 

0 

where N®(u,s) is the Planck function evaluated at w and T(s). The heat trans- 
fer to a point is obtained by integrating Eq. 13 over a spectral interval 
adequate to include the radiating species and a spatial region which includes 
the radiating gas. The result is 

q/A {/ N‘*(a),s)[dT(ca,s)/ds]ds}sinecos0dajd<}ide (14) 

where 6 and <p are the elevation and azimuth angles measured from the surface 

normal. In numerical evaluation of Eq. 14, Aw would correspond to the 
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spectral interval used in the band model representation. 



rL - atm-meter (for H^O at 2000K) 


Figure 3.3 Representative Optical Depth Behavior of 
H 2 O at 2.88 pm. 


Band Models for Inhomogeneous Gases 


Application of homogeneous gas band models to inhomogeneous gases is 
made using one of three techniques: multipling transmissivities, 
computing effective values of band model parameters using the modified 
Curtis-Godson approximation, or using derivative methods designed with the 
possibility of inhomogeneity included in the band model formulation. Multi- 
plication of transmissivities, which is valid for gray media, is not valid 
in all cases for band models. The mean transmissivity representation of line 
structure by a band model is valid when applied to gray sources or sources with 
line structure which is uncorrelated, such as a different gas. But for trans- 
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mission through the same gas at a different temperature, some of the lines 
can be expected to correlate and produce more absorption than predicted by 
simple multiplication of homogeneous band model results. 

The modified Curtis-Godson approximation uses the homogeneous band 
model formulation with effective parameters defined by 

s 


k (o),i ,s)u(i ,s)=/ k(w,i ,s‘)c(i ,s')p(s')ds' 


(15a) 


aLe(w,i,s)=[/ k{w,i,s')c(1,s')p{s')a|_(w,i,s')ds']/[kg(w,i,s)u(i,s)](15b) 

s__ _ 

aQg(w,i,s)=[ / k(ai,i,s')c(i,s')p(s')ap{w,i,s')ds]/[kg(a),i,s)u(1.s)] (15c) 


where c(i,s') is the mole fraction of the ith species and pCs") is the pressure. 
The product c*p represents the ratio of density at s' to the density at one 
atmosphere and a temperature equal to that at s'. This requires absorption 
coefficients defined relative to density at the gas temperature. In some 
applications, k is referenced to standard conditions (p=l atm, Tq= 273K), so 
the concentration must be represented by cp(To/T). 

Derivative models were summarized by Young (Ref. 7), and the Intuitive 
Derivative method he recommended on the basis of H 2 O predictions will be 
summarized here. The recommended formulation is a representation of the 
inverse or exponentially- tailed-inverse intensity distribution for Lorentz 
lines (equivalent to Eq. 11a for homogeneous gas). The procedure is compared 
with the Curtis-Godson approximation in Table 3.2 where the functional notation 
of wavenumber (w) has been omitted, and the parameter 3=2TTY|_/d is used rather 
than a]_=YL/‘^ in Eq. 11a, An efficient algorithim for approximating y(s) is 
provided by Young (Ref. 7), so the additional parameters are not difficult to 
compute. But the computation time can be greatly Increased because of the 
decrease in integration step size required in the Intuitive Derivative method. 


19 



F=3^K/T-r^O>-H IMO. 


RTR 014-9 


Table 3.2 


COMPARISON OF CURTIS-GODSON AND INTUITIVE DERIVATIVE METHODS* 
FOR AN INHOMOGENEOUS GAS WITH A SINGLE RADIATING SPECIES 




Same 

(16a) 

u(s)=/^c(s')b(s')ds‘ 

0 

(17a) 



Same 

(16b) 

(16c) 

1 S 

ke(s)=u(s) 1 c(s')p(s')k(s')ds'(17b) 

I s 



Same 

i Cc(s')p(s')k(s') 





G(s')ds'] 

(17c) 



Same 

(16d) 

Xe(s)=u(s)Fg(s)/3g(s) 

(17d) 



N/A 


u(s)kg(s) ^ [c(s* )p(s ‘ )k(s *) 





d(s')ds'] 

(17e) 



N/A 


r(s)=B(s)/Be(s) 

(17f) 1 



N/A 


q(s)=dg(s)/d(s) 

(I7g) 



N/A 


y(s)=y(iTXe{s),r(s),q(s)) 

(17h) 



N/A 


. 9 “ / 2xn^ 1 

y(x.r.q)-f / exp ^ 






1 ^2 
dnJl +22 

(17i) 



N/A 


=c(s)p(s)k(s)y(s) 

d ds 

(17j) 

w 

6e(s) 

7T 



K = / 1 dW(s‘) ■ . 
d 0 d ds 

(17k) 

d 

= [/ l+27rxQ(s)-l] (16k) 



N/A 


$L = _:p(s) 1 dW(s) 
dfs ^ ' d ds 

(171) 



Same 


T = exp{-W/d) 

(17m) 

N= 

0 

N° d7(s') 

(16n) 

N=-/^ N°(s') ds' 

0 ds' 

(17n) 


* The Intuitive Derivative model and Eq. 16k of the C-G method are 
approximations for an inverse or exponentially-tailed-inverse 
Lorentz line strength f"' ’stribution. 
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Numerical evaluation of radiance with the Curtis-Godson approximation 
takes the form 
_ *’’’ 

% = ^ ”n ( Vl - (16) 

n=l 

with Tq'!, and the size of the integration step can be arbitrarily large as 
long as the gas in the step can be assumed to be homogeneous. In contrast, 
the value of y(s) changes continuously with Xg even, in homogeneous sections, 
so small steps must be used to properly evaluate W/d (Eq. 17k). In the 
GASRAD code for the Intuitive Derivative method, an input parameter is 
provided to vary the step size to give a specified value of A(scp¥) which then 
allows As to vary depending on the strength of the absorption at each spectral 
and spatial increment. 

The Curtis-Godson approximation has proven satisfactory for inhomogenities 
in rocket exhaust plumes, and the long cold atmospheric transmission paths for 
which the Intuitive Derivative method was developed do not exist in the usual 
heat transfer problem. But the program user should be aware of the character- 
istics of the model he intends to use. An example of H 2 O radiation through a 
long atmospheric path is illustrated in Fig. 3.4. Isothermal predictions of 
both cold and hot paths using Eq. 11a are compared with measurements in Fig. 3.4(a) 
and (b) showing good agreement with the measurements. The hot through cold 
measurement. Fig. 3.4(c), shows representative behavior of the two band models 
where the Curtis-Godson method typically overpredicts cold gas absorption. The 
solution obtained from the product of cold-gas transmissivity and hot-gas 
radiance underpredicts the absorption caused by line correlation effects as 
indicated in Fig. 3.4(d). 
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Gas Path 

Temperature 

Pressure 

Mole 

Fraction 

Call 

(m) 

(K) 

(atm) 

H 2 O 

Nz 

Cold 

100 

293 

0.07 

0.014 

0.986 

Hot 

0.6 

1202 

0.1 

0.5 

0.5 



30 00 32 00 3<i00 3600 38 00 <«0 00 42 00 

a)(aa"^} 


(a) Cold cell measurement and 
Isothermal band model pre- 
diction. 


(b) Hot ceil measurement and 

isothermal band model pre- 
diction. 


(c) Hot through cold measurement, 

Intuitative Derivative ( ) , 

and Curtis-Godson ( ) * 


(d) Hot cell prediction (upper) , 
product of isothermal pre- 
dictions of cold cell trans- 
missivity and hot cell radiance 
(middle) and Intuitative 
Derivative prediction (lower) . 


Figure 3.4 Comparison of Measurements and Band Model Predictions for 

HjO using an Exponential-Tailed-Inverse Intensity Distribution 
(Ref. 11). 
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Radiation Models in the GASRAD Program 

The developments leading to the GASRAD program have used the expo- 
nential line strength band model represented by Eqs. 10 and 12 with the 
modified Curtis-Godson approximation for inhomogeneous gas effects. The 
code is easily modified to use the exponentially-tailed-inverse distribu- 
tion (Eq. 11) if it is desired. A modified subroutine is available using 
the Intuitive Derivative method for an exponentially-tailed-inverse 
distribution so comparitive studies of HaO radiation results can be made. 

Choice of the band model distribution is not significant in most 
heat transfer applications where Lorentz broadening is likely to dominate. 
However, the increasing emphasis on weak lines as the distribution function 
proceeds from Model A to D produces a significant change in the Doppler curve- 
of-growth as mentioned earlier (see Fig. 3.26). 

The program is coded for five radiating species (HaO, CO 2 , CO, HCl , and 
HF) , and provision exist in the code for adding a sixth species. There is also 
provision for Including carbon particles (soot) with the assumption that there 
is neglegible scattering so that they can be treated as a gas in the linear 
limit. 

3.2 BAND MODEL DATA 

Band model data are obtained by both analytical and empirical methods. 
Generally the line structure of diatomic molecules is simple and well defined, 
so analytical predictions are effective. As molecules become more complex, 
analytical descriptions become more difficult and empirical methods are more 
reliable. Much of the data available are in a processes of evaluation, so 
the most recent generally accepted tabulations are usually the best. 
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A reasonably complete set of data, compiled for NASA/MSFC, was summarized 
in Ref. 10. It included analytical predictions for many diatomic gases and CO 2 , 
and empirical values for high temperature HgO. Young (Refs. 11 and 12) made a 
significant improvement in the data by merging the high temperature H 2 O data 
of Ref. 10 with the low temperature data developed for atmospheric transmission 
predictions (Ref. 13). He also proposed improvements in COj data (Refs. 12 and 
14), but comparison with measurements (Refs. 15 and 16) indicate the CO 2 data 
may be incorrect at low pressures. A similar result was found by Lindquist et 
al. (Ref. 17), but recent work by Bernstein (Ref. 18) appears to offer improve- 
ments in the CO 2 data. The new Standardized Infrared Radiation Model (SIRfUi) 

Code (Ref. 19) will include a data bank of band model parameters which is expected 
to serve as a resource for the best current data. 

The GASRAD code uses formated card- image input for the band model param- 
eters, so the parameters can easily be changed. It is planned to furnish the 
program with data sets corresponding to data in Ref. 10 extended by data from 
Refs. 20, 21, 22, and 23. Data sets for HgO and CO 2 combining Young's data 
with the data developed for NASA/MSFC are also furnished, and these are recom- 
mended in preference to the other data sets. Much of the data filled in to give 
a more extensive spectral interval and temperature range is not considered high 
quality, but it is useful for heat transfer predictions. Sources used in pre- 
paring the data sets are summarized in Table 3.3. 

Absorption coefficient data for carbon particles were taken from Ref. 24 
(also available in Ref. 10) for 2500 to 10000 cm”^ and extrapolated down to 
1000 cm"^ as indicated in Ref. 25. The data are approximated in the program 
as polynomials in wavenumber using the coefficients in Table 3.4. 
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Table 3.3 

SOURCES OF BAND MODEL PARAMETER DATA SETS 


Element 
Name , 

Parameter 

Wavenumber 
Range (cm"^) 

Temperature 
Range (k) 

Source 

H20GDC 

k 

50- 

9300 

300- 

3000 

Tabulation of Ref. 10 converted 
to gas temperature ref, rather 
than STP. 



9325- 

11000 

300- 

3000 

Tabulation of Ref. 20 converted 
to gas temperature ref. rather 
than STP. 


1/d 

1150- 

7500 

600- 

3000 

Tabulation of Ref. 10 



50- 

11000 

50- 

1125 

7525- 

11000 

300 

600- 

3000 

600- 

3000 

Sinusoidal approximation of 
Ref. 21: 

^n(l/d)=0.7941 sin(0.0036w-8.043) 

D(T) = -2.295 + 3.004xl0-^T 
-3.66xlO-’T2 

H20CMB 

k & 1/d 

2500- 

4500 


Tabulation of Ref. 11 interpo- 
lated to match temperatures of 
H20GDC 



50- 

2475 

4525- 

11000 

300- 

3000 

300- 

3000 

Same as H20GDC 

C02GDC 

k & 1/d 

500- 

875 

i 


Tabulated F of Ref. 10 at 
Aio=5cm“^ were used for 
Aw=25cm"^ without averaging, 
and referenced to gas temp, 
instead of STP, Value for 2400k 
used also at 3000 k. Values of 
1/d approximated frcmi values 
in Ref. 22. 



900- 

1875 


Filled with zeros for continuity. 

1 


1900- 

2375 

300- 

3000 

TabulatejJ values of Ref. 10 at 
Aw=10cm" were interpolated for 
Aco=2^cm”’ without averaging. 

All k values converted to gas 
temperature ref. 
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Table 3.3 (cent) 


SOURCES OF BAND MODEL PARAMETER DATA SETS 


Element 

Parameter 

Wavenumber 

Temperature 

Source 

Name 


Range ( cm" 

Range (k) 




2400- 

300- 

Filled with zeros for continuity 



2975 

3000 



3000- 

3000- 

Same source as 1900-2375 cm”^. 



3750 

3000 


C02CMB 

k & 1/d 

2000- 

300- 

Tabulations of Refs. 12 and 14 



2400 

3000 

at 5cm" ^ intervals degraded to 



3100- 

300- 

Aw=25cm'^ as indicated in Ref. 



3775 

3000 

12. Values interpolated to mate!' 
temperatures in C02GDC. 



500- 

300- 

Same as C02GDC. 



1975 

3000 




2425- 

300- 




3075 

3000 


COGDC 

i R 

1025- 

300- 

Tabulations of Ref. 20 changed 



2350 

5000 

to gas temp. ref. instead of STf 


1/d 

— 

i 

1 

Approximation from Ref. 23 and 
errata: 





1/d = 0.29(cm)[(H-F(T))/(l-F(T); 

F(T) = exp(-1561.5/T) 

(Tabulated values available in 





Ref. 10). 

HCLGDC 

k & 1/d 

1000- 


Tabulations from Ref. 10 change< 



3225 

3000 

to gas temp. ref. instead of STl 
and interpolated to give data a' 
25 cm"^ intervals. 

HFGDC 

I & 1/d 

1400- 

300- 

Tabulations from Ref. 10 changei 



4450 

3000 

to gas temp. ref. instead of STl 
and interpolated to give data a' 
25 cm"^ intervals. 
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Table 3.4 

CARBON ABSORPTION COEFFICIENT PLOYNOMIALS 


Temperature 

K 

A 

xlO'* 

B 

xlO* 

C 

xlO’^ 

D 

xl0~^ 

E 

xl0"^2 

300 

-.13463853 

.38710213 

-.47055911 

.35084341 

-.79087507 

600 

-.19909966 

.42759743 

-.50848071 

.37144485 

-.83813311 

1200 

-.31886445 

.48804827 

-.52853813 

.36589761 

-.80385920 

1700 

-.38870225' 

.53648682 

-.50955662 

.33136839 

-.69866129 

2000 

-.60273281' 

.86118469 

-1.0340689 

.64889548 

-1.3474213 

2300 

-89695742 

1.3577194 

-1.8257443 

1.1016709 

-2.2102493 . 

2600 

-.41368281 

1 .8460052 

-2.9400371 
1 

1.9221795 

-4.1482373 

1 


~ A + Bw + Cw^ + Do)® +£(*)** 

= Absorption coefficient, cm^/gn 
w = Wavenumber, cm“^ 
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3.3 SPECTRAL AVERAGING 

Increasing the spectral integration interval beyond the 5 to 25 cm'’^ 
for which band model data are usually available cannot be recoirmended on any 
theoretical basis for inhomogeneous gases. However, experience in heat transfer 
prediction from rocket exhaust plumes has indicated that larger intervals can be 
used to give satisfactory results for spectrally integrated radiation with signifi- 
cant savings in computation time. Since most computing time in the radiation (as 
distinguished from flowfield interpolation) portion of the program is spent in the 
wavenumber integration loop, any factor increase in the wavenumber interval results 
in about the same factor decrease in the computation time. So increasing the wave- 
number interval from 25 to 100 cm~^ will reduce the radiation integration time by 
a factor of about 4. Applications have been made of 400 cm"* intervals on HgO 
plumes and 100 cm"* intervals on H 2 O/CO 2 /CO plumes, but it is necessary to verify 
the procedures on each problem to assure reasonable results. This section briefly 
surveys methods of spectral averaging, and presents some of the results which 
have been achieved. 

Selection of averaging methods is influenced by the conditions of the applica- 
tion. If the linear limit of the curve-of-growth predominates, an effective value 
of k is appropriate, and the value of 1/d is only significant if it causes an in- 
correct indication of the location on the curve-of-growth. In the case where the 
square-root limit predominates, an effective value of /F/d is desired, so average: 
based /F^and / 1/d appear to be appropriate. Representative expressions for use ii 
the linear (L) and sqiare-root (S) limits will be used here with the Planck functioi 
included to weight the coefficients based on the blackbody distribution. The resuV 
ing expressions are 
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f 


k|^(u,Ato»T) 


^^{ujAwjT) 

l/dgCwjAo),!) 


m 


m 


E T)k(o)' Aco'.T)/ I 

n=l " " / n=l " 


m 


m 


[ S N°(o);,J) /F(u);,.Aa);T)/ E N"K, T)]S 
n=l / n=l " 


(18) 

(19) 


c E N^(w' T) /k(u)* A«;t)/c1(w' A0)',T) / S N°«»T)/l<K,Aa)‘ ,T)]' 
n=l " " / n=1 " " 

( 20 ) 


where m = Aw/ Aw' 

and W|J| = (w-Aw/2+Aw'/2) + (n-l)Aw' 

If Eqs. 19 and 20 are used to represent the square root limit, the desired 
result Is achieved, but if ¥|_ and 1/ds are used, then k/d becomes 

m , m _ 

(k,/d„)^ = (k-/d„)^ 2 [( E N®k)^/ E ( 21 ) 

^ ^ ^ n=l n=l 

with the two representations of the square root parameter differing by the factor 
in brackets. A similar deviation from the desired average occurs in using Fg 
rather than for predictions in the linear limit. In either case, if k does 
not change significantly in the interval. Aw, then either method of averaging k 
will give essentially the same result. As the interval increases so that there 
is a significant difference between kj^ and k^, it is likely that wavenumber interval 
will be approaching the limit of acceptable accuracy. However, limited experience 
with heat transfer predictions for relatively large rocket plumes has indicated ¥5 
gives slightly better results than ¥|_. Results for typical 1 ines-of-sight using 
and l/d5 are given in Table 3.5. These indicate that predictions for H2O/CO2 
plumes become unreliable above 100 cm"^ while the decrease in importance of the 
gaseous radiation when soot is present allows very large intervals. The two 
examples of water vapor plumes indicate Aw up to 400 cm"^ may be satisfactory. 

This has also been shown in results for the Space Shuttle Main Engine plumes 
(Ref. 26) listed in Table 3.6. 
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Table 3.5 

RESULTS OF USING LARGE SPECTRAL 
INTEGRATION INTERVALS IN TYPICAL APPLICATIONS (Ref. 10) 


Line-of-sight 
(Engine type 
or stage ) 

Gas 

Constituent 

Mole 

Fract. 

Aw 

cm"^ 

Predicted 

Radiance 

watts/ 

cm^-sr 

Percent 

Change 

from 

Aw = 25 cm”^ 

Rocketdyne 

HO 

0.359 

25 

4.978 


H-1 

CO 2 

0.167 

100 

5.008 

0.6 

(No carbon) 

CO 

0.335 

400 

5.760 

16 


H2 

0.139 

800 

5.970 

20 




1600 

5.926 

12 




8000 

5.232 

5 

(With 1 per- 

H 2 O 

0.359 

25 

18.612 


cent mass 

CO 2 

0.167 

100 

18.599 

0.1 

fraction of 

CO 

0.335 

400 

18.837 

1 

carbon) 

Carbon* 

0.020 

800 

19.013 

2 


H2 

0.119 

1600 

19.227 

3 




8000 

19.172 

3 

Saturn S-II 

H 2 O 

0.693 

25 

1.407 


Stage 

H2 

0.307 

100 

1.411 

0.3 




400 

1.449 

3 




800 

1.512 

7 




1600 

1.495 

6 




8000 

2.096 

49 

Rocketdyne J-2 

H 2 O 

0.693 

25 




H2 

0.307 

100 


0.1 




400 


2 




800 

0.1509 

4 




1600 

0.1622 

7 




8000 

0.1701 

13 


Note: All cases integrated from 1000 to 9000 cm“^ 

* Carbon is treated as a gas with a molecular weight of 12. 
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Subroutine ACDATA of the GASRAD program is coded to produce a table of 
Fg and l/d 3 (Eqs. 19 and 20) for a specified Aw which is an integral multiple 
of the Aw' used in the band model data input. Change of averaging to Fl (Eq. 
18) or other methods may be made by modifying a few lines of code. 


Table 3.6 

RESULTS OF LARGE SPECTRAL INTERVALS FOR A WATER VAPOR PLUME 



Predicted Incident Flux - Btu/ft^-sec 

Aw 

ET Base Center 

Nozzle Lip Normal 

Fin Trailing 

cm“^ 




to Eng. Centerline 

Edge Side 


SL 

20kft 

40kft 

SL 

20kft 

40kft 

SL 

20kft 

40kft 

25 

2.80 

1.84 

1.04 

10.5 

3.16 

1.41 

1.60 

0.90 

0.53 

400 

2.84 

1.88 

1.06 

10.7 

3.26 

1.45 

1.63 

0.93 

0.54 


The usual method of integration for band models is to evaluate radiance 
(for a single line-of-sight) over a spectral interval using 


Wg m 


= I 2 N° (t „ , - T, ^ „)Aw 
^ ,, w,n' w,n-l w.n' 

(JU]^ 11“’ I 


(21a) 


where Aw is a small interval and N° is evaluated at the center of each Aw 
interval from Wj to Wj. When large values of Aw are used, the radiance may not 
be represented well by the central value, so the radiation code uses a table 
of integrated values of the' Planck function to evaluate Eq. 21a as 

w+Aw 


= 
’m 


w^ m 


^2 (J0T4i(O CO 

= 2 Z {/ N° dw-/ N'» dw}(T •, - t ) 
0 co,n 0 w,n w>” ' w»'* 


(21b) 


Since linear interpolation is used in the table of integrated values, a spectral 
plot by the program of a continuum spectra using small Aw will show a stepwise 
appearance rather than a smooth curve. This can be avoided by replacing the 
call to the PLANCK subroutine (and the ENUl, ENU2 definitions) by the following 
code: 
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SE l./ENU 

SEE = SE * EXP(0. 47933 * ENU/TW) 

EYE = (1.1905E - 12/(SEE ** 3 - SE ** 3)) * DENU 
The use of degraded band model parameters is recommended when a large 
number of predictions are required and spectral data are not necessary, but 
results should be checked periodically to assure the desired accuracy is being 
achieved. Where only a few predictions are required, it is unlikely that the 
checking overhead will be justified by the computation time saving. 

3.4 SPATIAL INTEGRATION 

The transfer equation (Eq. 14) is integrated numerically using the geometry 
and numerical procedure indicated below. The spectral interval is set at the band 
model parameter resolution or a larger interval with averaged parameters as 
described in Section 3.3. This section will discuss considerations in the 
selection and use of geometric intervals As, A<f>, and A8. 



Selection of As should be made considering the features of the plume, the 
program run time limitations, and the memory size of the computer. The flowfield 
interpolation routines determine properties along each line-of-sight at input 
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specified As intervals. The value of As should be small enough to properly 
represent significant features in the plume, but large enough that the number 
of points on a line-of-sight will be less than the dimension parameter LPTMAX 
(which is normally set at 300). Once As is selected, values of A0 and sin 9 A<t> 
are normally assigned to be consistent with As. Equal spatial increments can 
be defined for a particular value of s by 

As = s sin0 AiJ) = s A6 - (23) 

and the program input provides specification of sin0 A<J) through the input 
variable TANGLE. TANGLE defines an upper limit A(j>' used in detemining 
a A4) which provides an integral number of increments in the range to (f)f. 

The procedure used is 

At})' = TANGLE/sin0 (24) 

Atj) = (4>i-«j)f)/IFIX((<})^-(|)f)/Att)'+l) (25) 

In a typical problem for the Space Shuttle Main Engine (SSME) plume 
radiation to the base, values of As = 6 inches and At})= si 00 Acf) = 2° have been 
used. These result in equal resolution in all three dimensions at s = 172 inches. 

As distance between the plume and the point-of-interest varies, Eq. 23 can be 
used to estimate consistent intervals. For example, to maintain the resolution on 
on the SSME plume for a point 1000 inches from the base, A0 and sin0 A*}* would 
have to be reduced to 0.34''. 

After defining property variations along lines-of-sight using the 
resolution defined by As, A0, and sinO A<p (or A^), the radiation subroutines 
(SLGMLT and YNGH20) may increase the As integration interval based on tempera- 
ture variations along each line-of-sight. The input variable TDIFF specifies 
the temperature change permitted in a line-of-sight segment which is to be 
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considered honogeneous. In the program operation, points are examined at 
the As intervals along a line-of-sight segment until a temperature, T(L), 
exceeds the initial temperature on the segment, T(l), by an amount greater 
than TDIFF, and if iT(L)-T(l) |>1.1 TDIFF, the program retreats to L=L-1. 

Then the properties over the interval (points 1 thru L) are averaged, and 
it is considered to be a homogeneous region with a length of LAs. 

This allows an increase in integration speed while maintaining a small 
resolution to define plume features. In the usual plume problem any signifi- 
cant change in plume peoperties is accompanied by a change in temperature, so 
the procedure used is quite effective. However, caution should be exercised 
if the gas being analyzed has large changes in pressure or concentration which 
are not accompanied by a significant temperature change. 

Because of the characteristics of the derivative band model procedure 
(Section 3.1) used in Subroutine YNGH20, the use of TDIFF will not produce the 
time saving that will be realized for the Curtis-Godson method used in Subroutine 
SLGMLT. But time savings will be achieved since the program will select a As' 
interval at each spectral position consistent with the absorption at that position. 
This procedure is control ed by the input variable IRAD which selects As' by a 
procedure which gives 

As' (w,s)-l/[IRAD * l<(w,s)c^* (s)p(s)] (26) 

for large values of the term in brackets and As' - As for small values. So when 
the product (Fcp) is small, a large spatial integration interval will be used. 

Very limited experience with the derivative method indicates values for IRAD > 100 
should be considered, but each application should be evaluated based on the pro- 
perty variations and accuracy desired. 
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4 GEOMETRY MODEL 

Discussion of the geometry model components will include: (1) the 
geometry of the receiving surface, shading surfaces, and plume; (2) the 
procedures used for defining line-of-sight limits as set by input limits, 
plume limits, and shading surfaces; and (3) the options available in 
describing the plume gas properties. 

4.1 COMPONENT DESCRIPTIONS 

The problem geometry is stated in terms of three components: the point- 

of-tnterest (POI), the plume, and the shading (or blocking) surfaces. These 
components are located in either the central coordinate system (CCS) or in an 
intermediate coordinate system (ICS) which is in a specified orientation in 
the CCS. All coordinate systems are right-hand as illustrated for a sample 
problem in Fig. 4.1. The CCS is defined in any orientation convenient for 
the problem, and the origin and orientation of the other systems are specified 
by location coordinates and three rotations. The sequence of rotations were 
selected to match the RAVFAC program (Ref. 27) so identical surface input could 
be used in both programs. The rotation sense and order are shown in Fig. 4.2 
with some of the nomenclature used in defining the transformations. 

The point-of interest (POI) is the first system input, and it is followed 
by the plume origin(s). In each case, a location and three rotations are 
specified along with an ICS number if applicable. The final position of the 
POI axes are represented by U, V, W with W representing the surface normal 
and the elevation and azimuth angles defined as indicated in Fig. 4.1. The 
final position of a plume system is typically X|^, Y^, Z|^ with representing 
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Figure 4.1 Overall Geometry Arrangement 
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DIRECTION COSINES OF X. , Y. , Z. RELATIVE TO X . Y , Z 

ODD dad 


Xb to Xa ciii *= cosx cosi{> 

Yq &i2 =-sinx cosij; 

Za ai3 = sinx 

Yb to Xq &2i - sinif; cosoj - cosx sim{j sinw 

Yg d22 = cosx cosco - s1nx sinij) sinw 

Za d23 = cosij; sinw 

Zb to Xg 331 =-sinx sinw - cosx simf; cosw 

Yg d 32 ='COsx sinw + sinx simii cosw 

Zg d 33 = cosip cosw 


TRANSFORMATION MATRIX AND ITS INVERSE 


ail 

®12 

®13 

''ab = ''ba = 

^11 

321 

^31 

»21 

322 

^23 

®12 

®22 

®32 

®31 

^32 

®33 


®13 

®23 

®33 


Fig. 4.2 Definition of Rotation Angles, Direction Cosines, 
and Rotation Matricies 
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the plume axis. Plumes may be described as one or more identical axi symmetric 
plumes with different origins or as a single three-dimensional plume. Details 
of the plume descriptions will be presented later. 

The blocking surface descriptions following the POI and plume input require 
additional data to describe the surface type, shape, and limits as shown in 
Fig. 4.3. In general, the variable a specifies the surface offset or size, while 
g and Y specify the two dimensional limits. 

Each coordinate system (POI, plume, or surface) is described in the program 
by a transformation matrix and translation vector which are used to locate a 
point in any of the systems. The final matrix and vector defined in subroutine 
BLOKIN is a result of several steps. First, the transformation matrix .elements 
for are prepared as though the system(s) were defined in the central coordinate 
system (c), and the location vector for s in c is defined as 

^sc “ ^xs^c ^ ^ys’ic ^ ^zs^c* 

Next, the matrix and vector for each intermediate coordinate system (i) is 
defined, and matrices and vectors of systems (s) referenced to (i) are modified 

by 


and 


A 

sc 


'ic 


A, 


sc 



^xs 


f^xi 

«sc = 

RyS 

Pi 

+ 

> 

PiJ 


(28) 

(29) 


After all surfaces are referenced to the CCS, a transformation matrix 
(T(I)) is prepared for each surface which will be used to transform directions 
in the CCS to a SCS 

T(I) = (30) 

where I = 1 designates the POI, 

2 < I £ (1+N) designates the N plumes, and 

I > (1+N) designates the blocking surfaces. 
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Figure 4.3 Blocking Surface Geometry 
39 


< 360°, 


90° 


270° 



F=? ^ rs/1 T" E c: h-i I M O- 


RTR 014-9 



%ax 



Limi tations: 

0 < a < 90, 

^max ^ ^in» 

-360 _< Ymin Y^gx £ 
and Ymax 1 (Ymin +360°) 



Limitations: 
a > 0, 


® £ ^min ^max 5. ^^0 
-3601 Y„,„ < < 

a"<* fmK 1 ('frain ''Oeo") 


Surface Type 6 Sphere 

Figure 4.3 (continued) 'Slocking Surface Geometry 


K/1 tie: o I— f I rvj o- 


RTR 014-9 



Limitations: 
a > 0 

^ 1 ^min ^max 
-360 1 Ymin ^n]ax 1 


Useful relations: 

y2 . s 4nf7 

[dX/dZly^o = 2a/X 


Figure 4.3 (continued) Blocking Surface Geometry 
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A translation vector R(I) is also prepared giving the location of the POI 
in all other systems {I > 1) 



4.2 LINE-QF-SIGHT OPERATIONS 

Each line-of-sight (LOS) is defined based on the orientation of the 
point-of- interest (POI) as represented by the transformation matrix T(l), 
defined earlier, and the direction in the POI system obtained by incrementing 
the hemispherical angles and 0 (Fig. 4.1) between the input limits of 
integration (4>i, <|>f, 0i, 9f). Operations to produce the LOS are performed in 
sequence at the center of each angular element beginning at 0 = 6-j + A0/2 and 
(}) = 4)f + A(|)/2. Integration limits for distance (s) along the LOS are adjusted 
as required to satisfy blocking surfaces and plume limits, then plume properties 
are obtained by interpolation at the center of As intervals beginning at 
s=Sq+As/2. 

The direction of the LOS relative to the CCS is defined by the direction 
consines a, b, and g 

r T fsine coscH 

b = T(l) J sine sin^ , (32) 

_gj ^ Lcos0 

and the components of the LOS vector in each surface coordinate system are 
defined by 



where T'(I)is the inverse of the transfomration matrix for surface I and S^ 
represents the upper limit of the LOS. Since the vector ^(I) (Eq. 31) repre- 
sents the position of the POI in the surface system, the equation of the LOS 
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in the surface system (Xg,Yg,Zg) is 
Xs-Rl Ys-R2 Zs-R3 
Cl C2 C3 


(34) 


where 6 represents a fraction of the magnitude of the LOS upper integration 
limit So a solution is only applicable if 6 lies in the range 0 to 1. 

For the plane surfaces (Types 1, 2, and 3 in Fig. 4.3), the defining 
equation in the surface system is Zg = a which combined with Eq. 34 yields 
6 = (a-R3)/C3 (35) 

If 0 < 6 < l,the LOS may be blocked by the surface, so the position 

Xs = R1 + 6C1 (36) 

Yg = R2 + 6C2 (37) 

is computed and checked against the surface limits defined by 3 and y (Fig. 4.3). 
The conic shapes (Types 4, 5, 6, and 7 in Fig. 4.3) can be described by 
X| + Y| = A + B Zs + C Z| (38) 

where the coefficients for the surfaces represented are 

A_ -L. C 

0 0 

0 0 tan^a 

a" 0 -1 

0 4a 0 

Solution of Eqs. 34 and 38 yields 

AA 6^ + BB 6 + CC = 0 (39) 

where 

AA = (Cl^ + C2^ + C C3) 

BB = (2 R1 Cl + 2 R2 C2 - B C3 - 2 C R3 C3) 

CC = (Rl^ + R2^ - A - B R3 - C R3^) 

If either solution of Eq. 39 lies between 0 and 1, the position Is computed 

(Eqs. 36 and 37) and tested against the surface limits defined by 3 and y 


Cylinder 

Cone 

Sphere 

Paraboloid 
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(Fig. 4.3). If both solutions are between 0 and 1 and lie within the surface 
limits, then the smaller solution is chosen. 

As each surface is tested in sequence, the upper limit of integration 
is reduced to correspond to the resulting 6 rf it is smaller than any previous 
solution. The final solution then represents the distance to the nearest 
blocking surface, and the LOS upper limit of integration is redefined by S^j, 6Sj^. 

After checking for shading, the LOS is tested, to determine if it intersects 
a conical boundary around the plume. The boundary is defined during surface 
input using the same parameters as a blocking surface, but the parameters are 
defined by the program using the plume origin and the intercept (AP) and slope 
(BP) of the boundary at = 0. Solution for the intercepts of the LOS. and the 
plume boundaries is carried out in Subroutine SLIMIT using the same approach as 
for the blocking surfaces, but additional tests are required to set integration 
limits as a function of the LOS orientation. If the LOS lower limit computed in 
SLIMIT is greater than the input lower limit Sq, the lower limit for the LOS 
is adjusted. Similarly, if the final intercept is less than the upper limit 
(set by either input or blocking), the upper limit for the LOS is reduced. 

In all operations adjusting the LOS limits, the values used are multiples 
of As which fall just outside of the limit solutions. Once the limits are 
adjusted, the LOS distance is initialized at s = so + As/2 and incremented in 
As intervals until s > sm* At each value of s an interpolation is made in the 
plume property tables (described later) to define the temperature, pressure, and 
mole fractions of species identified by the input. Even though the portion of 
the LOS (Sq to Sjn) is within the cone bounding the plume, a point may not be 
within the plume property table since the cone does not usually fit the plume 
at all points. Points outside the plume are assigned a zero temperature, and 
after the LOS is complete, the temperatures are checked to determine if they 
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are non-zero and within the temperature range normally used in the band model 
parameter tables (300K to 3000K). If all temperatures are zero, the LOS is 
identified as having missed the gas, and if any points are found outside the 
300K to 3000K range, a warning message is printed at the termination of all 
LOS processing. 

The sequence used in 1 ine-of-sight testing was chosen to indicate first 
blocking and second missing the plume. This requires more computation than 
the reversed procedure since processing of each LOS missing the plume bounding 
cone could be immediately terminated. But the procedure used provides a more 
complete picture of the geometry to assure proper orientation of the POI and 
surrounding structure. 

4.3 PLUME DESCRIPTIONS 


Topics in the plume descriptions include the two basic fonns, axi symmetric 
and three-dimensional, along with methods of interpolation and recognition of 
discontinuities caused by shocks. Dimensional considerations are also discussed 
including methods of stacking plume sections. 

Forms of Gas Property Arrays 

Plumes use a cylindrical coordinate system as illustrated below with pro- 
perties specified as a function of Z, n* R in general. But for the axisymmetric 



Fig. 4.4 Basic Plume Coordinates 
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case, a single value of n is used and taken to apply for all n, and for 3-D 
plumes, optional parameters are available to modify the description to take 
advantage of engine cluster symmetry. 

The axisymmetric plume description is used for plumes which are assumed 
to be non-interfering. Up to six separate plume origins and orientations can 
be specified for a single set of axisymmetric plume data to simulate multiengine 
arrangements. (The limit of 6 plumes was chosen to fit existing input control 
flag limits and could easily be expanded.) Descriptions of the Space Shuttle 
low-altitude plumes (Ref. 26) used this procedure with three plume origins 
located at the three main engine exits even though plume interference occurs 
between plumes in this arrangement because of the close engine spacing. 

Operation of the program choses the closest plume for evaluation of gas properties, 
so adjacent plumes are effectively terminated at a surface equidistant from the 
plume axes. 

The form of 3-D flowfield input is governed by three parameters intended 
to facilitate the property description using symmetry in engine clusters. Brief 
descriptions of the parameters identified by the input names are: 

CTCD Represents the center-to-center distance between an engine 

and the axis of the engine cluster. Actually it is the 
distance between the plume cluster axis and an axis used 
for property input which is displaced along the axis. 

HANG Represents the half-angle of a symmetrical sector of a plume 

cluster. The sector is assumed to have symmetry about the 

X axis, 
w 

ISYMSC A flag to specify assumptions of symmetry about a plane normal 

to the X^ axis at a distance of CTCD/2 from the axis. This 
is used if the center and outboard plumes in a cluster are 
assumed to be mirror images within the limits defined by HANG. 

A few examples illustrating the use of these parameters are shown in Fig. 4.5. As 

noted in the descriptions above for cases of CTCD > 0, the axis used for property 

definition { 1 ^) is defined parallel to 1 ^. 
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(c) Symmetry about Yy^, A, 0-180 >0 45® 1 

and B. For a 4 engine 
cluster. ,, 



' (d)Symmetric about X^, Y^, 0-180* >0 45* 0 

A, B, and CD. For a 6 
engine cluster. 

Fig. 4.5 Illustration of optional plume descriptions using the parameter CTCD, 
HANG, and ISYMSC. 


47 



F=? e: rs^ nr ^ o h-H i r^j o. 


RTR 014-9 


For engine clusters which are not symmetrical or in which the engines 
are canted, special consideration needs to be given to the 3-D plume descrip- 
tion. Use of the general case (Fig. 4.5a) requires a relatively large internal 
storage capacity to obtain satisfactory detail in the plume, but the more 
effective descriptions represented by CTCD > 0 (Fig. 4.5c and d) may not fit 
well for non-symmetric clusters or canted nozzles. In applications for the 
Space Shuttle Main Engine (SSME) vacuum plume (Ref. .26), a symmetric 3-eng1ne 
cluster was assumed which approximated the actual arrangement, and a flowfield 
transformation program (DTRANS, Ref. 28) was used to put the flowfield developed 
for canted nozzles into the form required with the flowfield property axis parallel 
to the cluster axis (Z^). In the SSME vacuum plume some characteristics of this 
three-dimensional solution became unacceptable downstream, so the 3-D portion 
was terminated at = 200 inches. The remainder of the cluster flowfield was 
approximated by an extension of the axi symmetric central core of the flow which 
was expected to contribute most significantly to the plume radiation. To accom- 
modate this circumstance, an option was included in the program to switch from 
a 3-D plume in the nearfield to an axisytraietric approximation in the farfield. 

Due to the large internal memory requirement for the plume, dimensioning 
deserves special attention. Each location in the plume at which properties are 
specified requires three memory locations for radius, temperature, and pressure 
plus the number of species mole fractions dimensioned for the problem (parameter 
NK, data NKO) which is usually less than the allowable number of species in the 
input flowfield (parameter NKTP). The parameters specifying the maximum number 
of flowfield locations are outlined below: 

Axis.ymmetric 3-D 

lyi positions MZ2 MZ3 

T) positions at each 1 MN3 (40) 

Radius positions at each n MR2 MR3 

t 
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Combining these parameters with the number of values stored at each point 
gives an approximation of the flowfield internal memory requirement: 

Axisymmetric (3 + NK)(MR2){M22) (41a) 

, 3-D (3 + NK)(MR3)(MN3)(MZ3) (41b) 

An early goal in the program was to maintain the flowfield memory requirement 
in the IIK to 12K word range so the program could be run within a 32K word 
memory limit. Values chosen for the multiconstituent program (NK*6) were 
MR2=MR3=30, MN3=13, MZ2=40, and MZ3=3. For the SSME predictions, NK was 
reduced to 2 (for HgO and Nj) and the Zw positions were increased to MZ2=78 and 
MZ3=6. Even with this increase, the severe limitation on MZ3 is apparent. 

Because of this limitation, the program was designed to handle the plume in 
sections to prevent a restriction in length by the internal memory size. 

Procedures for stacking plume sections will be automatically applied by the 
program to either a single axisyranetric or 3-D plume. The flexibility allowed in 
positioning multiple axisymmetric plumes was not compatable with the stacking 
procedure. The procedure starts by reading in the allowed (dimensioned) Zy/ posi- 
tions, interpolating for propertie.s on all lines-of-sight, and storing the partial 
solution on Unit 3. The last Zy/ position in memory is moved to the initial loca- 
tion and additional Z^^ positions are loaded. Partial lines-of-sight are read 
from Unit 3, extended, and written on Unit 4. This procedure is repeated with the 
partial lines-of-sight shifting between Unit 3 and 4 unitl the plume processing 
is complete. At that time the completed lines-of-sight are output on Unit 10 as 
they are when no piecing together is required. This process is relatively straight- 
forward for lines-of-sight in the +Zy, direction, but it requires some maneuvering 
when a line-of-sight is in the -Zw direction because the last section of the 
line-of-sight is completed first. 

The large amount of input/output required when the plume sections must be 

stacked has become quite significant in the current computer cost assignment 
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system at MSFC, Therefore, it is recommended that the dimension parameter 
product (Eq. 41) be made as large as currently allowable and the relative 
parameter values be adjusted to the particular plume requirements so the 
plume does not have to be sectioned. 

Interpolation Methods 

The interpolation procedure used to define gas properties in the plume 
will be described first for axi symmetric plumes. Then the additional inter- 
polation in n required by 3-D plumes will be discussed. Plume data are stored 
as an array at a set of values and arbitrary values of An illustration 
of a pair of adjacent Z-cuts in Fig. 4.6 indicates the physical arrangement 
including the optional shock location. 

The property search begins with a location in the plume (ZW, RW) of the 
point desired on the line-of-sight. The selection of the plume (in multiple 
plumes) and transformation of the location to the plume coordinate system is 
performed in Subroutine PLUMPT. The values of Z in the plume property table 
are searched sequentially* to define the two values Z(I1) and Z(I2) containing 
the interval in which ZW is located, and the ratio for linear interpolation is 
defined (refer to Fig. 4.6 for nomenclature) 

ZRATIO = (ZW-Z(I1))/(Z(I2)-Z(I1)) (42) 

Next, the plume boundary is defined at ZW by 

RWMAX = RMXl + (RMX2-RMX1) *ZRATI0, (43) 

If RW > RWMAX the interpolation is terminated, otherwise, values of RWZ(l) and 
RWZ(2) are defined to represent RW in radial interpolation at ZZ(Il) and ZZ(I2). 

* Test using a uniform binary search in Z and R did not indicate a measurable 
saving in overall computation time for a typical problem. 
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Plume Boundary 


11 

(Normally the 
nozzle exit) 


<_R(I2,NRP)=RMX2 
9 NRP=NR(I2)-9 


^(I1,NRP)=RMX1 


,"lr"NRP=NR(Il) 


-f A(ZW,RW) 


^R(I2,LSP)=RS2 
,6 LSP=LS(I2)=7 

,*^RWZ(2) 


^R(I1.LSP)=RS1 
i5 LSP«LS{I1)=6 



Fig. 4.6 Illustration of Axisymmetric Plume Data Arrangement 
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The code is currently written to interpolate parallel to the plume axis 

RWZ(l) = RWZ(2) = RW, (44) 

but comment cards are included for code which can be used to vary the interpolation 
so that it becomes parallel to the plume boundary at the boundary. The revised 
procedure which was used in earlier versions of the code is 

RDUB = RW/RWMAX 

RWZ(l) = RMXl *RDUB (45) 

RWZ(2) = RMX2 *RDUB 

Although this procedure appears attractive, it produced unrealistic results for 
some portions of the SSME (3-D) vacuum plume due to large changes in the boundary 
and properties in the n-direction. 

The radial values of Z(I1) (and Z(I2)) are searched sequentially for the 
intervals containing RWZ(l) (and RWZ(2)). Search limits are determined by the 
optional shock specification. If a significant shock is in the plume, two radial 
points very close together should be specified to define the properties on both 
sides of the shock. The radial index of the smaller radius point is indicated by 
LS(I) where I is the index of the Z-cut. The shock locations RSI and RS2 are define 
at consecutive Z-cuts as shown in Fig. 4.6. If both radii are >0, the shock option 
is invoked, and the shock position at ZW is computed 

RS = RSI + (RS2-RS1) *ZRATI0 (46) 

This radius is tested against RW to determine if the point is inside or out- 
side the shock surface. The outcomes of the shock handling, and the resulting 
radial search limits for RWZl and RWZ2 are summarized belnv?: 

Condition Flag Radial Search Index 

RSI = 0 or RS2 = 0 LSHOCK =0 N = 1, NR(I) 

RW > RS LSHOCK =1 N = [LS(I)+1]. NR(I) 

RW < RS LSHOCK = -1 N = 1, LS(I) 

When the radial interval is located containing RWZ(II) (where II is 1 or 2), 
linear interpolation is performed to determine the corresponding properties 
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(temperature, pressure, and species mole fractions). If the upper limit of 
the radial search Is reached before finding the interval containing RWZ(II), 
the properties of the point at the end of the search are used for RWZ(II). 

After properties are determined at RWZ(l) and RWZ(2), a linear interpolation in 
is performed for all properties (similar to the procedure in Eqs. 43 and 46). 

As an illustration, consider Point A in Fig. 4.6 as the line-of-sight 
point located at ZW, RW. First, Z-va1ues in the table are searched to identify 
Z(I1) and Z(I2) and the interpolation fraction ZRATIO (Eq. 42) is computed. 

The plume boundary points at these Z-cuts are identified, the plume boundary 
is defined (Eq. 43), and the point is found to be inside the plume (RW < RWMAX). 
Shock radii at Z(I1) and Z(I2) are >0, so the shock radius is computed at ZW 
and the point is determined to be inside the shock (RW < RS) which sets LSH0CK=-1 
The desired radial values for interpolation at each Z-cut are defined (RWZ(l) " 
RWZ(2) = RW), and a loop is entered to define properties at RWZ(l) and RWZ{2). 
Since the search is limited by the shock position, the search at Z(I1) cannot 
reach RW, so the properties at R(I1,6) are used for RWZ(l). Properties for 
RWZ(2) are determined by linear interpolation between R(I2,5) and R(I2,6). 

Finally properties at Point A are determined by linear interpolation in Z 
between RWZ(l) and RW2(2). 

Without shocks, 3-D interpolation is a relatively straightforward extension 
of the axisymmetric procedure, but with shocks, special procedures are added 
which complicate the description. The values XW, YW, ZW returned by Subroutine 
PLUMPT for the location of a point in the plume are converted to ZW, SI6, RW, 
where SIG represents the angle n in the plume property description (refer to 
Fi,g. 4.4). All interpolation is linear and the sequence of interpolation is 
radial (RW in R^), angular (SIG in n)^ and axial (ZW in 


53 



i=?^rs/i-rE:o»— I n^o. RTR 014-9 

As in the axisymmstric case, 2 ^ is searched first to find II and 12 
which satisfy Z(I1) £ ZW < Z(I2). Then the values of n are searched at each 
of the two Zw locations to define n{Jl) £ SI6 < n(J2), and the maximum radius 
(RSIGMX(II)) and shock radius (RSHOK(II)) are determined (at each Z) by linear 
interpolation in n between the values at n(Jl) and ri(J2). An interpolation in 
Z^, is then made to obtain the maximum radius (RWMAX) and shock radius (RS) at 
ZW. If shocks are not defined for all four sets of radial data involved in this 
process, the interpolation for properties follows the same sequence as summzrized 
above for RWMAX with linear interpolation in radius to match RW at each of the 
four n-cuts (2 at each Z-cut). 

Existence of shock data (LS(I,J) > 1) for all four ri"Cuts invokes the 
shock option which has a special procedure useful in flowfield descriptions for 
clustered engines. Typical examples of the two cases considered are shown in 
Fig. 4.7. In each case the procedure is designed to treat the properties as 
axi symmetric in the region between the nozzle centerline and the shock. This 
is accomplished by using radial interpolation on the n-cut with the largest 
shock radius. For the case in Fig. 4.7a the radial properties at n(01) would 
be used, while the case in Fig. 4.7b would use properties at ti(J2). For 
example. Point B in Fig. 4.7b lies between n(Jl) and ri(J2) and inside the shock 
surface, so its properties will be best represented by the properties on n(J2) 
at the same radius. As a result the program determines properties at Point B 
as though it were in an axisymmetric plume described by n(J2). 

This feature can be disabled by setting control flag ISHOCK = 1. The 
result in Fig. 4.7b would be a linear interpolation in n between properties 
at RSI (on the side of the shock closest to the nozzle) and those at a radius 
RW along ri(J2). 
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Fig. 4.7 Illustrations for the 3-D Flowfield Geometry With a 
Shock for a Typical 4-Engine Cluster. 
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The 3-D code is written to interpolate in parallel to the flowfield 
axis as described in Eq. 44 for axi symmetric plumes. However, code is 
included as comments which can be used to interpolate as in Eq. 45 for 
axisynmetric plumes. 
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A1 INPUT GUIDE 

The input guide describes aspects of the program set-up and mapping, 
the card input formats and order, and the unformated binary input. Nearly 
all of the binary files listed in Table Al.l are used for both input and 
output, so they cannot be clearly assigned to either category. However, the 
program is usually run with the gas property data (flowfield) input so it will 
be described as an input file although it can also be generated as an output 
file from card input. The other binary files (problem data and line-of-sight 
properties, radiance data, and spectral data for plotting) will be described 
as output since they are normally generated by the program. 

Al.l PROGRAM SET-UP AND MAPPING 

The job control cards currently used on the MSFC UNIVAC 1108 are illustrated 
with the Sample Problems (Appendix A3). A typical sequence involves: assignment 
of input tape(s) containing the program source code, object code, map element, 
band model data elements (card image), and flowfield property input (binary); 
Increase in the default size of the temporary program file; and assignment of 
file units to be used by the program (3, 4, 8, 10, and occasionally 11). The 
code and band model data elements are copied into the temporary program file 
while the flowfield is copied to Unit 8 (or 11 if it is to be modified as noted 
on Input Card 2), and the object code is executed if no updates are required. 

The normal program map element in Table Al.2 shows the segmentation and 
usual deck names. All deck and subroutine names are identical with the excep- 
tion of Subroutine SLG. Here, the normal radiation deck using the statistical 
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UNIT 

NUMBER SUBROUTINE 

3 FLOWAX 
FL0W3D 


Table Al.l 

EXTERNAL FILE UNITS FOR THE GASRAD PROGRAM 

FUNCTION 


4 FLOWAX 

FL0W3D 

SLG 

VIEW 

SPCOUT 
PLOT IT 

8 FFPREP 

FLOWIN 

FLOWAX 

FL0W3D 

10 FLOWIN 

FLOUT 

FLOWAX 

FL0W2D 


11 


SLG 

FFPREP 


Used as temporary storage when it is necessary to 
load the plume in sections. If the dimension param- 
eter is larger than the number of Z-cuts in the 
plumes, it is not required, but it is still referenced 
by rewind statements. 

Same as Unit 3 


Used to store and transfer line-of-sight radiance data 
from Subroutine SLG to Subroutine VIEW. 

Used to transfer spectral data from Subroutine SPCOUT 
to the plot program PLOTIT. 

Stores gas property data read from cards or Unit 11 
Provides flowfield title and constituent input 
Provides flowfield gas property data. 


Stores input data so it can be used for program restart 
with abbreviated input. 

Uses stored input data for problem output. 

Used to store line-of-sight gas property data for 
transfer to Subroutine SLG. The property data are 
stored after the problem specifications used in 
Subroutine FLOUT. 

Used to input line-of-sight gas property data. 

Input gas property data which is to be modified 
(scaled) and output on Unit 8 for use by the program. 
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Table A1.2 

NORMAL PROGRAM MAP ELEMENT 


q!MAP,S A,GASRAD 

M AP ..29R I SL 73R1 .01/03/80 19 15.9,:.2 2 i 4 , ) 


1 . 

LIB SYS$*MBFC$. 


2 • 

SEG mains 


3. 

‘i . 

IN 

SEG L0S6EN«, (MAINS ) 

MAIN 

5 . 

IN 

F LOW, SL IMIT, PLUMPT.B LOCK M.NOGO 

6* 

SEG FL0WIS«. (LOSGEN) 


7. 

IN 

FLOWIN .BLOKIN.FFPREP 

8. 

SEG FLOUTS#, ILOSGEN 1 

* 

9. 

IN 

FLOUT, SKIP 

1C. 

SEG FL0WAS«t (LOSGEN) 


11. 

IN 

FLOWAX 

12. 

SEG FLOW 3S#, (LOSGEN) 


13. 

IN 

FL0W3D 

14. 

. ... SEG RADS*. (MAINS) 


15, 

IN 

PAD 

16. 

SEG_,...ACD/YTS#, ( R ADS) 


'"“’'"'17."' 

’ ' ' 1 

ACDATA 

18. 

SEG SLGS#,(9ADS> 


Y"9 . 

IN 

SLGMLT .PLANCK 

20. 

SEG SPCOTS*. (RADS ) 


■'21." 

IN 

SPCOUT 

22. 

S.E .G„„ .V IE.K:.S «. « MAINS). 



... 
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band model with exponential line-strength distribution and the Curtis-Godson 
approximation for inhomogeneous gases has the deck name SLGMLT. If a prediction 
is desired for H 2 O using the Intuitative Derivative Method (Section 3.1) to 
account for inhomogeneous gas effects, the program should be mapped with the 
deck name YNGH20 used in place of SLGMLT, and Subroutine PLANCK can be omitted. 

All declariations are included in program subroutines through reference 
to procedure entries in the (UNIVAC 1100 Series EXEC 8) POP (procedure defini- 
tion processor) element COMMLT. Reference to the procedures in a subroutine 
causes them to be inserted in the source code. Use of the program on another 
computer system will require changes to this system and may require duplications 
of the declarations for insertions at appropriate locations. 

A1.2 CARD INPUT 

The card input formats are described in Table A1.3. Since the program 
was used primarily for the Space Shuttle Main Engine (SSME) plumes, the input 
is designed to utilize program data statements for SSME plume locations unless 
IFL0IN=1 (Card 2). However, if the program were to be used extensively for 
another engine arrangement, the data statements could be changed accordingly to 
simplify input. 

Card input requirements for the flowfield were omitted from the normal 
sequence and inserted at the end of Table A1.3 to simplify the input card 
sequence listed to that normally used with the flowfield input as a binary 
file. 

In earlier versions of the program, (Refs. 1-6) the band model parameters 
were included in data statements, which is convenient from an input standpoint, 
but it limits flexibility in evaluating performance with other input sets and 
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limits the ease with which the program dimensioning can be varied for appli- 
cations requiring less data. The input data elements summarized in Table 3.3 
are usually included in the program storage file in the proper input format 
specified by Cards 17 through 19 so they may be added to the runstream with 
system control statements (@ADD) as illustrated In the sample problems. The 
coding is set to read band model parameters data for the first NKA (parameter) 
constituents identified on the program constituent data list IDENTC(I) (where 
there are 2 six-character words per constituent). In normal form, NKA=6, so 
the program will attempt to read six data sets. If only one data set is to be 
loaded, as in a run for H 2 O only, the H 2 O data would be loaded followed by 5 
blank cards to indicate no data are loaded for the other 5 gases. 

A1.3 UNFORMATED BINARY INPUT 

Plume property data input on Unit 8 or 11 (Table Al.l) are normally 
considered as input, so the format is described here in Table A1.4. Other 
binary files are described in Appendix Section A2.2. For further clarifica- 
tion of the number of Z-cuts (Group 2 Record 1) allowed by the program the 
reader should refer to Section 4.3. 
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Table A1.3 
INPUT GUIDE 


CARD 1 (12A6) TITLE CARD 

1-72 Heading for the problem - Blank terminates input. 


CARD 2 
1-6 


7 

8 
9 


12 

13-24 

26 

27 


NOTE: 

29,30 

33 

36 

39 

42 

45 


(16,2013) CONTROL CARD 

KASE =0 Repeat previous case with scaling - See Note C at end of 
Card 2. 

>0 Case number identifying the flowfield on tape. 

<0 Flowfield input on cards and transferred to unit 8 with 
scaling if desired. Case number on tape will be {KASEj. 

IFLOTP =0 Flowfield not modified. 

>0 Flowfield input on unit 11, modified as indicated by scale 
factors (Card 8) and output on unit 8 for use in this problem. 
IFLOIN =0 Nozzles (flowfield origin) located using data statement in 
subroutine FLOWIN (SSME flight null). 

=1 Nozzle location input on Card 10, 

IFLOW <7 Number of identical axisymmetric plumes 

=7 Use previously generated line-of-sight data. See Note B. 

=8 Three-dimensional flowfield. 

=9 Three-dimensional flowfield followed by an axisymmetric flow 
flowfield after IZF Z-cuts. (See note for IZF Col. 50, 51) 

NOKON Number of species to be identified on Card 6. (The program is 
normal dimensioned for a maximum of 6.) 

Not used. 

NSPEK =0 Spectral output per unit wavelength- 
=1 Spectral output per unit wavenumber. 

=0 No spectral output. 

=1 Spectral radiance per steradian - short form 

~2 Spectral radiance per steradian - long form 

=3 Spectral radiance for region specified on Cards 3 and 4 - short. 

=4 Spectral radiance for region specified on Cards 3 and 4 - long. 

All spectral parameters other than radiance are for the last line-of-sight 
processed. 

NPLOT Number of spectral plots specified on Card 16* (Requires NSPEK > 0) 
NFLOW ==0 Short flowfield output with only centerline and boundary. 

=1 Complete flowfield output. 

NABS =0 Short form absorption coeff. output 
=1 Complete absorption coeff. output 
ISHOCK =0 Special 3-D shock procedure for clustered engines Refer to 

=1 Disable special 3-D shock procedure Section 4.3 

lUNITS =0 Flowfield and geometry in English units (inches, R, Ibf/ft^). 

=1 Flowfield and geometry in Metric units (cm, k, atm) 

ISCALF =0 No plume scaling. 

=1 Plume scale factor if KASE ji^ 0. Entire geometry scaled if 
KASE = 0. 


* The program is currently configured so that the plotting is done by PLOTIT 
which is compiled and executed separately but has all input passed from the 
radiation program on unit 3. 
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47,48 

50,51 

54 

60 


IZI =0 Read first Z-cut in input flowfield. 

>0 Index of first Z-cut to be used. 

IZF =0 Read entire flowfield. 

>0 Index of final Z-cut to be included, or last Z-cut in the 
three-dimensional flowfield if IFLOW = 9. 

IRAD Spatial integration control for Deck YNGH20 only. See Eq. 26. 
ISCALA =0 No scaling of band model parameters (k and 1/d). 


NOTES FOR USING PREVIOUSLY GENERATED LINE-OF-SIGHT DATA 

A. Input line-of-sight data on unit 10. 

B. If a change is desired in problem limits: 

1. Set IFLOW = 7 and KASE equal to the case number used in 
generating the line-of-sight data. 

2. Input Cards 3 thru 6 using geometric integration limits and 
intervals which will produce line-of-sight data available on 
Unit 10. For example A0 can be doubled but not halved. 

3. Omit Cards 7 and 9-14. Use remaining cards as required by 
options chosen. 

C. If only scaling is desired with no change in limits; 

1. Set KASE = 0 

2. Omit Cards 3-7 and 9-14. Remaining cards are used as 
required by options chosen. 

CARD 3 (8E10.6) ELEVATION ANGLE INTERVAL - DEGREES 

I- 10 THETAI Initial value (6,). 

II- 20 THETAF Final value (Of). 

21-30 DTHETA Integration Interval (A9). 

CARD 4 (8E10.6) AZIMUTH ANGLE INTERVAL - DEGREES 

I- 10 PHII Initial value (<j>i). 

II- 20 PHIF Final value 

21-30 DPHI Integration interval (A4>). 

31-40 TANGLE Used to increase A(j> at small 9. A(|) = TANGLE/sinG will be 
used if it is larger than DPHI. 

CARD 5 (8E10.6) WAVENUMBER INTERVAL - CM"^ 

I- 10 ENUI Lower spectral limit - w-j. 

II- 20 ENUF Upper spectral limit - wf. 

21-30 DENU Integration interval - Aw 


CARD 6 
1-10 
11-20 
21-30 
31-40 


41-50 

51-60 


(8E10.6) LINE-OF-SIGHT INTERVAL - UNITS PER COL 42 ON CARD 2 
SZERO Distance from point-of-interest to initial LOS point 
SMAX Upper limit of LOS from the point of interest. 

DS Integration increment on the LOS. 

TDIFF Temperature step size used to increase LOS integration 

increment. Assumes gas with a temperature change of less 
than TDIFF is isothermal. 

AP These values for the plume limiting cone intercept (AP) 

BP and slope (BP) will overide plume input value from binary or 

card input (Card 2F) if either AP or BP are >0 on this card. 
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CARD 7 

(126A) 

1-6 

KON(l) 

7-12 

K0N(2) 

13-18 

K0N(3) 


CARD 8 

(9E8.5) 

1-8 

RSCAL 

9-16 

TSCAL 

17-24 

PSCAL 

25-32 

FSCAL(l) 

33-40 

FSCAL(2) 

65-72 

FSCAL(6) 


CONSTITUENT IDENTIFICATION 

First Constituent - first word 

First constituent - second word 

Second constituent - first word 

Continue with a list of NOKON (Card 2) two word names 

in the form specified in the IDENTC data list: 

H201(G). C102(G), ClOl(G), HCL, HF. C1(S), N2(6), H2(G) , 01(G), 

FLOWFIELD/GEOMETRY SCALE FACTORS - USE ONLY IF ISCALF = 1 
Length scale. Notes: 1. If KASE f 0 (Card 2) the flowfield will 

be scaled, but if KASE = 0, the entire 
problem will be scaled. 

Temperature scale. 2. Zero entries will be set to 1. 

Pressure scale. 

Mole fraction scale for first constituent on Card 7. 

Mole fraction scale for second constituent. 


Mole fraction scale for sixth constituent 


IF KASE < 0 INSERT CARD FLOWFIELD DATA HERE. 
SEE INPUT REQUIREMENTS AFTER CARD 19. 


CARD 9 

9,10 

21-30 

31-40 

41-50 

51-60 

61-70 

71-80 


(8X, 12, lOX, 6F10.3) POINT-OF- INTEREST LOCATION AND ROTATION 
ICS(l) Intermediate coordinate system number. 0 if none is used. 

RX (1) 1 

RY (1) [ Coordinates of the point in the ref. or inter, coord, sys. 

RZ (1) j 
CHS(l) X] 

PSS(l) ij) } Rotations of the axes relative to the ref. or inter, coord. 
OMS(l) w{ sys. 


CARD(S) 

10 (8X, 

12, lOX, 6F10.3) 

9,10 

ISC(I) 

Inter, coord. 

21-30 

RX (I) 1 

1 

31-40 

RY I 

[ Coordinates 

41-50 

RZ (ij i 

1 

51-60 

CHS(I) 

X 1 

61-70 

PSS(I) 

} Rotations 

71-80 

OMS(I) 

0)1 


PLUME ORIGIN (S) LOCATION AND ROTATION 
system. 


Notes: 1) USED ONLY IF IFLOIN = 1 

2) CARD(S) 10 I FLOW IF IFLOW < 7 

3) CARD(S) 10 = 1 IF IFLOW > 7 


CARD(S) 11 (15, 45X, 5A6) BLOCKING SURFACE IDENTIFICATION 
1-5 NSF(I) Surface number USE -1 TO END SURFACE INPUT. 

51-80 COMMET(I,J)Surface identification. 

CARD(S) 12 (9X, II, 20X, 5F10.3) BLOCKING SURFACE SHAPE PARAMETERS 


10 

ISF(I) 

Surface 

type 

31-40 

ALF(I) 

a 


41-50 

BMIN(I) 

^min 



51-60 

BMAX(I) 

^max 


Geometric surface limits. Refer to Fig. 4.3 

61-70 

GMIN(I) 

Ymin 


71-80 

GMAX(I) 

Ymax , 
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CARD(S) 13 (8X, 12, lOX, 6F10.3) SURFACE COORDINATE SYSTEM ORIENTATION 
9,10 ICS(I) Intermediate coordinate system. (0 if none is used.) 

21-30 RX(I) \ 

31-40 RY{I) f Surface system origin in reference or inter, coord, sys. 

41-50 RZ(I) ] 

51-60 CHS(I) X y toward x 

61-70 PSS(I) X toward z 

71-80 OMS(I) w y toward z 

REPEAT CARDS 11 THRU 13 IN SETS OF 3 CARDS FOR EACH SURFACE. 

END SURFACE INPUT BY A CARD 11 WITH NSF{I) = -1. 

CARD(S) 14 (8X, 12, lOX, 6F10.3) INTERMEDIATE COORDINATE SYSTEMS. 

9,10 NCS(J) Coordinate system number. USE -1 TO END ICS INPUT. 

21-30 RXI(J) I 

31-40 RYI(J) } Coordinates of intermediate coordinate system origin. 

41-50 RZI(J) 

51-60 CHC(J) X y toward x ) 

61-70 PSC(J) X toward z [ Rotations relative to the reference system- 

71-80 OMC(J) w y toward z J 

REPEAT CARDS 14 UNTIL ICS REQUIREMENTS ARE SATISFIED. 

END ICS INPUT BY A CARD 14 WITH NCS(J) = -1. 

CARD 15 (12E6.3) BAND MODEL PARAMETER SCALE FACTORS - USE ONLY IF ISCALA = 1 

1-6 ABSCL(l) Abs. coef. scale factor 

7- 12 DINSCL(l) Line density scale factor 

13-18 ABSCL(2) Repeat up to 6 pairs depending on the number of radiating species 

19-24 DINSCL{2) which are being used. The order must follow the IDENTC{I) 

list with unsealed species left blank. 

CARD 16 (1615) PLOT CONTROL CARD. USE ONLY IF NPLOT > 0. 

3-5 JPLOT(l) A series of NPLOT three digit numbers coded to define the 

8- 10 JPL0T(2) plots as follows: 

: First digit - abscissa control 

I = 1 wavenumber (cm“‘) 

I = 2 wavelength (ym) 

JPLOT(NPLOT) 

Second digit - ordinate control 

= 1 irradiance NSPEK 

1 or 2 watts/cm^-y-sr 
3 or 4 watts/cm^-y 
11 or 12 watts/cm-sr 
13 or 14 watts/cm 

= 2 transmissivity 
« 3 linear optical depth, X* 

= 4 collision fine structure 
parameter, ac 

= 5 Doppler fine structure 
parameter, ap 


For the final 
line-of-sight 
only. 


Rotations relative to reference or 
inter, coord, sys. 
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Third digit - constituent according to IDENTC(I) data 
= 0 used for second digit = 1 or 2 
= 1 first 

= 2 second Not applicable when second 
= 3 third digit = 1 or 2. 
etc. 


NOTE: Loop to read six sets of band model parameters in the order specified 

by the identc list: H20, C02, CO, HCL, HF, undefined. Blank cards are 
used for omitted species. 


CARD 17 (1615) 
lABSCO 


IFINST 


NACT(I) 

NUL(I) 

NUU(l) 

IDNU(I) 


BAND MODEL DATA CONTROL PARAMETERS AND LIMITS 

= 0 No data input for this constituent (blank card). Skip 
cards 18 and 19. 

> 0 Data input for F. 

= 0 No line density (1/d) input. This assumes program will 
compute a value as in the case of CO. 

> 0 Data input for 1/d using the same table format as for IT. 
Number of temperatures in the tables 

Lower wavenumber in the table (cm“^) 

Upper wavenumber in the table (cm“M 
Wavenumber interval for the table (cm"M 


CARD 18 (7E10.5) TEMPERATURE VALUES IN THE TABLE 
TEMAC(I,1) Temperatures (K) in ascending order. 

TEMAC(I,2) 


TEMAC(I,NT) NT = NACT(I) If NACT(I) > 7, two cards will be required. 

CARD 19 (7E10.5) BAND MODEL PARAMETERS FOR TEMAC(I,1) 

C0EF(J,K) K = 1, Absorption coefficient at NUL(I), TEMAC(I,1) 

NOTE; NT = NACT(I) 


C0EF(NT,K) K == 1, Absorption coefficient at NUL(I), TEMAC(I,NT) 

REPEAR CARD 19 FOR K = 2 TO (NUU(I)-NUL(I)/IDNU(I)+1) TO READ THE ABSORPTION 
COEFFICIENTS. THEN IF IFINST > 0 REPEAT THE RANGE OF K FOR THE LINE DENSITY 
1/d. 

RETURN TO READ CARD 17 UNTIL IT HAS BEEN READ 6 TIMES. 

PROGRAM EXECUTION COMMAND - Program begins execution on the first case and 

then reads data for subsequent causes starting 
at CARD 1. Execution terminates normally when 
CARD 1 is blank. 
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CARD FLOWFIELD INPUT BELOW TO BE INSERTED 
BETWEEN CARDS 8 AND 9 IF KASE < 0 ON CARD 2 


CARD IF (12A6) FLOWFIELD HEATIND 
1-72 HDG Heading 

CARD 2F (4A6, 12, 5X, 4E10.6) PROPELLANT OR GAS ID AND GEOMETRY (Refer to F1g. 4.5) 
1-24 BETA Propellant or gas identifier for information only 

25 ISYMSC = 0 Axisymmetric plume or a cluster with a center engine 

(assuming symmetry between center and outboard engines about 
the normal bisector of a line connecting the exit centerlines) 

= 1 A 3-D plume with no symmetry assumptions about a center engine 
31-40 HANG Included angle of a 3-D plume symmetrical sector. 

41-50 CTCD Distance from the plume property input axis to the Z-axis of a 

3-D plume system. 

51-60 AP = 0 No limiting cone used around the plume. 

> 0 X intercept of a limiting cone about the plume Z-axis. Used 
to speed interpolation by eliminating property searches for 
points outside the plume. 

61-70 BP Slope of the cone bounding the plume (dR/dZ). Must be > 0. 

CARD 3F (12A6) FLOWFIELD CONSTITUENT IDENTIFICATION 
1-6 CON(l) First constituent - first word. 

7-12 C0N(2) First constituent - second word. 

13-18 C0N(3) Second constituent - first word, 

19-24 CON (4) Second constituent - second word. 

Repeat up to a maximum of 6 constituents. 

REPEAT CARDS 4F THRU 6F AS REQUIRED WITH ZZ, ETA, AND R IN NUMERICAL SEQUENCE 

Card 4F (E10.6, 110) AXIAL POSITION - Units set by lUNITS on CARD 2. 

1-10 ZZ(I) Z coordinate of the data to follow 

20 NETA = 0 Terminates flowfield 

> 0 Number of n values for this Z. Use 1 for axisymmetric 

CARD 5F (E10.6, no, 15) ANGULAR POSITION 

1-10 ETA(I,J)n position (use 0 for axisymmetric flowfield) - degrees 
20 NETA = 0 Terminates flowfield 

> 0 Number of values for this Z. Use 1 for axisyironetric 
24,25 LS(IsJ) Index of the radial point adjacent to and at a smaller radius 

than a shock discontinuity. 

CARD 6F (9E8.5) RADIAL PROPERTY DATA - Units set by lUNITS on CARD 2. 

1-8 R(I,J,L) Radius - the first value (L=l) will be set at 0 

regardless of input. 

9-16 T(I,J,L) = 0. Return to Card 5F for next ETA or Card 4F for next ZZ. 

> 0. Temperature 
17-24 P(I,J,L) Pressure 

25-32 F(I,J,L,1) First constituent mole fraction. 

33-40 F(I,J,L,2) Second constituent mole fraction. Use the same order 

Continue for constituents on Card 3F. as on CARO 3F. 


All 




F=3 ^ is/1 nr e o h-i irv4c::. 


RTR 014-9 


Table A1.4 

FLOWFIELD GAS PROPERTY FILE 
GROUP I GENERAL INFORMATION (3 RECORDS) 


Record 

1 (HDG(I),I=1,12) 

ICASE 


Heading or title 
Case number 


2 (BETA(I),I=1,4) 
ISYMSC 

HANG 

CTCD 

AP 

BP 

3 NOCON 


Gas identification 


3-D geometry parameters (Ref. Card 2F, Table Al,3) 


Plume boundary cone parameters (Ref. Card 2F, Table A1.3) 
Number of constituents defined 


(C0N(2*I-1), Constituent name (2 words) and molecular weight 

C0N{2*I), (not used in program) for the NOCON constituents, 

CONWT(I), 

1=1, NOCON) 


GROUP II GAS PROPERTY DATA (NZ(1+2*NETA(I)) RECORDS) 


Record 

1 ZZ(I) Value of Z^ for data in records 2 and 3 in units 

specified by lUNITS (Card 2, Table A1.3) 

NETA(I) Number of ETA values used at this Z. This is the 

upper limit of a loop on J for reading Records 
2 and 3. For axi symmetric plumes NETA = 1, and 
the dimension limit is set by the parameter MN3. 


A12 



f=3 ^ rs/l "T ^ O HH I rvl c. 


RTR 014-9 


Table A1.4 
(concluded) 

FLOWFIELD GAS PROPERTY FILE 


Record 


2 ETA(I,J) 


Value of ri (degrees) for data in record 3. 


NR(I,d) Number of radial points in record 3. This is the 

upper limit (NOPS) for a loop on L for reading the 
radial points in Record 3. The dimension limit is 
set by parameter MR2 or MR3 . 

LS(I,J) Index (L) of the radial point in Record 3 which is just 

smaller than the shock radius RS. If LS(I,J) = LL, 
then R(I,0,LL) < RS < R(I,J,LL+1) 


3 (R(I,0,L) Radius Note: Units as specified by lUNITS 

T(I,J,L) Temperature (Card 2, Table A1.3) 

P(I,J,L) Pressure 

(F(L,N), N=1,N0C0N), Mole fractions in the same order as specified in 
L=l, NOPS) Group I Record 3. 


Records 2 and 3 are repeated NETA( I) times for each Record 1. The 
number of Records 1 (NZ) allowed for a single plume is unlimited, 
but for multiple axisymmetric plumes NZ is limited to MZ2. 
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A2 OUTPUT FORMATS 

The output formats to be described include the normal printed output, 
diagnostic messages, and binary files used for data storage and transfer. 

A2.1 PRINTED OUTPUT 

Printed output is provided In up to seven sections in separate subroutines. 
Output of these subroutines will be summarized in the paragraphs below in 
sequence of output, and examples of the usual output forms are illustrated in 
the sample problems (Section A3). 

FFPREP 

This subroutine is used only for reading and flowfields or scaling. When 
it is used, a portion of the flowfield data are output to assure proper operation 
and indicate the point that operation ceased if a difficulty is encountered. 

The data output is described in Table A1.4 as Group I and Records 1 and 2 of 
Group II. 

FLOUT 

The output in this subroutine is normally the first printed output. It 
contains general problem identification, geometry, and input limits. The data 
are obtained from the initial group of records on Unit 10 (Table A2.1) which 
are prepared by Subroutine FLOWIN. The general output features in sequence are: 

• Problem identification and data. 

• Flowfield identification, plume orientation, relative to the 
central coordinate system, and plume conical limits (intercept 
and slope). 


PRECEDING PAGE BLANK NOT r»iMF,D 
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• Statement of the units used for Input and output. The internal 
operation of the program is metric, and the radiation output is 
always metric as indicated by output formats. 

« Tabulated limits and increments for integration. 

• Location and direction cosines for the point-of-interest system 
(U,V,W) relative to the central coordinate system. 

• List of constituents in the input order (Table A1.3, Card 7) 

• Geometry input data for the point-of-interest and plumes. 

• Input data for blocking (shading) surfaces. 

• Input data for intermediate coordinate systems. 

• Identification of flowfield scaling performed on the input. 



Output of the flowfield properties is made in one of these subroutines 
depending on the flowfield geometry used. Normal output omits all interior 
points on each set of radial data to limit output while assuring proper flov/- 
field operation, but the complete set of property data can be obtained by 
setting the control flag NFLOW = 1 (Table A1.3, Card 2). The output is labeled 
by headings of Z, R, T, P for axial position, radius, temperature, and pressure. 
The mole fractions of the constituents are listed form FI to F6 (as required) 
according to the input order (Table A1.3, Card 7) printed in FLOUT (above). 

If the temperature range encountered in preparing the line-of-sight data is 
outside the range of 300K to 3000K, a message is printed at the end of the 
subroutine as a warning that the user should evaluate the temperature range of 
the band model data. If the gas temperature is above (or below) the table 
limits for the band model data, the upper (or lower) table value is used. A 
calculation time is printed for the flowfield interpolation, but because of a 
change in the MSFC UNIVAC 1108 system the number listed is currently SUP-seconds 
rather than actual time or CPU time. 
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ACDATA 

Band model parameters are output in this subroutine. First, the 
broadening coefficients for unidentified species (61 through 66) are listed 
for the first six species in the IDENTC data list order (IDENTC(I), 1*1,12). 
These correspond to the coefficients (Eq. 8) for the fraction of the 
plume gases which are not identified by the problem input. Following the 
broadening parameters, band model parameter scale factors are output if they 
are used. Next, the band model data being used are output after scaling or 
averaging over larger wavenumber intervals. Only the portion of the spectra 
to be integrated (ENUI to ENUF) Is listed regardless of the range of data 
which may be input. The wavenumbers listed are at the center of the v/a'venumber 
interval represented by the parameters, and this is the interval specified for 
spectral integration (Table A1.3, Card 5, DENU). Normally only the first and 
last wavenumber entry in each parameter set are listed to verify operation while 
minimizing output, but a full listing can be obtained by setting the control 
flag NABS=1 (Table A1.3, Card 2). 

SL6 

If Deck YN6H20 is used, the optical depth integration step size set by 
IRAD (Table A1.3, Card 2) is output. The step size is equivalent to 1/IRAD. 

The bulk of the normal output is a one-line summary of results for each line- 
of-sight. The values tabulated are defined below; 

• THETA - Value of the hemispherical elevation angle 6. 

• PHI - Value of the hemispherical azimuth angle <j). 

0 SHAPE FACTOR > This is not a true shape factor. It is defined by 

0 (j) 

SHAPE FACTOR - i i sin0 cos0 A0 Acfi/ir (A2-1) 

e-j 4’j 
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where the summation Is performed only for 6 and <|) defining 
lines-of-sight which intersect the plume. 

• FLUX - Incident radiation for the line-of-sight (DELTA) and 
the cumulative value (SUM) 

FLUX-DELTA = E Z N*(a>,s)[T(w,s-As)-T(w,s)]sin0cos(J) A0 Acj) Aw(A2-2) 
(jO s 

0 (f) 

FLUX-SUM = Z E (FLUX-DELTA) (A2-3) 

01 4>i' 

• RAD - Radiance for the line-of-sight (LOS) and the average for all 
lines-of-sight through this one (AVG). 

RAD-LOS = Z N“(w,s)[t(w»S-As)-t(w,s)] (A2-4) 

s 

0 (J) 

RAD-AVG = FLUX-SUM/ (Aw Z Z sin0 cos0 A6 A<1)) (A2-5) 

01 (f)^ 

The denominator in RAD-AVG is summed over all lines-of-sight 
rather than just those intersecting the plume as in Eq. A2-1, 
so RAD-AVG represents an average radiance for the space integrated. 

• AVG TRANS - This is a straight spectral average of the trans- 
missivity over all wavelengths for the line-of-sight defined by; 

AVG TRANS = Z t(w,s„,)Aw/2 Aw (A2-6) 

w w 

where sm indicates the upper limit of integration along the line- 
of-sight. 

• COMB OPT PATH - This is a sunmation of the optical paths in 
(cm-atm)c;-rp for all radiating species along the line-of-sight. 

It is denned by: 

COMB OPT PATH = Z Z [p- (s)c^ (s)273/Tj (s)]As (A2-7) 

s i 

where c represents the mole fraction and i represents the radiating 
species only . In Deck YNGH20 the title is changed to H20 PATH since 
H 2 O in the only radiating species. 

• LBLOCK - This is the identifing number (NSF(I) on Card 11) for the 
shading surface blocking this line-of-sight. If no blocking occurs, 
zero is printed. 

• SZERO - This is the lower limit for integration along the line-of- 
sight, Sq. It is set either at the input value (SZERO on Card 6) or 
by Subroutine SLIMIT. 
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• SMIN - This is the upper limit for integration along the line- 
of-sight, Sm. It is set either at the input value (SMAX on Card 
6) or by Subroutine SLIMIT. 


After all lines-of-sight are processed, a summary of "shape-factor”, 
radiance, and incident radiation flux is printed. The plume "shape-factor" 
was defined above in Eq. A2-1. The total "shape-factor" can also be expressed 
by Eq. A2-1 if all values of (j> and 0 are included rather than only values 
representing plume intersections. The exact totaV shape factor is the exact 
value of the Integration which is 

0f (j)f 

EXACT TOTAL SHAPE FACTOR = / / sin0 cos0 d(J)d0/Tr (A2-8) 

6-j <j>'j 

= (sin^0f-sin^0-j ) (^f-<j)j ) 

Overall radiance is defined by Eq. A2-5 using numerical - total - shape-factor 
in the denominator, while the plume radiance is defined similarly using the 
numeri cal - pi ume- shape-f actor in the denominator. The calculation time out- 
put for the SLG subroutine is in SUP-seconds for the MSFC UN I VAC 1108. 

SPCOUT 

Spectral data can be output per unit wavenumber or wavelength in either 
long or short formats as specified by NSPEK (Card 2). In the short fonn, only 
the first 5 values described below are output, but even in this case, the 
fourth and fifth values are for the final line-of-sight processed. In the long 
form, additional band model values are listed for the final line-of-sight. These 
values are generally valid only for the Curtis-Godson approximation using the 
exponential line strength distribution. They are not generally valid if Deck 
SLGMLT is converted to the exponential ly-tailed-inverse line strength distribution 
or if Deck YNGH20 is used for the Intuitative Derivative method. The values out- 
put are: 
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• WAVENUMBER - Wavenumber at the center of the Internal. 

• WAVELENGTH - Wavelength at the center of the interval. 

• IRRADIANCE - Radiation per unit spectral interval as specified by 

NSPEK (Card 2). The units in the format Indicate the option 
chosen. It should be noted that the conversion to radiation per 
unit wavelength is a function of the wavenumber interval used. For 
example, 

N;^{watts/cm^“Uim-sr)=(w^-(Aoi)/2)^)10“‘*N^(watts/cm^-cm"*-sr) (A2-9) 

• G - Transmissivity (x) over the final line-of-sight 

• OPTICAL DEPTH - The value of optical depth defined as 

X = Z Xi = -In r (A2-10) 

i 

is printed on line with the wavenumber. If the long form output is 
selected, additional lines are printed for each radiating constituent 
with Xi for the species identified or each line. 

• CONSTITUENT - Identifies the radiating species associated with each 
row of output. 

9 X* - The value of the optical depth for the linear limit of the 
curve of growth for the species identified (Eq. 15) 

X* = ke(a),i,Sm) u(i,Sm) {A2-11) 

• X/X* - Ratio of the optical depth for combined collision (Lorentz) 
and Doppler broadening (X) to the linear limit optical depth. 

• XC/X* - Ratio of the collision {Lorentz line) broadened optical depth 
{Eq. 10a) to X*. 

• XD/X* - Ratio of the Doppler broadened optical depth (Eq. 10b) to X*. 

• AC - Collision broadened fine structure parameter, aLg {Eq. 15b). 

• AD - Doppler broadened fine structure parameter, apg {Eq. 15c). 


VIEW 


This subroutine produces a printer plot of the line-of-sight radiance {the 
result RAD-LOS in the line-of-sight output). The plot is in a polar fonnat with 
the elevation angle (e) used as the radius and the azmiuth angle {({>) used as the 
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angle. The radiance is represented as tenths of the peak radiance by the 
digits 0 to 9. So if the line-of-sight radiance is N and the maximum radiance 
computer for a line-of-sight is Np, then an integer D defines a radiance interval 
DNp/10 < N £ (D+l)Np/10 (A2-12) 

The value of Np/10 is printed above the plot. 

Lines-of-sight which are blocked and do not enter the gas are represented 
by the character B, but if radiating gas is encountered before blocking occurs a 
digit is printed and no indication of the blocking is given. Blanks in the plot 
generally indicate lines-of-sight missing the gas, but they also may occur locally 
due to the lack of precision in converting the polar plot format to a rectangular 
grid with one symbol/line resolution. If lines-of-sight are closer than the 
resolution of one symbol, the results are averaged. But if they are widely spaced, 
no interpolation is performed, so blanks will occur. 

A2.2 DIAGNOSTIC MESSAGES 

Diagnostic messages are provided in 7 subroutines of GASRAD and the plot 
program (PLOTIT). The messages and corrective actions usually implied are 
listed below for each subroutine in the usual sequence of subroutine use. 

MAIN 

NUMBER OF CONSTITUENTS REQUESTED (KKK) EXCEEDS THE NUMBER ALLOWABLE (LLL). 

ACTION: Reduce input (NOKON, Card 2) or increase dimension parameter NK. 

(ENUF-ENUI)/DENU+1 IS GREATER THAN THE MAXIMUM NUMBER OF SPECTRAL INTERVALS 
DIMENSIONED. 

ACTION: Reduce input interval, increase DENU, or increase dimension parameter 

MNU. 
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******S0RRY, CASE LLLLLL IS NOT ON THIS TAPE. RUN TERMINATED. 

ACTION: Input case LLLLLL (KASE, Card 2) could not be matched to any flowfield 
stored on Unit 8. Check card and tape input. This message will only be written 
if the final record on Unit 8 is a 72 character string of *. Otherwise an end- 
of-file will be indicated with no diagnostic. 

******CASE NUMBER LLLLLL DOES NOT CONTAIN THE CONSTITUENT 'AAAAAABBBBBB' 

ACTION: Check input case numbers (KASE, Card 2) and spelling of input constituent 

names on both Card 7 (KON) and the flowfield (KONAME). Note that the names are 
all two words (2A6), and the flowfield tape (Table A1.4) has a list of 3 word 
items since the molecular weight was included (although it is not used in the 
program) . 

THE CONST ITUENT-AAAAAABBBBBB- IS NOT RECOGNIZED BY THE PROGRAM. CHECK AGAINST 
NAMES IN THE IDENTC LIST. 

ACTION: Check as directed with the data list specifing the constituent names 

(IDENTC(I)) in Subroutine MAIN. 

NO. OF CONSTITUENTS ON TAPE (IT) EXCEEDS STORAGE AVAILABLE (NKTP=K) 

ACTION: Increase dimension parameter NKTP. The storage required by this 

parameter is only (3+MR2[or MR3])*NKTP. 


FLOUT 


All messages in this subroutine concern inconsistencies detected when using 
data on Unit 10 to restart or run a repeat case. The action in all cases is to 
check the input indicated. 

LOS TAPE CASE NUMBER ICASE DOES NOT MATCH INPUT CASE NUMBER KASE. CASE TERMINATED. 
INPUT SMAX XXXX.XX EXCEEDS LOS TAPE VALUE OF XXXX.XX 
INPUT DTHETA XXXX.XX EXCEEDS LOS TAPE VALUE OF XXXX.XX 
INPUT DPHI XXXX.XX EXCEEDS LOS TAPE VALUE OF XXXX.XX 

THETA LIMITS DO NOT AGREE - INPUT VALUES ARE THETAI = XXXX.XX AND THETAF = XXXX.XX 
PHI LIMITS DO NOT AGREE - INPUT VALUES ARE PHI I = XXXX.XX AND PHIF = XXXX.XX 

FL0W3D 

NUMBER OF ETA- CUTS EXCEEDS MN3. NETA(I.O) = N Z(I) = XXXX.XXXX 
ACTION: Correct input flowfield or increase dimension parameter MN3. 
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NUMBER OF RADIAL POINTS EXCEEDS LIMITS. NRMAX = MN3 NOPS = N 

ACTION: Cosrect input flowfield or Increase dimension MR3. 

NO. OF LINE-OF-SIGHT POINTS EXCEEDS LV, NEED LARGER DELTA S. 

ACTION: Increase DS(Card 6)» reduce [SMAX-SZERO] (Card 6), or increase 

dimension parameter LV. 

NEGATIVE PROPERTIES ON THE LINE-OF-SIGHT AT THETA = X.XXXE±XX RAD AND 
PHI = X.XXXE ±XX RAD 

ACTION: This message is followed by a labeled listing of line-of-sight 
properties. Examine it to determine the cause of negative properties and 
correct. Note that 6 and <|) are in radians and all parameters in the line- 
of-sight property list are in the cm/K/atm. 

RADIUS POINTS NOT IN ORDER OR ARE EQUAL 

ACTION: This is a warning message. No action is taken to correct the 
situation. The input should be examined to determine if the error will 
have a serious effect on linear interpolation. 

***WARNING - GAS TEMP RANGE OF XXXX.X TO YYYY.Y IS OUTSIDE THE RANGE OF 
COEFFICIENT TABLES 

ACTION: This warning message is based on the usual table range of 300K to 

3000K (actual ranges are not defined at this point in the program). No 
corrective action is taken and subroutine SLG does not extrapolate outside 
;the band model parameter tables, so the nearest table value will be used for 
k and 1/d if the gas temperature is outside the table limits. 


FLOWAX 


Messages are the same as FL0W3D except: there is no test on ETA, the 
parameter MR2 replaces the parameter MRS, and a test is made on the number 
of Z-cuts for multiple axi symmetric plumes. 

******Z-CUTS EXCEEDING MZ2 CANNOT BE LOADED WITH MULTIPLE PLUMES, INCREASE 
PROGRAM DIMENSION MZ2 OR SET IZF = MZ2. 

ACTION: Remove Z-cuts from the flowfield to reduce it to the MZ2 limit or 

follow the actions recommended in the message. 


ACDATA 


**********UPPER LIMIT HAS BEEN ADJUSTED TO GIVE AN INTEGRAL NO. OF INTERVALS 
**********NEW ENUF = XXXXXX. 
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ACTION: Evaluate results based on altered spectral range. Since program is 

not designed to handle a fractional wavenumber interval, the limits of inte- 
gration (wf-w-j) must be an integral number of intervals (Aw). 

NPLOT MUST NOT BE GREATER THAN 30. 

ACTION: Correct input or increase the dimension of JPLOT. 

NO. OF TEMPS EXCEED MAX. NACT(I) = L NTMAX = MT = K 

ACTION: Reduce the number of temperature values in the input band model 
parameter table for constituent I (of the IDENTC list) from L to K; or 
increase the dimension parameter MT to L. 

ABS COEF DATA BEING LOADED FOR IDENTC LIST SPECIE I EXCEEDED ALLOWABLE ARRAY. 
INDXCO = J NCOMAX = K 

ACTION: The dimensioned number of vmvenumber/species combinations (MCO or 

NCOMAX) was exceeded by constituent I (of the IDENTC list) which would have 
required INDXCO positions to be loaded. Reduce the total wavenumber/species 
positions to less than the NCOMAX specified or increase the dimension parameter 
MCO. 

WAVENO. INTERVAL IDNU(I) = J EXCEEDS DESIRED IDENU = K 

ACTION: The input wavenumber interval J specified for constituent I (of the 

IDENTC list) is larger than the integration interval. Since the program is 
not designed for this condition, the integration interval must be increased 
to 0 or the input table of band model parameters needs to be corrected. See 
next note also. 

WAVENUMBER INTERVAL MUST BE AN INTEGRAL NO. OF THE CONSTITUENT INTERVAL. 

IDENU = J IDNU(I) = K 

ACTION: If the wavenumber integration interval (DENU) is larger than the table 

interval (IDNU(I)), the procedure for averaging spectral parameters is invoked, 
but it requires that DENU be an integer multiple of IDNU(I). Modify either the 
input for constituent I (of the IDENTC list) or the integration interval. 


SLG 


The diagnostics are the same for either deck (SLGMLT or YNGH20). 

******THE SUM OF SPECTRAL INTERVALS FOR ALL CONSTITUENTS IS KKKKK. THE 
MAXIMUM ALLOWED IS LLLLL 

ACTION: Reduce the spectral interval or number of species, or increase the 

dimension parameter MOD. Since all input has been read at this point, the 
program returns to MAIN to read the next case instead of stopping at the error 
message. 
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******NEGATIVE DM0 AT WAVENUMBER = JJJJJJ. DM0 = -X.XXX±EXX 

ACTION: This is a warning of a negative optical depth. It is treated as 
a transparent (At^O) interval and the program continues. The cause and 
importance of the negative value should be evaluated. It can be caused 
by a negative pressure or mole fraction, but this should have been detected 
causing an error message in FLOWAX or FL0W3D. 


PLOTIT 

Messages based on 3 digit value for OPLOT identified here as UK. 

ABSCISSA INDICATION MISSING 

ACTION: Correct input. Either JPLOT < 0 or I = 0 for the missing plot. 

ORDINATE INDICATION MISSING 

ACTION: Correct input. Missing plot had J = 0. 

FINE STRUCTURE PARAMETERS NOT USED FOR CARBON. PLOT LL DELETED. 

ACTION: Correct input. Designed to notify deletion of plots of ac and ap 

If they are requested for carbon (147 or 157), but this message will also 
be printed for any case with K ^ 7 and J ^ 4. 

A2.3 BINARY FILES 

Binary files are used on Units 3, 4, 8, 10 and 11 to specify flowfield 
properties (8 or 11), to store line-of-sight property data (3, 4, and 10), to 
transfer line-of-sight radiance from Subroutine SLG to Subroutine VIEW (4), 
and to transfer spectral data from Subroutine SPCOUT to the plot program 
(PLOTIT) (4). These functions are summarized in Table Al.l and the format 
of the flowfield (Unit 8 or 11) was described in Section A1.3 as input. Formats 
for the line-of-sight data and the two other data transfers will be described here. 

LINE-OF-SIGHT DATA 


Unit 10 stores line-of-sight data for transfer from the flowfield Subrou- 
tines (FLOWAX and FL0W3D) to the radiation subroutine (SLG), but it also stores 
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input data to serve as a restart resource for input with IFLOW = 7 (CARD 2). 
Unit 3 and 4 may also be used in processing lines-of-sight when the flowfield 
must be processed in sections, but these units do not have the input data 
(Records 1 through 4) written on Unit 10 for output (in FLOUT) and restart. 
Contents of the line-of-sight data file are described in Table A2.1. 

RADIANCE DATA 

At the end of the PHI loop in subroutine SLG, a record is written on 

Unit 4 for each line-of-sight. Each record contains the following values: 

PTHETA - Elevation angle 6 in degrees 

PHI - Azimuth angle in radians 

RADIOS - Radiance for the line-of-sight (Eq, A1.4) 

MBLOCK - Missed/Block flag 
=0 Missed the gas 
= 1 Blocked before reaching the gas 
= 2 Radiation computed 

The variable MLINE is incremented as each record is written, so that it can 
be used in Subroutine VIEW to read the correct number of records. 

SPECTRAL DATA 

At the end of Subroutine SPCOUT a single record is written on Unit 4 
to provide spectral data to the plot program (PLOTIT). The arrangement of 
the record is described in Table A2.2. 


A26 



(=? ^ ts/l T- ^ O I— ! i rvl o. 


RTR 014-9 


Table A2.1 

LINE-OF-SIGHT DATA FILE 


GROUP I GENERAL PROBLEM SPECIFICATION (4 RECORDS) 

BE -Q.QBa 

1 DATE DATE OF TAPE GENERATION 

KASE, IFLOW, NOKON, ISHOCK, UNITS, ISCAEF, IZI, IZF, THETAI, THETAF, 
DTHETA, PHI I, PHIF, DPHI, TANGLE, ENUI, ENUF, DENU, SZERO, SMAX, DS, 
TDIFF - Input data from cards E through 6 (Table A1.3) 

HD6(12), BETA(4), ISYMSC, HANG, CTCD, AP, BP - Input from flowfleld 
(Table A1.4) 

NENG Number of plumes 

NEPl NENG + 1 

(NUMCON(K), K=1,N0K0N) - Order in IDENTC list of the K-th constituent 
in the input list (card 7). 

IBLOCK = 0 If no blocking surfaces input 

= Index of initial blocking surface in the surface input 
array (should be NENG + 2) 

(TX(N), TY(N), TZ(N), ((TCS(J,K,N),J=1,3), N=l, NEPl) - Coordinates 
and direction cosines of the point-of-interest (N=l) and 
and the plumes in the central coordinate system. These 
are computed in BLOKIN. 

2 RSCAL, TSCAL, PSCAL, FSCAL(NKTP) - Scale factor input from card 8 

(Table A1.3) 

3 (KON(I),J = 1, N0K0N2) - Constituent identification input from card 7 

(Table A1.3). Two word names are used, so NOKON2=2*NOKON. 

4 NBLOCK Index of final surface. 

JBLOCK Index of final intermediate coordinate system. 

(NSF{I),(C0MMET(I,K),K=1,5)ISF(I),ALF(I),BMIN(I),BMAX(I),GMIN(I), 
GMAX(I),NCS(I),RX(I),RY(I),RZ(I),CHS(I),PSS(I),0MS(I),I=1,NBL0CK) - 
Surface input from cards 9 through 13 (Table A1.3). 

(ICS(j), RXI(J), RYI(J), RZI(J), CUC(J), PSC(J), OMC(J), J=l, JBLOCK) - 
Intermediate coordinate system input from card 14 (Table A1.3) 
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GROUP II LINE-OF-SIGHT DATA (See note below) 


RECORD 



ALL 

EXCEPT 

LAST 

LSET 

Line-of-sight sequence number is [LSET|. A negative 
value of LSET indicates the line-of-sight missed the 
gas. 


LPT 

Number of points on the LOS. 


LBLOCK 

Identification (NSF(I)) of blocking surface. If none, 
LBLOCK = 0. 


SINKO 

(sin0 COS0) 


THETA 

Elevation angle 6. 


PHI 

Azimuth angle <t». 


DSFACT 

(sine COS0 A0 A^j/ir) 


DX, XY, 

DZ - Dimension increments in the central coordinate system 
corresponding to DS. 


SMIN 

Upper limit of the LOS length. 


(XT(L), 
LPT) - 

YT(L), ZT(L), TWT(L), PWT(L), (FWT(L,K), K=l, NOKON), L=l, 
Location in the CCS and properties (Temperature, pressure, and 
mole fractions) for points on the LOS. 


DPMI 

A4) for this LOS. This may be different from the input DPHI 
because of the operation of the input value of TANGLE (Card 
4 of Table A1.3) 


SZERO 

Lower limit of the LOS length. This may be different from 
the input value of SZERO because of fitting the LOS length 
to the plumes boundaries in Subroutine SLIMIT. 

LAST 

LSET 

= 0 


All other data identical to the previous LOS. 


NOTE: If TANGLE = 0, the number of lines-of-sight is 


[ (THETAF-THETAI )/DTHETA3[PHIF-PHI I )/DPHI] . 
The number will be less for TANGLE > 0. 
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Table A2.2 

SPECTRAL DATA FILE 
(1 RECORD) 

L Number of wavenumber intervals = (ENUT-ENUI)/DENU + 1 

I Number of optical depth values. 

NKC Number of radiating species for which the program is set up. 

It is equivalent to the dimension parameter MKA. 

NODC Number of spectral intervals for carbon. It is equivalent 

to the dimension parameter MODC. 

NPLOT Number of plots from Card 2 (Table A1.3) 

(ENUV(LL), Spectral parameters •* Wavenumber (cnrM 

ELAIWCLL), - Wavelength (pm) 

SFLXTO(LL), ■■ Spectral flux in units set in ACDATA 

for output. Determined by NSPEK (Card 2) 

GNEW(LL,LL=1,L) ■■ Transmissivity of final line-of-sight 

(FO(II), Spectral parameters for each species - optical depth in weak line 

limit. 

ACO(II), - Lorentz line fine sturcture 

parameter, 

AD0(II),II=1,I), - Doppler line fine structure 

parameter, aQg 

(JPL0T(N),N=1, NPLOT) - Plot specifications from Card 16. 

(TITLE(N1),N1=1,12) - Title from Card 1. 

(NUMC0N(N2),N2=1,10) - Position on the IDENTC list of constituent N2 
in the input list (Card 7). 

(NUL(N3), Wavenumber interval of data for constituent N3 - lower limit 
NUU(N3),N3=1,NKC) - upper limit 

(F0C(N4),N4=1,N0DC) - Optical depth for carbon. 

NOKON Number of constituents for this case from Card 2. 
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A3 SAMPLE PROBLEMS 


2’. 


u, 

6 . 


1 INCREASE DEFAULT SIZE OF TEMP. PROG. FILE 


ASSIGN EXTERNAL STORAGE DEVICES 


8 . 

— 

in. 


1?! 


16. 


iREriirGlsR^^ ASSIGN AND REWIND INPUT TAPE 



3ASG,T TPF».,r/l/POS/5 

dASG,T 4,F///?500 

e-f- T - 

dASG,T 10. ,F///2500 



aFREE GASRAD 

- a xer e ASRA t> 

SAMPLE PROBLEM 1. HEAT TRANSFER -- SSME PLUME RADIATION TO ET BASE CTR. 

^ j 0 — ^ ^ 


COPY IN AND FREE INPUT mPE 









IP. 

t^T“ 

20,. 


8. 

1000 . 
-Q- 


7*1. 

"«T OV 

9800. 
1 8 00 ^ ■ 


<1 . 



400. 








H2ol (G) N2(G) 

a ADtr^STtt 


22! 


0 

6 


FLOWFIELD CARD IMAGE DATA ELEMENT FROM INPUT TAPE . - 

1<I3<I.S 0. 63.5 leo. 90. Q. POINT-OF-INTEREST 

- 0. O'? 1B7 • 0'« ' D. 0. \ V __ 

0. 0. 157. 0* 0* 0. i PLUME ORIGINS 

Ott-J 


"•©T- 


157 


-Or- 


-ax- 


ENG 2 BELL 

0^- 


-seoir 


H 2 O BAND MODEL PAI^ETER CARD IMAGE DATA ELEMENT FROM INPUT TAPE 



BLOCKING SURFACES 


jlNTERMEDIATE COORDINATE SYSTEMS 


5 BLANK CARDS TO INDICATE NO DATA FOR OTHER GASES 


HF-AT TRANSFER FOR A 0.04 SCALE MODEL \ 

1 i REPEAT CASE BY SCALING DATA ON UNIT 10 


j5 BLANK CARDS FOR BAND MODEL INPUT 
ONE BLANK CARD TO TERMINATE RUN 


INPUT LISTING FOR SAMPLE PROBLEMS 1 AND lA 
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SSME FLIGHT 

N02ZLF SEA 

LEVEL FLUME ’ 

BETA LH2/L0?/AIR 

ISYI 

N0C0N= 4 

H201IG) 

H2C6) 1 

• UO 

1 



• on 

20 

5.0D 

1 



.no 

- 23 

1 0 . 00 

1 



.00 

21 

IS. CO 

1 



• oc 


2C.0D 

1 



♦ no 

20: 

3L.00 

1 



. GO 

19 

3s.no 

1 



.GO 

20 

ML. 00 

X 



UD 

17 

45.00 

1 



.GO 

17 

5f .on 

1 



.CO 

18 


52.00 
52. XC! 
6 L . CIO 
7C.U0 
8C.P0 
9r.00 


inc.oo 

1 


i CP 

l"'5.Gn 

1 


.00 

115. rn 

1 


.uc 

i 2 c.cn 

i 


.on 

UC.DO 

1 


.lo 

i4r,0G 

X 


.no 

15 ,.U0 

1 


;CP 

16 l. in 

1 


.no 

ISL.GO 

1 


.L-n 

ioj.cn 

X 


9 jO 


02IGI 

i 


C78 - IC»SE: 1 

.00 CTco= 6.25o«ol BP: l.250-t;X 

N2.GI-- — CONSTITUENT ORDER MUST MATCH MOLE 
FRACTION ORDER ON INPUT 


ZZ(I), NETA(I) 

ETA {ETA(I,J),NR(I,J),LS(I.J),J=1,NETA(I) 
RADIAL POINTS AND PROPERTIES AT EACH ZZ(I) 
ARE OMIHED 


PARTIAL ECHO PRINT OF FLOWFIELD 
INPUT IN SUBROUTINE FFPREP. 

ONLY OCCURS WITH CARD INPUT OR 
FLOWFIELD SCALING 


13 


1 
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CO 

CO 



1 

1 

1 



1 

1 



I 

V 

1 

1 

1 

1 

I 

t 

1 _ 

1 

1 



I 

1 

1 

I 

X 

X 

X 

1 

^ SAMPLE PROBLEM 1 
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.0 


7nn*U0 

1 


• 10 

72r.Qn 

1 


.00 

760*00 

1 


*L0 

8nu*oo 

1 


• UO 

8 4 j:.ud 

1 


*cr 

86w.uu 

1 


.00 

9rr;.n0 

1 


i tiO 

92C.0P 

1 


iCO 

94r.oo 

1 


iilV 

lono.ao 

1 


i 00 

uun*OD 

J 


• 00 

l2Q0.i3O 

1 


.on 

*00 

n 


--IS- ^ — 

16 - 1 

16 1 

IS 1 — 

15 i 

15 1 

-16 r-- 

16 1 

16 1 

- 16 ^ 1 

16 1 

16 I 



SAMPLE PROBLEM 1 
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3 AXI<;YHHEtRIC plumes using cose 


t WITH exit CTRLINE LOCATIONS *ND DlR^rTIONS tMntrATED 


PUN date 03Q68Q 



padiaTion problem identification - sample I 

PROBLEM 1, HEAT TRANSFER - SSME PLUME RADIATION TO 

ET BASE CTR. 


•c 

95 


FLOW FIELD IDENTIFICATION SSHg; FLiGHT ^OZZLE 

GAS IDENTIFICATION 

SEA level plume - 

- LhZ/L0?/AIR 

TKE 5/4/78 




t% 


ENG^ 

no/ 

LOCATION 

X 

1 OF EXIT CENTER 
Y Z 

OiRECTlO,^ 

X 

CqSINES 

Y 

OF AXIS 
Z 


C9 


1 

_ 2 

1595.92 

1622.78 

.□0 R85.28 

-53.00 369.90 

.96126 

.98481 

.00000 

•00000 

.27564 

.17365 


tK 


' 3 1622.78 

FLOW FIELD CONICAL 

53.00 369.90 

Limits - Intercept 

.98461 
= 6. 250*01 

.00000 

SLOPEr 

,17365 

1.250-01 


•» 

•* 


~ 1 rucTU t u Tururc _ rrMorDATtior tm nirro 

rrc DAUi/Tur — ddtc 

Clior TM DCfi-A 

liMTTC CnO TMD IT 







'SUrfL AN 

— VM X 1 

o 1 xm V 1 




LIMITS OF INTerRATtON VARIABLEc INITIAL VALUE 

FINAL VALUE 

INCREMENT 




C» 


z-index 

THETA 

1 

8.00 

55 

7N.0D 

4.0000 




0» 

m 


PHI 

NU 

• OQ 
iCOD« 

40.00 

A M M A 

7 O VI U • 

R‘0000 

4uO. 




it 


s 

TEMP step Size 

• 00 

1800.00 

6.0000 

200.0000 




« 

:s> 

Ca> 

Yangle (6eg) 



4.0000 




ft 


X y 

point location mSN.SO .00 

z 

1 nrATTAM TW rrMTDfli rnnonTMATc cvctcm 


a 



&.wwril AVM Xil 


./ixwxnni t. 

^ i L.ri 


tr 


COSINES OF U -•LOOOO -iODODD 

i.OOQOO 






t>t 

» 


9 nrroiss -i.oaoes 

W l.OOOOn .00000 

TffMISQ 

.00000 

y DIRECTION COSINES RELATIVE lU CENIKAL EUUKU. 

It 

• 


CONSTITUENTS H20KG) N2(GJ—. 

— NO PARTTCUI AR ORDFR RFOUTRED 




« 

tx 









tL 

M 


COORDINATE DATA INPUT FOR THE POINT OF INTEREST AND PLUMpS 

ICS NO. X Y Z ' CHI PSI 

OKr^^A 



- - - 

OC 

«i 


pojNr OF Interest 

0 IRJR.Sn .00 

65.50 

180.00 90.00 

• 00 




•t 


PTOHE NO. I 

5 .00 .00 

X57.00 

.00 .no 

• 00 


. 


M 


PLOft'C f(6. 2 

6 .00 .00 

157.00 

.00 .no 

• 00 






pXumI NO. 5 

7 .00 .00 

157.no 

.00 .no 

• 00 


■■ 



O 

_ _ . - - __ . . ___ __ _ _ 

t| 

ho 

t 

* 


SAMPLE PROBLEM 1 



BLOCKING SURFACES 
NO. TITLE 

TYPr 

_ ALF 

jimiiii^^s^^viiiini 

BMAy . 


6MAX 



TCS 

■mfii 

RY 

RZ 

CHS 

PSS 

OMS 

S2r)l ENG 2 BELL 

7 

4^23^^♦00 

2.7l7<**01 

1.2884*02 

• Q 

360.0 



6 

■ LILIl — 

O.nOQD 


*0 


- .0 _ . 

1200 ORBITER LOHfR FUSELAGE 

3 

D.OOQD _ 

D.0000 

1.2407*03 

-6,D 

6.n 



n 

2.7700*02 

0*0000 

2.7670^02 

9o;o 

.0 

• 0 

36'11 BODY FLAP 

3 

0.0000 

4.9373402 

5*9333+02 

-12.6 

X2.6 



0 

2.1038*03 

D.DQOO 

3.3922+n2 

270.0 

.-p 

-.--6,. 3 


intermediate coordinate system for blocking surfaces 

ICc Rvl RYI RZI CHC PSC 

OHC 

5 1.4450*03 

0.0000 

4.4300+02 

180.0 

74,0 

• 0 

6 1.4682*03 

-5.3nOD*01 

3.4264+02 

180.0 

eolo 

_ _,_*0 _ , , ,, .. . _____ , 

7 I. 4682+03 

5.3000*01 

3.4264402 

len.o 

80.0 

• 0 


on F IELD .S C »L£ .£ A<JORa..-M £ff C, .USja.. 


I 


K 

CC 


oc 
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> 

CO 

■VJ 


FLOW FIELD 


z 

R 

~T 

r-" — 

-FI 

-- F2 

.COQ 

• ndO 

l4C6* 4 

39i6400 

. 7702 

iunoo 


45.39D 

558.0 

2116.2200 

.COlO 

.7900 

5.00G 

.COO 

1366. C 

34.8700 

.7702 

.aoDO 


64.570 

558.1 

?116.2?00 

.0000 

.7930 

in.DOO 

.POO 

1326TV9 

307 6900” 

17702 

-.aooQ 


62.580 

558.7 

2116.2200 

.0002 

.7898 

15.QOO 

.noo 

1866. 1 

ne.scno 

.7702 

.anaa 


61.160 

562. 4 

2116.2200 

.0014 

. 7885 

2Q.D00 

.nao 

1866.2 

129;C000 

.7702 

.0000 


60.180 

574. 3 

2116.2200 

• COSO 

.7849 

3D. JOO 

.POO 

"isiar? 

risjcToo 

17702 

iornn 


73 .760 

558.1 

2116.2200 

.unoo 

.7900 

35.cno 

.POO 

1343.9 

128.8900 

.7702 

.0000 


72.200 

558.3 

2116.2200 

.0001 

. 78 99 

«*O.OOD 

• POO 

1776.** 

110.0700 

.7702 

. onno 


72.630 

558.9 

^116.2^U0 

I.r003 

.7897 

**5.000 

• 000 - 

" lBf)l .“ l* 

' - ir7764fT0 

- 17702 

iOCOO 


71.870 

559.9 

2116.2200 

• 0006 

. 7894 

SO.OOQ 

.000 

1767.5 

108.7700 

. 7702 

.0000 


70.610 

561.7 

2X16.2200 

.0011 

• 7890 

52.000 

.noo 

1765.7 

108.3100 

• 7702 

.GOOD 


70.610 

561.7 

2116.2200 

• 0011 

. 7890 

52.irO 

.noo 

" 56n0. 

”9-35<JlOrOCt-- 

* 17700 

icooa 


70.610 

561.7 

2116.2200 

• UOll 

.7890 

60.000 

.POO 

57oq.X 

5400.0000 

17700 

.0000 


67.920 

566. 1 

2116.2200 

.C02 3 

.7880 

70.000 

.noo 


4712.0000 

17700 

iOOOO 


64 .630 

574.2 

2116.2200 

.C04 3 

.7865 

eo.noo 

.noo 

5 6rn^a 

4602 i 0000 

-.7700 

iOPQO 


62.380 

587.1 

2116. 22CO 

.0072 

. 7844 

90. GOD 

.POO 

5 6pf3i 0 

4686.0000 

.7700 

.GOOn 


79.440 

558.3 

2116.2200 

• cool 

. 78 99 

100.000 

.POO 

555c. C 

4o33icnno 

.7700 

.0000 


80.700 

556.8 

2116.2200 

• CC03 

.7893 

105.000 

.POP 

5 530Vtr 

-"35r3;cnoa 

17700 

- iOOOD 


80.930 

559.6 

2116.2200 

.0005 

. 7896 

115. ono 

.POO 

53nO*D 

2386. CCOO 

. 7700 

.0000 


81 .540 

561.1 

2116.2200 

• OHIO 

. 7893 

120.000 

.000 

5170. C 

1779. coon 

. 7700 

.0000 


82.P6U 

563.6 

2116.2200 

• 0016 

• 7888 

13C.C00 

.PQO 

4 970.0“ 

* X U79 5C000 ” 

17700 

1 coun 


82.300 

567.3 

2116.2200 

• 0026 

.7881 

i«*a.Qoo 

.noo 

507n.0 

1402. onnn 

.7700 

.0000 


82.510 

572.7 

2116.2200 

.0039 

. 7871 

150. ODD 

.PflO 

5210.0 

1958,Cin00 

. 770 1 

.0000 


96.750 

558.1 

2116.2200 

• tool 

.7900 

160.000 

• non 

SZlpiO 

- 19181 C'^CO 

17700 

iOCOO 


83.350 

590.4 

2116. 22C0 

.U080 

. 7843 

180. onu 

.000 

5210.0 

2U35.0Cn0 

. 7700 

iDOOO 


96 .660 

558.9 

2116.2200 

.noo 3 

.7898 

190.000 

.'^nn 

519jjiC 

17981 coon 

. 770U 

• an 00 


SAMPLE PROBLEM 


SHORT FORM OF FLOWFIELD 
OUTPUT IN SUBROUTINE 
FLOWAX. INTERIOR POINTS 
ARE NOT PRINTED. 
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OJ 

CO 


p O 
O fci 


‘ w 

p 
L*^ H:^ 
? H C?:-i 

c2 


PLOW FIELD ICOnTINUED) 

2 R T P -ft ¥2 F3 f4 F5 F6 


190. uno 

97.570 

S59i6 

2116i2200 

iCODS 

*7896 

ano.ono 

.TOO 

5210. C 

1932. QOOO 

.7700 

.0000 


97.780 

560i7 

2116i22m 

*0009 

, 7894 

210. UOO 

• nuQ 

5480.0 

3123.0000 

.7700 

• uoon 


98i470 

562rr 

— tnev22C0 * 

’*0013— 

“ * 78 91 

220.000 

.000 

5 55q. C 

3760. UDOO 

.7700 

.0000 


98.860 

567i5 

21l6i2200 

i C02 7 

.7881 

230.000 

.000 

492Q.0 

2182 . anoo 

. 7700 

.0000 


99.040 

571^ 7 

2116i220n 

*0038 

*7873 

250.000 

.000 

4915.0 

2184.0000 

. 7700 

.0000 


99.060 

— 2ri6T^2aa~" 

—Toer- 

-"*^78 5 3 

280.000 

.000 

492o* D 

2286. OCOO 

.7700 

• 0000 


123.220 

558.-8 

21l6i 2200 

*0003 

*7898 

290.000 

.000 

465p#0 

1332.0000 

.7700 

.0000 


I23i660 

559i2 

2116^2200 

*0004 

i7897 

300.000 

. 000 

45Sq. C 

1104.0000 

.7700 

• 0000 


123;870 



”21t«rv2?P0— 

— ?000<r 

*7896 

iio#ano 

.000 

4808.0 

1802.0000 

. 7700 

.JO no 


123i940 

56pi6 

2116^2200 

iC'008 

*7894 

320.000 

.000 

99ol.C 

2159. QPOO 

.7700 

«0000 


129.130 

561^.6 

2116i2200 

icon 

*7892 

350.000 

.000 

4929.0 

2218. UOOO 

. 7700 

.00 00 


I23i280 

^6tFTT^ 

'-MlfrT2200~ 

*0024- 

*7883 

360.000 

.000 

486Q.C 

1911.0000 

. 7700 

• ucoo 


122.920 

569s5 

21l6i2200 

i0030 

*7879 

370.000 

.000 

5q49.0 

2798. COOO 

. 7700 

.onoo 


122 .780 

572 i 7 

2116i220n 

*0038 

.7873 

380.000 

.000 

4755.0 

3438.0000 

.7700 

• CJOGO 


122 i87Q 

576 r7 

-2116*2200- 

“i0047 

*7867 

390.000 

.000 

4755.0 

3483.0000 

. 770P 

.0000 


I22i950 

58U3 

21l6i22Q0 

.6057 

*7860 

400. UOO 

.000 

4532.0 

2063.0000 

.7700 

.0000 


123.260 

586i9 

2116*2200 

*0068 

*7852 

420.000 

♦ 000 

4463.0 

1807.0000 

• 7700 

.0000 


124.540 

59359 

- 2116.2200 

suoei 

*7843 

440.000 

.000 

4 540. 0 

2119. troD 

.770U 

• QOOO 


170.720 

558 i5 

2116*2200 

iU002 

*7899 

460.000 

♦ con 

4282.0 

1210 . cnoo 

. 7700 

.0000 


170.920 

558i9 

2116*2200 

*0003 

.7898 

480.000 

.nou 

4 54 0. t: 

2L26.G000 

. 7700 

• onoa 


171.140 

559i6 

*2116*2200 

*0005 

*7897 

500. cno 

.f^no 

4644.0 

2637. GOOD 

.770U 

.0000 


170.640 

561 i 2 

2116.2200 

.0009 

.7894 

520.000 

.000 

4532.0 

2142. OCOO 

. 7700 

.0000 


17D.33P 

563ir 

2116.2200 

.0013 

• 7891 

540. ODU 

.POO 

4q0. C 

3276. crOO 

. 7700 

.upnn 


171.310 

565^3 

2116*2200 

.L‘018 

*7887 

s6o.orc 

.000 

4152. r 

1835. LOOP 

. 7700 

• UODO 


172.700 

568.5 

2116*2200 

.0025 

. 788 3 

580.000 

.r'OO 

4t;56.C 

1 563. ?r^nO 

. 7700 

.unuo 


AK/RP\t f“ *8 
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FLOW field ICONTINUEOI 

2 R P- Ft- F7 F5 F*i F5 F6 


580. OOU 

174, r^O 

572.7 

2116,2200 

i0033 

.7877 

bro.ooo 

.000 

4 168.C 

1932,0000 

. 77 OO 

.0000 


175.150 

581.2 

2116.2200 

• 0050 

.7865 

620.000 

.non 

4U4. C 

1335. GOOD 

.7700 

• 0000 


174i3i0“ 

- “588T6 — 2X18^2200 — 

~i0063 ” 

;7856 

640.000 

.000 

«»120*C 

2355.0000 

.7700 

• 0000 


173.430 

597.6 

2116.2200 

.0079 

.7845 

680.000 

.TOO 

•*[,20.0 

2901.0000 

.7700 

.0000 


217.180 

559,1 

2116i220Q 

iC003 

• 78 98 

*<««WARNIM 6 - TEMP PANgF 

OF 316. 

0 To 3130,2 IS OUTSIDE THE 

RANGE OF COEFFICIENT TABLES 


BAND MODEL 

PARAMETERS 

Calculation TiMe 

^‘'•aOOOO SECONDS 

BROADENING 

COEFFICIENTS 

FOR UNIDENTIFIED SPEcIfS g* G3 G4 

G5 G6 



,u5 ,u8 ^06 ,uS 

• 05 *01 


H2Dl(Gi 

ABSORPTION COEFFICIENTS ~ 




-- ■ 

WAVE NO 

. TEMPERATURES 

(DEG H» 


k 



i l/CH) 

3 rc. 

600. 

lOpo* 

I 5 CO. - 

20005 

2500. 

30005 

12 EU. 

1.476-03 

2.232-02 

4.463-U2 

5.496-02 

5.304-02 

l♦•934-02 

4.649-U2 

16U0. 

5.463 -qI 

2.367-01 

1.252^01 

6.934-C2 

4.406*02 

35244-02 

2.956-02 

ZuuU. 

2.935-02 

3.065-02 

2. 529-^2 

2.036-02 

1.698-02 

1.413-T2 

1.237-02 

24C0i 

l. 2 sn~a<» 

1,591-04 

4iU13-L4 

7. 840-04 ■ 

1?774-03 

-2i486-M3 

2i522-03 

28 no. 

1.143-U4 

2.160-04 

4.792-04 

1.331-03 

2.647-03 

4.025-03 

4.154-03 

32t;0^ 

3 . 510 -OJ 

4 . 297-03 

1. 126-02 

1.695-02 

2.122*02 

2i084-q2 

li727-02 

36C0* 

1.2<ll-01 

1.357-01 

1.023-01 

5.646-C2 

3.386-02 

2,366-02 

1.8*21-02 

4CC0. 

1.867-rjl 

6.958-02 

3. 931-02 

2 . 310-02 

1. 574-02 

1.228-02 

1.090-02 

44C0* 

2. 144-05 

5.886-05 

2. 33i-r4 

4.487-04 

8.348-04 

1.315-03 

1.993-03 

48C0« 

5.262-aS 

1.644-04 

6iC71«“C4 

15777-03' 

35t37-OS 

3i862-03 

3i958-03 

5200 , 

1.937-02 

1.377-02 

1.244-02 

8.892-03 

6.527-03 

5.187-03 

3.930-03 

560n. 

2.127-C2 

8.697-03 

6.071-Q3 

4.477-03 

35679*03 

3.311-C3 

25831-03 


1.495-05 

2.230-05 

6.891-05 

1.731-Q4 

4.158-04 

8.396-04 

9.445-04 

64C0« 

4,118-05 

1.462-04 

3.187-04 

6.157-04 

15246-03 

1.858-03 

1.674-03 


1.3Dt-n3 

1. 965-03 

3.261-03 

3.333-03 

3.169-J3 

3. 108-03 

2.918-03 

77C0* 

1.724-02 

1.8 12-02 

1.095-02 

4.9$6^D3 

3. 26 4*03 

25595*0^ 

2.192-03 

7600. 

3.293-1*3 

5.191-04 

8 , 334-04 

7.129-04 

6.241-04 

5.565-04 

6. 153-04 

anoo. 

3.34C-U6 

6.010-06 

2.121-L5 

7.211-05 

1.910^04 

25497-04 

3.377-04 

84LC. 

l.925-n5 

3 .107-115 

9.499-rS 

2.174-04 

3.799-04 

3.896-04 

4,434-0<i 

6800. 

1.484-05 

1.396-05 

6.584-^5 

1.982-04 

3.498-04 

3. 179-04 

3.054-04 

92t0. 

2.092-06 

1.627-06 

3.826-F6 

1.713-oS 

8.846-05 

1.006-04 

1.473-U4 

9600, 

o.dcc 


U.UCO 

0.000 

15076-04 

1. 152-04 

1 5 196-04 

H20K6) 

FINE 

STRUCTURE PARAMETERS 

l/cT 



WAVE NO 

I. TEMPERATURES 

(DEG K| 




U/CM, 

310. 

6Q0. 

lOnu. 

1«;00. 

2000. 

2500. 

3000. 

12C0. 

2.739-rjl 

4.1S1-01 

1 .026*C0 

3.058 400 

i,c^5*ai 

3.999411 

1.265*02 


LONG FORM k AND 1/d OUTPUT 
(Parameters have been 
averaged over 400 cm”^ 
intervals and the output Is 
for the center of the interval. 
Integration interval is 
1000-9800 cm"*) 


SAMPLE PROBLEM 1 


RTR 014-9 



WAVE NO. TEMPER*TU»ES (DEG K) 

g/CM> 3CD. 600. IDOO. 1500. 2D00. 2500. 3000. 

1600. 1 . 461 -ni 3.855-01 7.975.01 j, 997*00 9 . 701+00 1.115+01 2. 963+01 

2t,nn. 1.161-01 2.992-01 6.709-01 1.861+00 9.803+Qff t.3Bi+Ql tt.DlS+Dl 

2900. 3.061-01 6.083-01 1.521+00 9.123+00 9.153.00 9.058+01 1.939+02 

2bCi3. 1.632+00 2.379+00 3.166+00 9.636*00 1.099+01 3.850+01 l.De9*Q2 

32po. 8 . 93 f ,-01 1.823+00 2.263+UO 2.689+00 5.221+po 1.352+01 2.558+01 

3600. 1.009*00 l , 92n+nn 1 . 8 . 1+00 2.650+00 3.865+00 7.887+00 1.369*p t 

9000. 7.882-01 1.170.00 1.512+00 2.29j+0Q 5 . 325+00 1.151*01 1.999+01 

■ MOn. 9.908-01 1.103 + 00 1.756 + 00 9.899 + 00 2.852+01 7.5n»*01 I. 7 OO + O2 

9800. 2.771-01 l.OSl-pl 3.131-01 1.230+00 5.808*00 2.222*01 8.313*01 

S 200.. 1.163-01 8.919-02 2.616-01 1.139 + pr, 9. 662 + 00 1 . 978* 01 1.082 + 02 

5610. 1.288-01 2.239-01 3.879-01 1.557+00 8.065+00 3.536+j.l 2.799*02 

MOO. 3.897-01 1.239-01 8.185-p, 5.961+00 7.212+01 2.787*62 1.396+03 

69o!J. 3.930-01 8.677-02 5.201-01 3.382+00 9.519+oj 1.709+02 9.368+02 

6800. 1.529-01 5.9n7-02 3.075-01 1.917+00 1.517*01 6.210+01 2.731+p; 

7200. 1. 257-01 3.229-02 1.952-01 1.965*00 7.937*00 3.556*01 1.300*02 

76nn. ?. 339-01 1.097-01 1.D82+DC 6.028+00 2.B92+q1 1.230+02 9.115+02 

scon. 9.992-Dl l,o65*00 2 . 76 O+UO 7.888*00 1.878+Sl 3.709*01 6.068+qI 

89aa. 2.292-01 9.850-01 1.279+00 3.5|n+0n 8.29 a*00 1.693*01 2.727+Dl 

880„. 1.161-01 2.585-01 6.799-01 1.931*00 9.577*00 R.pjfc+OO l.WaJ+Ol 

9200. 1.803-01 9.2Q8-nl 1.039+00 3.108*00 8.686*00 1.765+01 2.913+n i 

9503. 0.000 0.000 0.000 fi.poo 1.995*01 3.826 + 01 6.273*01 
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-CM-SR AVG cOMB OPT PATH lBLoCK 

VG TRANS CM-ATMISTP> 

52ER0 

SMIN 

DY= -4 .595-OlCM 

, 02= 2.606*0OCM 

0. 

0. 

OY:: -l*323^naCM 

, 02= 2.292*0ncH 

0* 

t' • 

OYr -M*280-nlCM 

, 02= 3.662*0rCM 

0. 

0. 

OYr -I,26l^nOCM 

, 02= 3.465*0nCH 

0. 

c. 

DY= -2*n2S*n0CM 

, 02= 3.080*OnCH 

0. 

Oi 

^75-03 ,964 

1.309*nu 0 

912. 

930. 

0Y= -l.219*nOCH 

f 02= 4,549*orCM 

0. 

0. 

OY- -1.990*O0CM 

, 02r 4.268-»O0CM 

0. 

0. 

0Y= -2.70t+anCM 

, 02= 3.858*OnCM 

0. 

c. 

H8-01 *753 

2.079*01 0 

726. 

942. 

.805 

1.797*01 0 

744. 

942. 

165-OX ,944 

2 . 649*^0 0 

912. 

946 • 

0Y= -3.275^00CM 

, 02= 4.677*onCM 

0. 

0. 

ipl^nO .477 

4.309*01 0 

612. 

96q. 

94*00 .488 

4.222*01 0 

6X2. 

96C. 

85*00 .754 

2.890*01 Q 

6 66. 

966. 

07= -3.136*00CM 

f 02= 5.899*0ncH 

0. 

0* 

0Y= -3.927*00CM 

, 02= 5*405*00CM 

0. 

0. 

27*ro .393 

6.126*01 0 

5 28 . 

984. 

41*n0 *370 

6.917*01 n 

5 22e 

9 8 4. 

28*00 .567 

4.938*nX 0 

546. 

984. 

53*00 .8 34 

2.924*01 0 

630* 

990. 

DY= •3,8jrj*00CM 

, 02= 6.599*0nCH 

0. 

0. 

0Y= -4.550*00CM 

, 02= 6,lx2*0nCM 

0. 

0. 

22*00 .403 

6.837*01 0 

4 68. 

IOO 8 • 

53*00 .385 

7.866*01 0 

462. 

1C08. 

5X*00 .596 

5.183*01 0 

480. 

966. 

65*00 .874 

1.600*01 0 

5 64. 

766. 

0Y= -4i26l*00CM 

» 02= 7.380*0PCM 

0. 

n« 

Dy= -5.089*n0CM 

, 02= 6.836*0DCM 

0* 

0. 

43*00 .356 

8^140401 0 

4 26. 

1032. 

49*00 .400 

8.206*01 . 0 

4 14. 

1032. 

47*00 .498 

4.787*01 0 

420. 

822. 

59*00 ,812 

2.582*01 0 

4 56. 

738. 

Dyt -4.071+00CM 

« 02= 8i454*orCH 

Pi 

0« 

DY= -4,892*0UCM 

, 02- 8.U06400CM 

n. 

0« 

3Y= -5i665*00CH 

» 02= 7.479*0PCM 

0* 

0. 

16*00 .355 

8.932*01 0 

390. 

1062. 

13*00 .404 

8.760*01 0 

3 78 i 

1062. 

51*00 .527 

3.895*01 0 

3 84. 

702. 

^8*?0 i859 

1.789*01 0 

4 14 • 

630. 

)Y= -4.425*OOCM 

f 02= 9.188*0rCM 

n. 

V* 

)Y= -$.317*0aCM 

• 02= 8.701*0PCM 

0. 

0. 

)Y= -6.157*nOCM 

♦ 02= 8.129*0nCM 

0. 

0. 

19*00 .395 

8.163*01 0 

360« 

990. 

I5*n0 *.475 

7.150*01 0 

346. 

924. 

)2*^0 -536 

4.180*01 0 

348. 

624. 


THESE HAVE BEEN MODIFIED FROM 
INPUT LIMITS OF SZERO AND SMAX 
BY SUBROUTINE SLIMIT. UNITS ARE 
THE SAME AS INPUT 

M l1 
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tmeta 

PHI 

SHAPE rCUX-EATTSrSTJ’»Cf1 


RAD"*WATTS 7 SQ*CM»SR AVG comb OPT 

Path lplcck 

SZERO 

SMIl 

DEGREES 

factor delta sum 


LOS AVG TRANS CM 

-ATMiSTP) 



•♦ 6.00 

17.500 

3.P7I-P2 2.16-03 li 

15^00 

7in96-01 6i597^00 .784 2 

.098+01 0 

366 . 

588* 

46.00 

22.500 

3,i68-P2 1.83-06 !• 

IS^OO 

6.025-04 6.484>n0 .969 1 

.444400 0 

4 32. 

486. 

4 6.00 

27.500 

LINE OF SIGHT MISSED THE 

GASi 

0 X = 

I.059»ClCMj oF= -5iC62*00CM. 

DZr 

9.72440PCM 

0 . 

C. 

46.00 

32.5Q0 

LINE OF SIGHT HISSED THE 

GAS. 

ox = 

i.059+0lCMt OYr -5.890+00CM* 


9.246 + OnCf^ 

0 * 

0 . 

46.00 

37.503 

" LINE OF SIGHT HlSSEtr^Hr”G«TJ-t>X^ 

"1 .059f0lCMr'0T 6 ,674^0CCMi 

02 r 

8*697+0nCM 

Oi 

0 • 

60.00 

2.500 

3.263-02 4.67-02 1. 

19^00 

1.558+01 6.315 ^nai .445 6 

.960+01 0 

3 36. 

852. 

5 C.no 

T.sro 

3 . 359-02 9 . 66-02 li 

29^00 

3i222+0l 6i71?^no .400 7 

.477401 0 

3 24 . 

798. 

50. DO 

12.500 

3.454-02 6.23-02 !• 

35 + 00 

2.n77«0I 6.93S«0O .452 4 

.891401 0 

3 24. 

558. 

5o*00 

17. 5cn 

3.550-^02 5i93»04 li 

3S+PC 

li977-01 6*833-00 .832 1 

.747401 0 

342. 

522. 

50.10 

22.500 

LINE OF SIGHT MISSED THE 

GAS. 

OXr 

9. 796*aCCH, OV: -4 .u 68 ♦t'OCM, 

02 = 

1.079401cm 

0 . 

0 . 

50.no 

27.5C0 

line of sight HISSEO^HE*- 

G^ASv 

Oyr- 

9i796»0nCM; OY = -Si391-00CM* 

02 r 

li03640lCM 

0 . 

0 * 

50.00 

32.SC0 

LINE OF SIGHT MISSED TmE 

GAS. 

OXr 

9.796*0nCM, Dy- -6.273+OOCM, 

02 = 

9.846+OOCM 

0 . 

n. 

50.00 

37*500 

LINE OF SIGHT MISSED TME 

CASi 

DX = 

9i796»0PCMi 0Y= -7ilo7*OOCM, 

02 = 

9.262-OPCM 

0 . 

0. 

54.00 

2.222 

LOS BLOCKED BY BODY FLAP 





3601 

0 . 

D. 

54. OC 

6.667 

LOS BLOCKED BY BODY FLAP 





3601 

0. 

c. 

54.00 

11 . Ill 

LOS BLOCKED BY BODY FLAP 





3601 

0 . 

0. 

54.00 

15.556 

LOS BLOCKED BT BaOY“FrAP~ 


— 

- 


3601 

Oi 

Oi 

54.00 

20.CU3 

LOS BLOCKED BY BODY FLAP 





3601 

n. 

D. 

54.00 

24 .444 

LOS BLOCKED BY BODY FLAP 





3601 

0 . 

n. 

54.00 

28.889 

line of sight missed the 

GAS. 

ox = 

8*958^0nCMt 0^- -$.957+00CHf 

02 = 

I.U6Q401CM 

0 . 

a. 

54.00 

33^333 

LINE OF SIGHT MISSED THE 

GASi 

DXr 

8.958+orCM» OYr -6i775+00CMi 

oZ= 

U03Q401CH 

0* 

a. 

54.00 

37.778 

lInE of sight MISSED THE 

GAS. 

DXr 

8.958-aOCM, OYr _7,553-00CM, 

0 2= 

9*745+00CH 

0. 

a* 

58.no 

2.222 

LOS blocked BY BOOY~ruAp“ 

— 

— 

. 


36P1 

Oi 

0* 

58.00 

6.667 

LOS BLOCKED BY BODY FLAP 





36^1 

0. 

0 . 

58.00 

11.111 

LOS blocked by body flap 





3601 

Gi 

0 . 

58.00 

15.556 

LOS BLOCKED BY BODY FLAP 





3601 

0. 

0 . 

58.00 

20.000 

LOS blocked by body FLAP 





3601 

0. 

D. 

58.00 

24.444 

LOS BLOCKED BY BODY FLAP 





3601 

0. 

0. 

SSiOO 

28.889 

LINE OF SIGHT MISSEO-TME 

OAS»' 

-DX2 

8r076-00CMi qY: -6«244*00Cm» 

02 = 

1.1 3 2- 01 CM 

n« 

0. 

58. OD 

33.333 

LINE OF SIGHT MISSED THE 

GAS. 

DX- 

8,g76-onCM, OY: -7.102-OOCM, 


1.08U-01C,i, 

n* 

u. 

58. OC 

37.778 

LiNt Of" sight missed THE 

GASi 

OXr 

8.U76-00CM, OYr «7.9j7«aCCM* 

D2= 

1*022+01CM 

0 . 

0. 

6 2 • 00 

2.222 

LOS blocked by BODY F LAf> 





3601 " 

0 . 

0. 

62.00 

6.667 

LOS BLOCKED BY BODY FLAP 





3601 

0. 

0. 

62.00 

11.111 

LOS BLOCKED BY BODY FLAP 





3601 

0 . 

0. 

62.00 

15.556 

LOS BLOCKED BY BODY“FLAP- 

— 


- - 


3601 

G» 

D« 

62. JO 

20.000 

LOS BLOCKED BY BODY FLAP 





3601 

0. 

C. 

6 2.00 

24.444 

LOS BLOCKED BY BODY FLAP 





3601 

0. 

G. 

62 .no 

28.689 

LOS BLOCKED BY BODY FLAP 





3601 

0. 

c. 

62 . on 

33.333 

LINE OF SIGHT MISSED THE 

GAS. 

OXr 

7.155-mCM, oY:: -7.394-OOCih. 

02: 

1.124401CM 

0« 

0. 

6 2.nc 

37. 778 

LINE OF SIGHT MISSED THE 

GAS. 

DXr 

7.155-OPCM, 0Y= -e.243-00CM, 

qZ= 

1.064401CM 

0* 

D. 

66.no 

2 .ono 

LOS BLOCKED BY BODY FLAP 





3601 

0. 

c« 

66.00 

s.aro 

LOS blocked by body Flap 





3601 

0* 

0. 

66.00 

10.CC3 

LOS BLOCKED BY BODY FLAP 





3601 

0. 

c. 

66.00 

I4.cnn 

LOS BLOCKED BY BODY FLAP 





3601 

0. 

G. 

66. CD 

I8.CG0 

LOS BLOCKED BY BODY FLAP 





3601 

0. 

C. 

66. OC 

22. JGQ 

LOS BLOCKED BY BODY FLAP 





3601 

a. 

0. 

66.00 

26. COO 

LOS BLOCKED BY BODY FLAP 





3601 

0. 

0. 

66. ID 

3o.uno 

LOS BLOCKED BY BODY FLAP 





3601 

c. 

G. 

66.00 

34.000 

LINE OF SIGHT MISSED THE 

GAS* 

OXr 

6«199‘D'^CH, qY: -7.785-OOCm, 

0 2: 

1.154401cm 

0 . 

0. 

66.00 

38.CC3 

line of SIGHT MISSED THE 

GAS* 

CX = 

6,199-OOCM, 0Y= -8.572-POCM, 

0Z= 

I.Q97401Cm 

0 * 

C. 
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THETA PHI 

OEGKEES 

SHAPE 

EACTOR 

ELUX=WArTS7S0=Xr« RA(r=WATTS/SQ*C«»SR AyO COMB OPT PATH 

DELTA SUM LOS AV6 TRANS Cm-ATM(STPI 

LB LOCK 

5 ZERO 

SMIN 

7n*.io 

2 .rno 

LOS 

BLOCKED 

BY ORBITER LOWER EUSELA6E 

1200 ^ 

C* 

L . 

To*OQ 

6. ano 

LOS 

BLOCKED 

BY ORBITER LOWER EUSELAGE 

120o 

0. 

0. 

7n*D0 

10 . coo 

LOS 

BLOCKED 

BY ORBITER LOWER FUSELAGE 

1200 

0. 

u. 

TC*00 

14.00C 

LOS 

BLOCKED 

BY ORBITER LOWER FUSELAGE 

1200 

0. 

c* 

00 

18. COD 

LOS 

BLOCKED 

BY BODY rxAP^ - 

3601 

0. 

04 

7n.oa 

22 . 0 no 

LOS 

BLOCKED 

BY BODY FLAP 

3601 

0. 

c . 

70.00 

26 .cro 

LOS 

BLOCKED 

BY BODY FLAP 

36C1 

0. 

G. 

7’".0CJ 

3U.U0U 

LOS 

BLOCKED 

BY BODY FLAP 

3601 

a. 

c. 

7n,ao 

34. coo 

LINE 

OF SIGHT HISSED the GAS. DX= 5.212»0nCM, oT= -8.P08*0aCHt 02= 1.187»01CH 

0. 

0. 

70.00 

38.000 

LINE 

OF SIGHT MISSED THE GAS. OX: 5.2j2*anCH, DY: -e.817«00CM, q^= I.129*01Cm 

n* 

c. 


NUMEPIcAl PLUMP Sh^PL PACTqP* iD^SSO PLUME RADIANCE r U2l«i3>ni W A ITS / SO-CM TER 

NUMERICAL TOTAL SHAPE EacTOR: *09843 OVERALL RADIANCE^ 4*379l*00 WA TTS/ SO-CM-STER 

EXACT TOTAL SHAPE EACTOR t .10052 TOTAL FLUX = 1. 3541400 WATTS/SO^CM 


> 

CO 


CALCULATION TIME 


12 SECONDS 



NOT ACTUAL TIME 



SAMPLE PROBLEM 1 



TRANSMISSIVITY 




WAVE 

1/CM 

WAVE 

IRP AOIANCE 


OPTICAL 


L E NG T H 
MICRON 

W A t TS / 
CM 


DEPTH" 

- FOR LAST LINE-OF-SIGHT WHICH MISSED 
THE GAS ON THIS PROBLEM. 

1200* 

8.333 

3. 191-05 

l.CUOti 

r.ooo 

160C* 

6.250 

9. 102-05 

x.noon 

O.OOCi 


2C00. 

2400, 

b.OOO 

4.167 

1.787-04 

1.675-04 

1 .oono 

0.000 


1 .pooa 

L . 00 G 


280r. 

3.571 

2.5^4-C4 

i.naoo 

0.000 


3200. 

3.125 

2.992-04 

l.OODU 

0.000 

- • • - - 

3600. 

2.778 

1.534-04 

l.ntjoc 

o.uoo 


4C0C. 

2*600 

2. 198^04 

i.rjoo 

0.000 


4400. 

2.273 

1.676-04 

1.0000 

O.GOU 


4600. 

2 . 0 8 3 

3. 142*' w4 

1 .0000 

oicno 

^ . 

5200. 

1.923 

3.rOQ-C4 

l.OOOO 

0.000 


S60Q. 

1.786 

2.636-C4 

i.rooo 

c.ooo 


6COO. 

1.667 

1.090-04 

i.nono 

0.000 


6 400. 

1 .563 

1.908-04 

i.naoo 

c.uoc 


6600. 

7200. 

1.471 

1.389 

2. 352-C4 
1.826^04 

i.nno 

i.rooc 

u.oou 


' OitlOC 


7600. 

1.316 

6.436-05 

i.nooL 

0.000 


8C00. 

1 .25U 

2.9nx-05 

i.noao 

D.UOO 


8400. 

1 . 1 9 n 

3.943-05 

1 .oDon 

o.oon 


8800. 

1.136 

2.890-G5 

i.raon 

CsOOO 


9200. 

960r. 

1.C87 

1.042 

9.605-C6 

9.030^06 

1.0^00 

i.oooc 

o.uoo 


C i CO 0 


.P^CTRAL 

OUTPUT 

time 

.orucu 

SECONDS 



SAMPLE PROBLEM 1 
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PEAK RADIANC E 

(9x3.22 < N < lOxNMAX)" 


B 8 B B 

B B B B B B 

B B B BB B B B 
BB9BBB6B 

JB B INDICATES BLOCKING* 

X(?) 6 86 B 

6 5 0 

7 8 30 * 

_ 

75 0 
7 71 0 
77 10 
55 0 

0 OP 

0 


W INTO PAPER 



OUTPUT POLAR PLOT OF 
RADIANCE IN SUBROUTINE 
VIEW 



<180 DEG» 
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RADIATION problem I DEN T I F I C A 1 1 ON - SAMPLE PROBLEM lA. HE*T TRANSFER FOR A 0.04 SCALE MODEL 


ft 


FLc“ EIELO identification - - - - SSME 

GAS 

FLIGHT NOZZLE SEA LEVEL PLUME - TKE - 5/4/78 
IDENTIFICATION - LH2/L02/AIR 



rl 

ft 

tun 

. LOCATION 

OF EXIT 

CENTER 

DIRECTION 

COSINES 

OF AXIS 

M 

NO. 

X 

_y __ 

Z 

X 

Y 

Z _ 

19 

1 

63.84 

• 00 

19.45 

.96126 

•OOQOO 

.27564 

09 

2 

64.91 . 

-1.12 

14,80 _ 

c9J481 _ 

.00000 

.17365 


3 

64.91 

2.12 

14.80 

.98481 

.00000 

.17365 

9» 

FLOW 

FIELD CONIrAL LIMtTS •• 

ImTERCEpT.: 

L -2A5g.0i.qiL- 

-?L.pP^r_ 


If 



English units - length 

IN INCHES - TEMPERATURE IN DEGREES 

RANKINE - 

PRESSURE IN PSFA 

** 


LIMUS or INTEGRATION 

vartasi FS 

initial VALUE 

FJNAL VALUE 

INCREMENT _ _ _ . __ ^ _ 

** 



Z-INOEX 

1 

55 





THETA 

8.00 

74.00 

4. nano 

«r 



PHI 

• 00 

40, DO 

4.0000 

fC 



NU 

»1000. 

9800. 

400. 

ir 



S 

• 00 

72*n0 

• 2400 

« 



TEMP STEP StZE 



2L00.U000 _ _ __ 

«c 



TANGLE «DEG) 



4.0000 

•c 

cc 


POINT location 

X Y 

57.38 .00 

2 

2.54 


ALL GEOMETRY IS SCALED BY 0.04 

B 






FROM SAMPLE PROBLEM 1 

ec 


COSINES OF U 

-.00000 -.00000 

l.OQOOO 



a 


V 

.00000 -1.00000 

•00000 



K 


w 

_UP.oqgg_ _».D.atiflfl. 

•00000 



tz 


■CONSJIT.uenT? N2t.l 5? 


J. 


a 

a 


COORDINATE DATA INPUT FOR THE POINT OF INTEREST AND PLUMfS 

^ ICS NO. X Y g CHI PSI OH^GA 

pojNT OF Interest 

0 57.38 iSfi 2.54 180. 00 90.00 .00 

PLUME NO. 1 

5 iUfi tJlfl tills tXm diQ tiUL 

PLUME NO. 2 

6 !. 0 g s_Qfl 6.26 iilQ *JIQ sM. 

PLUME NO. 3 

7 jPD iOg 6.28 .00 .00 .00 


•t 

tt 


RTR 014-9 







BLOCKING SURFACES 

NO« title TYPf ALF BHIN BHAk GHIN GNAX 

ICS RX RY R2 CHS PSS OHS 

6201 ENG 2 BELL 7 1,7138-01 i.bg70*00 5,1538+00 7o 360.0 

6 0,0000 0,0000 1 .1?62 + D0 ,0 j_0 ,0 

1200 ORBITER LOMER cUSFLAGE 3 0.0000 O.OOOD ».9628+0i -6,0 6.n 

D 1.1080+01 O.Oqqo 1.1068+01 90.0 .0 .0 

36ll BODY FLAP 3 0.0000 i.97<l9+01 2.3733 + 61 -12.6 1276 

0 8.9151+01 0,0000 1. 3569 + 01 270,0 ,0 -6.3 


INTERMEDIATE COORDINATE SYSTEM FOR BLOCKING SURFACES 

ICs Rxl ill RZI CHC PSC OHC 

5 5.780U*0l 0.0000 1. 7720+01 180.0 79.0 ,0 

6 5.8727+01 -2. 1200*00 1.3706+01 180.0 80.0 .0 

1 5T577T^TJI ;^.1206+o0 i.37p6+ol IBO.O eoTo TcT 


NO FLOWFIELD OUTPUT 


SCALE FACTORS RSrAL TSCAL PSCAL FI F2 F3 F9 ^ 

• ui4uuO l«OriODO X«D0OQO ieOOQOO 1 «OqODQ 

wi'm TODEi ThnjHmLni 


•e 


H 


feROAOENING COEFFICIENTS FOR UNIDENTIFIED SPEcI^S $1 G2 G3 G9 G5 G6 

.05 .08 .06 .05 .05 ^1 


H20KG) ABSORPTION COEFFICIENTS 

wave no. TEHPERAtURES JDEG KJ 

(1/CM, 300, 600. lOno. I5OO. 2000. 2500. 3000. 

l200. 1.47B-(J3 2.212-02 9.963-02 5.996-02 5, 3o9-02 9,939-02 9.699-02 

9600. 0.000 Q.QOO 0.000 0.000 1.076-09 1.152-09 I.i96-Qg 

H20 KGJ FINE structure PARAmETERS 

T7XVT“R0r rrHPERATuRR ' TDEg' k j 

(I/CM) 300. 6pn. 1000. 1500. 2000. 2500. 3000. 

12D0. 2. 73^-01 9.i51-0l 1.026+00 3,n58+0D 1.065+01 3.999+01 1.265+02 

96oq. g.OOD O.DOO 0.000 0.000 1.995*ni 3.826+01 6«g73*01 


SHORT FORM OF PARAMETER OUTPUT 








r-- "i ' „ -I, 




u 


RTR 014-9 


S/SQ“CH»SR AVG COMB OPT PATH LBLOCK 

SifRO 

SMIN 

AVG TRANS CM-ATH(STP» 



JICH, OY: "li838-n2CM 

, DZ= 1.042-0JCH 

D. 

0. 

JICH, Ov= -5.293-02CH 

, DZ= 9.167-02CM 

c. 

a. 

)1CH» or: -Ii712-02CH 

, D2= 1*465-01CM 

n. 

n. 

llCM, t)Y: -5.04«-f?2CM 

, 02= 1.386-OlCH 

D, 

0 . 

IlCM; DY: «3ilO'»-n2CM 

i 0?= 1*232-01Ch 

r* 

Q* 

( 1.447,-04 ,997 

5.235-02 0 

36* 

37 . 

JlCHi OY: -4.876-n2CM 

, 02= 1.820-OlCM 

C. 

0 . 

JICM, OY: -7.961-02CH 

. 02: 1.707-01CH 

0. 

0 . 

JlCMi DY: -Ii080-01CM 

, 02: 1*543-01CM 

0. 

a. 

L 2.752-02 .951 

8.318-01 f) 

29. 

38. 

! 3 ;4 18-02 i969 

7*107-01 0 

3Ci 

OP 

*■ 

! 3.125-02 .995 

1.060-01 0 

36. 

38. 

JlCMi OY; -Ii310-01CM 

, 02= U871-D1CM 

0. 

D. 

J 2.461-01 .905 

1.723*00 0 

24. 

38. 

) 4il81-01 *908 

1*689*00 0 

24 i 

38. 

1 4.029-01 .953 

1.156*00 0 

21* 

39. 

nc«5 DY: sj.jSS-OlCM 

i 02= 2i36a-QlcM 

Oi 

Oi 

31CM, OY: -1.571-OlCM 

, 02= 2 . 162 -OlCM 

n. 

n. 

3 5. 815-01 .871 

2.451*00 0 

21 . 

39# 

) 8.073-01 .858 

2.767*00 0 

2i« 

39. 

3 8.332-01 .910 

1.975*00 0 

22. 

39. 

Z 7.990-01 .966 

1.170*00 0 

25. 

40* 

3ic«; OY:- -1 i524-giCH 

* 02= 2!640-01CM 

0* 

Oi 

31CM, DY; -1.820-01CH 

, 02= 2.445-q1CH 

P. 

0. 

3 8.851-Cl .866 

2.735*00 D 

19. 

40. 

3 1.034.00 .853 

3.145*00 0 

16. 

40. 

3 1*035.00 .915 

2,073*00 Q 

19. 

39. 

5 9.997-rjj .979 

6.4C2-01 0 

23. 

31. 

31CMi OY: -Ii704-01CH 

j 02= 2*952-01CH 

Oi 

0* 

3 ICM, Py: -2.n35-01CM 

, 02= 2.734-OlCH 

0. 

0. 

3 1*070.00 .845 

3.255*00 0 

\1* 

m* 

3 1.161.00 .854 

3.283*no D 

XI* 

4 1. 

g 1.201*00 *904 

1*915*00 0 

XI* 

33. 

2 1.170. 00 ,964 

1.033*00 0 

18. 

30. 

31CM} OY: -1*628-01CM 

, 02= 3.381-OlCH 

Oi 

P* 

31CH, DY: -1.957-niCH 

, 02: 3.202-01CM 

0. 

0. 

3tCH* OV: -2.266-OlCH 

, 02 : 2.992-OlCH 

0. 

0. 

0 1.192*00 .838 

3.573*00 0 

16. 

42. 

0 1.267*00 .8^3 

3.504*00 0 

15« 

42* 

0 1.291*00 .917 

1 . 558*00 0 

IS. 

28. 

3 1 *262*00 *975 

7. 156-01 0 

i?* 

25. 

DKM, DY: -1.77n-01CH 

, 02= 3.675-OlCH 

0. 

0* 

OlCM, OY: -2.127-OlCH 

, 02= 3.481-OlCH 

0. 

0. 

OICH, DY: -2.463-OICM 

, D2= 3.252-OlCM 

0. 

C. 

D 1 .238 *00 .8 55 

3.265*00 0 

14. 

4 0. 

0 1.263*00 .831 

2.860*00 0 

14. 

37. 

0 1.275*00 .914 

1*672*00 C 

14i 

25* 


OBLEM lA 


RTR 014-9 



ThETA 

PHI 

SHAPE FLUX -(TAT nrS^J's-Cfr ' 

'17A0-’‘WATTS7S0=CM^SB AVG C0M8 OPT 

Path lplock 

SZEPD 

SMIN 


DEGREES 

factor delta sum 

LOS AVG 

TRANS CM 

.ATMlSTP) 




46,00 

17.500 

3.n71-02 4.20-P4 2.18-01 

Ii378*01 1.255*no .964 8 

.393- 

oi 0 

J 5 * 

24. 


46. nn 

22.500 

3.158-02 2.44-07 2.18-01 

0.OUS-O5 1.234+00 *’97 5 

.776- 

H2 0 

17. 

19. 


46.00 

27.50Q 

LINE OF SIGHT HISSED THE GAS. DXr 

4,235-OlCMi DY: 

-2.025-01CM, 

OZ: 

3.89C-01CM 

0. 

0. 


4 6.00 

32.sna 

LINE OF SIGHT MISSED THE GAS, DX= 

4.235-OlCM, OY:: 

-2.356-OlCMt 

DZ: 

3.698-OlCM 

0. 

0. 


4S.00 

37.500 

LINE OF SIGHT HIS SED“rRE”GA-S-; ' OX- 

“-4;235-eiCMf OYr 

-2;669'‘D1CM, 

DZ: 

3.479-OlCM 

Oi 

Hi 


50.00 

2.500 

3. 263-02 9.35-03 2.28-Dl 

3.117400 a .204*no .874 2 

.784+00 0 

13. 

34 . 


50.00 

7*500 

3.359-02 1.54-02 2i43*0l 

5;I244DQ 1.265+00 .858 2 

•991^n0 0 

13. 

32. 


50.no 

l2.5Pn 

3.«,54“C2 1.07-02 2.54-01 

3.550+00 1.300+00 .892 1 

.956+00 0 

13. 

22 . 


50.00 

17.5P0 

3.55C-02 1.17-04 2i5*»-01 

3.905-02 1.28I+P0 .972 6 

.988- 

01 0 

14. 

21. 


.00 

22.5r0 

LINE OF SIGHT MISSED THE GAS. DXr 

3.918-OlCM, 0Y= 

-1.787-OlCM, 

DZ: 

4.314-OlCM 

0. 

0. 


50. DO 

27.5C0 

LINE OF sight MISSED' THE-GASrOX=: 

- 3. 918+‘ClCMi DY: 

— 2 *156 *0 1C M t 

OZ: 

4;142-01cM 

0. 

Oi 


50. 00 

32.500 

line OF SIGHT missed THE GAS. Ox= 

3.9l8-DlCMt 0Y= 

-2.509-.0ICM, 

DZ: 

3 . 938 -QlCM 

0. 

0. 


5G.00 

37.500 

LINE OF SIGHT HISSeD ThE GAS. DX= 

3.918-OlCH, Oy: 

-2.843-OiCM, 

OZ: 

3.705-01CM 

0 . 

0 . 


54.']fJ 

2.222 

LOS BLOCKED BY BODY FLAP 




3601 

0 . 

c. 


54,00 

6,667 

LOS BLOCKED BY BODY FLAP 




36D1 

0. 

0. 


54,00 

11. m 

LOS BLOCKED BY BODY FLAP 




3601 

0. 

0 . 


54,00 

15.556 

LOS BLOCKED BY BODY TfAP - 

..... 



36P1 

n. 

Oi 


54.00 

2C.uro 

LOS BLOCKED BY BODY FLAP 




3601 

0 . 

G. 


54, ?0 

24 ,444 

LOS BLOCKED BY BOOT FLAP 




3601 

Oe 

0, 


54, 00 

23,889 

LINE OF SIGHT HISSED THE CAS. DX: 

3.583-OlCM, ovr 

-2.383-niCMf 

OZ: 

4.318-OlCM 

o« 

D. 


54,00 

33.333 

LINE OF SIGHT HISSED THE GAS. 0X= 

3.583-OTCH, oy: 

-2.710-niCM, 

d 2 = 

4.120-OICm 

0* 

C. 


54.00 

37.778 

LInF OF sight missed the GAS. OXz 

3.583-OICH, oy: 

•+3.Q21-OICH, 

D^: 

3.898-OlCM 

0. 

0 , 


sa.oo 

2.222 

LOS blocked by boot FrAp 

... 



3601 

Oi 

Ci 


se.oo 

6,667 

LOS BLOCKED BY BODY FLAP 




3601 

n. 

a. 


5 9.00 

n*ui 

LOS BLOCKED BY BOUY FLAP 




3601 

Oi 

D* 


58.00 

15.556 

LOS BLOCKED BY BODY FLAP 




3601 

0 . 

0, 


58.00 

20.000 

LOS BLOCKED BY BODY FLAP 




3601 

0 . 

Oi 


58.00 

24 .444 

LOS BLOCKED BY BODY FLAP 




3601 

0 . 

0. 


58.00 

28.889 

LINE OF SIGHT HIS SE0“-THE-GAS F 0X = 

“ 3;230~01CKf oY= 

-2.498-DlCHf 

DZ: 

4i5Z6-QlCM 

Oi 

(Ji 


58.00 

33.333 

LINE OF SIGHT HISSED THE GAS. DX= 

3,230-OICM, oy: 

-2.841-OlCH, 

dZ= 

4.319-QICh 

0* 

0. 


58.00 

37.773 

line OP sight HISSED THE GAS. DXr 

3.230-01CM, DY: 

^3,167-aiCM, 

02 = 

4.086-01CM 

n. 

Oi 


62*00 

2.222 

LOS blocked by BODY FLAp 




3601 

0. 

0. 


62.00 

6.667 

LOS BLOCKED BY BODY FLAP 




36P1 

0. 

0. 


62.00 

11. m 

LOS BLOCKED BY BODY FLAP 




3601 

0 . 

D. 


62. OC 

15.556 

LOS BLOCKED BY BOOT FrAP " ’■ 




36P1 

Oi 

Oi 


6 2 . or 

20.000 

LOS BLOCKED BY BODY FLAP 




3601 

0 * 

0. 


62.00 

24.444 

LOS BLOCKED BY BODY FLAP 




3601 

Oi 

0. 


6 2.00 

25.889 

LOS BLOCKED BY BODY FLAP 




3601 

0 . 

0. 


6 2.00 

33.333 

LINE OF SIGHT MISSED THE GAS; 0X= 

2.862-OlCM; qY : 

-2.958-OICm. 

DZ: 

4,497-OICM 

0 . 

D. 


6 2.00 

37.778 

LINE OF SIGHT MISSED THE GAS. 0X= 

2,862-OlCM, oy: 

-3 .297-OlCM, 


4,254-OICm 

n* 

C. 


66.00 

2 .crc 

LOS BLOCKED BY BODY FLAP ‘ 




3601 

0. 

c. 


66.00 

6*000 

LOS blocked by body FtAp 




3601 

0 , 

Cl. 


66 . CO 

10*DCo 

LOS BLOCKED BY BODY FLAP 




36P1 

0. 

r. 

yo 

66.00 

14.000 

LOS blocked by body FLAP 




3601 

a. 

0 . 

— 1 

66.00 

is.ora 

LOS BLOCKED BY BODY FLAP 




3601 ' 

0. 

0 . 

pp 

66.no 

22.000 

LOS BLOCKLD BY BODY FLAP 




36D1 

0 . 

0 , 

0 

66. nn 

26.000 

LOS blocked by body flap 




3601 

Oi 

Oi 

p-y 

66.U0 

30.000 

LOS BLOCKED BY BODY FLAP 




3601 

c. 

0 . 

1 

66.0C 

34.000 

LINE OF SIGHT MISSED THE GAS. DX= 

2.479-01CM» oy: 

-3.114-OICm, 

DZ: 

4.617-aiCM 

n. 

c. 

CO 

66.00 

38 . OPO 

LINE OF SIGHT MISSED THE GAS. DX- 

2.479-OlCM, oy: 

-3.429-OlCM, 

0^' 

4.388-OICh 

0* 

0, 



SAMPLE PROBLEM lA 



theta 

PHI 

SHAPE 

FLUX-aArT$7SD-CFr T?A0-WATTS/S0»CM-$R AvG COMB 

OPT PATH 

LPLOCK 

SZERO 

SMil 

DEGREES 

FACTOR 

DELTA SUM LOS 

AVg TRANS CM- 

A TMIstP) 




v^.no 

2.00Q 

LOS 

BLOCKED 

BY ORBITER LOWER FUSELAGE 



1200 

0« 

n. 

7o«OQ 

6.000 

LOS 

BLOCKED 

BY ORBITER LOWER FUSELAGE 



12Do 

0. 

0. 

ToOO 

iQ.noo 

LOS 

BLOCKED 

BY OPbITER LOWER FUSELAGE 



i20n 

0. 

fu 

7o*OD 

14.C0Q 

LOS 

BLOCKED 

BY ORBITER LOWER FUSELAGE 



I 20(5 

0. 

c. 

Tjj.oa 

18. 'ICO 

LOS 

BLOCKED 

BY B^OOr-FXAP 

— 


3601 

r* 

0. 

70. no 

22.UC0 

LOS 

BLOCKED 

BY BODY FLAP 



3601 

0. 

0. 

70.00 

26.0C0 

LOS 

BLOCKED 

BY BODY FLAP 



3601 

0* 

0* 

70.00 

30. COO 

LOS 

BLOCKED 

BY BODY FLAP 



3601 

0. 

D. 

7^. no 

34.cna 

LINE 

OF SIGHT MISSED THE GASi 0X= 2iD85 

-OICM# qT= -3.2O3-01Ch» 

D2- 4*749 

-OICM 

0. 

D. 

70.00 

38. noo 

L INE 

OF SIGHT MISSED THE GAS* OXr 2.^85 

-OICM, OY: -3.527-DlCM, 

□2- 4.514 

-01 Cm 

0* 

L • 


NUMERICAL PLUME SHAPE FACTOP^^ *03550 PLUME PAOIANCE - 2*2771400 WA TTS / SQ -CM-S TER 

numerical total shape FAC^OP- .098*43 OVERALL RAOIANCEr 8.2119-01 W A TT 5/ SQ-CM-S T ER 

EXACT TOTAL SHAPE FACTOR r .10052 TOTAL FLUX r 2i5393-0l WATTS/SO-CM 


calculation TIME 11 SECONDS 


cn 

o 


70 
— * 
70 

0 
*-* 

1 

LO 


SAMPLE PROBLEM lA 



EACH symbol represents A RADIATION STEP OF .567148 WATTS/SO-CM. 

THETA 5CALE“GOES“TD 9TT70OT OEG^RtES 

<0 OEG> 


BBSS 
B B B B B 6 

8 B B 8B B B B 
B B B B e BB B 

p -g .„ .g .p ~ g 

5 9 6 BB B 

« 740 

9 7 30 ♦ 

9 4 0 

74 0 

-?■ 81 — 0 “ 

89 20 
55 0 
0 GO 

0 


cn 

K-* 


DEGI « 




« (9n DEG) 




♦ 


* 




♦ 

1180 DEG) 


FINISH 


SAMPLE PROBLEM lA 


RTR 014-9 



7CW 


I 


A 


-- tr 

2 . 

-Tr 

6. 

R. 

10 . 

Ttr 

12 . 

TTr 

't5T 

16. 

“TTr 

le. 

-t^n- 

20 . 

2Tr 

22 . 

■«r 

24. 

■rsT 

26. 

"2TT 

2P. 

•m- 

30. 

-mr 

32. 

-^Tr 

34. 


3 irstrfT;J “t3ASl? A07 tJ9S7l7tr?0" i 

3REWIN0 GASPAP 

L._ .. 

— srRrc" Tprs. ■ " ■ 

SASG,T TPFt. .F/l/POS/S j 

1 

uASu , T 3TF77T7STT0 f 

aASG,T i*,F///?500 ; 

> SAME AS SAMPLE PROBLEM 1 

8 A s G , T"B 7 Tr r / strtr i 

dASG,T 10.tr///2500 \ 


'^OPTfT^irsPAcr ^ 

8FREE 6ASRAD 

l.._ _ 

"tnnn — gaspao 1 

SAMPLE PROBLEM ?. ANALYSIS 

— -r— 12 — 5 - 1' 

*♦0. SO. in. 

— rrtr. - T70. 2rrr 

2000. 4000. 2S,. 

rrr r r 

OF BAND MODEL PARAMETERS ON A SINGLE LOS. 

— 0 1 TJ- ■ 13 a~~ U- 

m-fYr 


J — 

H20i(6> C102(G) CIOHG) N2CG) 02(6> 


tm^D •TlTAtr 


0 6?. 0. 0. 0. 0* 0. 

tr — 3 i"j-5 cn cn xn vn — 

0 -31.5 0. 0. 0. 0. 0. 

1 


-1 

-rrcr — rzv — i — 1-4 r" i s i t3 2 ^ i 4^ — rs2 — t3T“"1 ?t3 — t5T^ 

HpO, CO 2 , and CO BAND MODEL PARAMETER 

^AOD COGDC / CARD IMAGE DATA ELEMENTS FROM INPUT TAPE 

IL 3 BLANK CARDS FOR BAND MODEL INPUT 

ONE BLANK CARD TO TERMINATE RUN 

3xcr -rrtrr EXECUTE PLOT PROGRAM 

8FIN 


INPUT 


LISTING FOR SAMPLE PROBLEM 2 
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LBH CTjTAN first STAGE) 2qK LEvEt PLUME <8/30/79). 

BETA H20/C0/c02/N2/02 

NOCONr 5 H 2 OKG) 


<jn 

<x> 


.00 

5*00 

10.00 

15*00 

20.00 

25*00 

30.00 
35i00- 
<t0*00 
45*00 

-50.00 

55.00 

50.00 

-70.0a 

®o*oo 

90«00 

ioo«oo 

11 O.QO 
120.00 
130*00 
130.10 
1^0*00 

150.00 
150 .QO 
ITq.OO 

180.00 


1 

• 00 
- 1 
• 00 

1 

• DO 

1 

• 00 

1 - 

.00 

1 

.00 

I 

.00 

1 

.00 

1 

• 00 
1 

.00 

1 

• 00 

- 1 - 

• 00 
1 

• 00 

4- 


ClOl <S) 
30 
30 
30 
30 
30 
30 
30 


ICASE= 1 

ISYHSG= 0 HANOsf .00 - GTGOs .00 AP= S.OOO^ot BP=~ 3^000-01 

C102C6’ N2<G) 02<6) 


30 

30 

30 


30 

30 

30 


PARTIAL ECHO PRINT OF FLOWFIELD INPUT 
IN SUBROUTINE FFPREP 


• 00 
1 

30 

15 

X 

• 00 

~ 1 ■ 

• DO 

30 

30 

l3 

1" 

— i” 

• 00 

30 

l2 

1 



• 00 

. 1 

30 

lO 

• 00 
t 

30 

5 

A 

• 00 

30 

2 

1 


— 

• 00 

30 

3 

1 



• 00 

1 

30 

1 

• 00 

30 

lo 

1 



• 00 

29 

l3 

1 



• 00 

J 

30 

if 

• 00 

26 

1 








SAMPLE PROBLEM 2 
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>• 

cn 

4 i» 



l^OiOO 

200*00 

220*00 

2«lOiOO 

260*00 

280.00 

3o0*00 

320*00 

3 < i 0.00 

340*10 

343*75 

360«00 

380*00 

398*30 

420*00 

440.00 

460*00 

462*00 

- 502*80 

520*70 

541*60 

- — - 562*50 

580*30 
600*50 
- 621*40 

641.70 
662*40 
682*60 
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2 AXISYHHETRIC pLtJHEs USING CASE 


1 irlTH EXIT-tTRtrNE-COCATHJNS -ANO DIRE CTl OKS INttlCjfTEo RUK TfSTr OSCJOfl 


RADIATION PROBLEM IDENTIFICATION - SAMPtE PRoBLEF 2. ANALVS^^F BAND MODEL PARAMETERS ON A SINGLE LOS. 

FLOW FIELD IDENTIFICATION .... L6M ITITAN FIRST STAgE> 20k e^VEl RLUME IB/30/79). 

- - -GAS lOENTTFICATlOK -=^-«20/COFC02/K2/02 


ENG. LOCATION OF EXIT CENTER DIRECTION COSINES OF AXIS 

NO. X-- If 2 -X- V 2 

1 31.50 .00 .Do .00000 .00000 1.00000 

2 - -3U50 iOO iOO - iOOOOO .00000 1.00000 


FLOW FIELD CONlcAL LIMjTS - InTERCErTs S.OOO*01 SLOPE: 3.000-01 


English units - length in inches - temperature- in degrees- rankinc -- pressure in psfa - 


LIMITS OF integration 


variables initial value final value 

2-lNDEX 1 LAST 

THETA 40.00 50.00 

PHI - - 170.00 170.00 

NU 2000. 4000. 

S - - - —.00 600.00 

TEMP STEP S12E 


INCREMENT 


20*.SSSS-}-SET to give a single LINE-OF-SIGHT at 45,180 
2sl 

--3r0000 - - - 

SQ.OOOO 


X Y 2 

POINT LOCATION * 63«0C - »00 ^00 — 

COSINES OF U liOOOOO vOOOOO #00000 

V #00000 1.00000 #00000 

- N -#0000o VOODOO 1 #00000 


CONSTITUENTS H2oit65 CI02IG) — cioitcr N2iei “ ~ “ 02 (6 > — NO PARTICULAR ORDER REQUIRED 


COORDINATE DATA INPUT FOR THE POINT OF INTEREST AND PLUMeS 

ICS NO. X Y CHI PSI OMgSA 

POjNT OF Ij^TEREST 

0 S3. 00 *00 -,0Q -,Qo-- .00 .00 

PLUHE NO. 1 

^ 0 3 Is SO -9 00 §00 ^00 sOO ^00 ~ 

PLUHE NO. 2 

0 •^Jl.SO “ *00 “~«Q0 »00 .00 .00 " 


NO FLOW FIELD SCALE FACTORS WERE USED 
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71 .0<I39 .3592 .0000 

00 sOOOO --^7900 *2100 

7i .0 <i39 .3592 .0000 

00 *0000 t 7890 *2099 

71 .Q<l39 .3592 .0000 

03 *0002 *7863 *2092 

71 .0939 .3592 .0000 

00 ^0000 *7900 *2100 

71 .0939 .3592 .0000 

Oa *0000 *7899 *2099 

71 .0939 .3592 .0000 

01 — *0000 *7897 .2099 

71 .0939 .3592 .0000 

02 .0001 — *7890 .2095 

71 .0939 .3592 .0000 

90 *0000 *79o0 *2100 

71 .0939 .3592 .0000 

DO *0000 - *7900 *2100 

71 .0939 .3592 .0000 

30 .0060 *7899 *2100 

71 .0939 .3592 .0000 

30 .0000 *7899 .2099 

71 .0939 .3592 .0000 

31 *0000 *7896— *2098 

71 .0939 .3592 .0000 

32 *0002 .7891 .2095 

71 .0939 .3592 .0000 

30 .flOoo .7900 .210& 


71 .0939 .3592 .0000 

10 .0000 *7899 .2100 


n .0939 .3592 .0000 

J1 .0000 .7897 .2098 

rj ,0939 .3592 .0000 

12 .0000 *7892 .2095 


'1 .0939 .3592 .0000 

10 -.0000 *7900 ,2100 

1 —.0939 .3592 .0000 


10 

• 0000 

.7900 

• 2100 

0 

.0990 

.3590 

• 0000 

10- 

— .oodo — 

-.7899 

-.2100 

0 

^0670 

.3390 

•0000 

10 

•0000 

.7899 

• 2100 

10 

•06S0 

• 3«»00 

• 0000 



F2 

F3 

F9 

F5 

0001 

• 0000 

• 7898 

.2099 

0620 

i06i|D 

♦ 3930 

50000 

0001 

.0000 

• 7 895 

.2097 

0690 

.0630 

• 3950 

•0000 

0003 

^VOOOO 

• 789^ 

52099 


0650 •0620 iOOOO 

OOOQ .0000 .7900 *2100 


0660 •0610 •3960 *0000 

0000 ^0000 i7900 i2l00 

q700 *0580 *3990 *0000 

0000 ^0000- ,7900 i2100 “ “ 

0790 •0590 •3520 *0000 

rOODO •000(^—5 7900 ¥tlOO — 

0750 •0530 *3530 .0000 

rOOOt — 50000 57698 i2099 

.0750 *0530 •SSBO *0000 

VVfn •0000~~5t7896 52097 

.0760 .0520 .3590 «0000 

ruOOO '•OOOO #7900 »2100 

.0790 *0500 «35S0 #0000 

r000(7— •0000 57900- - »2100 

F0B2O »098(7 »3970 - ^0000 

*0000 #0000 *7900 •2100 


»0670 .0990 «3590 *0000 

*0000- *0000 ^7699- .2099 

iO870 50990 ^3590 *0000 

•0001 *0000 *7897 *2098 


.0570 

50001” 

sOBSO ~ 
• 0001 

• 0670 

.3390 

• OOOO 

.... 

• 0000 

“50670 
• QDOO 

■ .7897 

i3390 
• 7897 

52098 — 

50000 - 

• 2098 

• 0580 

• 0670 

• 3390 

•OOOO 


• 0000 

iOOOO 

- .7900 

- *2100 — 

■' “ 

50580 

i0630 

5 3390 

iOOOO — - 

- - -■ 70 

■! 

• 0000 

.ODOD 

• 7900 

.2lOO 

70 

^0600- 

»C6S0 

^3900 

- .0000 _ . . 


• 0000 

•oooo 

.7900 

•2100 


*0600 

*O6S0 

*3900 

• QOOq 

- - 

t 

• 0000 

.0000 

• 7900 

• 2100 

1 

KO 

.0600 

>06SO 

53*^00 

*0000 

■' - - - 



F2 

F3 

r«» 

FS 

• OOOO 

•OOOO 

.7899 

.2100 

•0610 

50690 

iJOlO 

.OOOq 

• OOOO 

• OOOO 

.7899 

.2099 


i09SB *0356 *3597 ,0000 

•0001 •OOOO *7896 *2099 

— 5 0955 - 50356 i3598 »0000 

•0001 *0000 .7896 •2096 

”•0955 •0356 *3598 #QOOo 

•0002 *0000 47899 •2096 

1^0955 ^0355 s 36GI »0000 

•0000 .0000 #7900 *2100 

-50^56 -^0351 i3616 •0001 

•0000 *0000 .7900 .2100 

»0959 •0338 • 3656 .0002 

•0000 4OOQO *7900 #2100 

.0969 50306 ^iT3762 sOOlO 

•0000 .0000 47900 42100 

50961 ^0226 59089 - *0066 

•0000 •OOoo •7900 •2100 

^0950^ 50190 59 259 — 5OIU 

•OOOO #0000 •7899 •2X00 

•6932 40156 5 9990 *0169 

•OOOO *0000 *7899 *2100 

40909 50«2 49611 •0231 

•OOOO #0000 .7899 *2099 

— 40896 5^0118 -^9705 #0267 

•0001 •OOOO #7898 *2099 

-69 W6E~0F-€06r FI €I 


t 5 5OO00O SECONDS 


— 59 65 66 

•05 *05 .01 


00 . 

38-02 

79-02 
10 - 02 “ 
36-03 
99-03 
97-03 
15-03 
23-03 
90-03 



iOOw 3000.=- — - - 

?87-02 2.748-02 

rS2-02 2^803-02 * ” 

>29-02 2.730-02 

10 9-02 -2 # 8 39-0 2 — 

►76-02 2.475-02 

t45-02 2«4U-02 - - 

>59-02 2.193-02 

100-02 2 # 020-02 - 

>12-02 1.893-02 

59 1-02 1 r829-02 — — — 

559-02 1.911-02 
591-02 I .938-02 - 
581-02 1.820-02 

i 07-02 1 r3 65 -G2 

343-02 1.028-02 

l68-O2^“^^828-03 

256*02 9*919-03 

223— OB—l-rl-l^j— C2~ 

387-02 1.238-02 

583-02^1.401-0 2- 

824-02 1.611-02 

097-02-2.793-02 

206-02 1.902-02 

206-02-3^84-^02-^ — 

097-02 1.729-02 

835-02-1.465-02 


500, 3000.^ 

SIO^OO 2.440>0i 

I30^00-2533o^0l 

010401 2.610401 

1 5 04 03 8 1 04^01 

240401 3.4804Q1 

100401-3.580401 

680^01 6.540401 

440401 ^ 8.900401 

900401 1.1704Q2 

09C4o1"1^2S0402 

300401 1.330402 

590401 6.410402 

790401 1.51040? 

080401 1.610.02 
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Oi 3000 i- - - - ■ — 

0^01 1.710^02 
0^01 li«20402 
O^Ol 1*9<I0^02 

©♦0l~""2»O7O^O2 

0+01 2.210402 

O^Ot 2.350^02 - 

0401 2*51p402 

O^Olr-e'^670402 - - : 

0401 2.840402 

04^01 "2 • 6Oq402 

0401 2.370402 

0401: 2^170402 - - ~ - - 

0401 1.9804Q2 

040 h i 10402 — - r — - 

04q1 1.650402 

04 0 l~t7r51O402 r 

0401 1.370402 

0401 1»2504q2 

0401 1.150402 

G4G1 i^»05G4S2 - - - - - 

0401 9.6004q1 

D4Olr-0r9604Ol 

Q4Q1 a.00Q40l 

0^01—7*300401 

0401 6.670401 

D4|j1 -^.^9040 1 — 

04dl 5«560401 

0401-5*080401 

O 4 QI 4.64o4q| 

0401 - 4*240401 - - 

0401 3.870+01 

(J40l~3 .-5-30401- 

0401 3.230401 

(?4Q 1 - 2« 9sO401 — 

0*01 2.690401 

0401-2*460401 ^ 

34q1 2.660401 

3400-2*290401-' - - — 

3400 1.730401 

3*0 0~ 1V7 Op407 

3400 1.370401 

34 00 1»360401 - - 

3400 1.29q4q1 

)400-1»070401-- - - - - 

3400 9.500400 
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00* " 3000^ ' 

a5-03 3.185-Q3 

131- 03 2^839-03 
01*03 2*839-03 
^ 37-03 2 ^ 502-03 
63-03 2.430-03 

132- 03 2^357-03 
;88-Q3 2.293-03 
(81-03 2i284-03 
184-03 2.266-03 
^98-03 2^248-03 
113-03 2.257-03 
r04-03 2i266-03 
»60-03 2.284-03 
^38-03 2.302-03 
>60-03 2.348-03 
^ 55 - 03-2 5 ^ 393 - 03 - 
Ul-03 2.457-03 
?28-03 2.548-03 
>23-03 2.684-03 
108-03 2^894-03” 
?21-03 3.121-03 
J4^C3-3t840-03 - 
>64-03 3.795-03 
573-03 4 i 377-03- 
r07-03 4.550-03 
>25-o3 55005-03 
788-03 5.478-03 
B53-03 65051-03 
B69-03 6*688-03 
273-03 75289-03 
539-03 8.053-03 
708-0385 8 27-03 
190-02 9.646-03 

1 90-02^-15 065-02 
452-02 1.174-02 
551-02~1.292-02- 
835-02 1.410-02 
933-02 15547-02 
020-02 1.702-02 
250-02 15838-02 ' 
490-02 2.002-02 
774-02 2.193-02 
872-02 2.384-02 
090-02 2.548-02 
221-02 2.666-02 
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30* 3000# 
SO^GO 1«010«01 
50*^00 JsG 30^01 - 
»Q«00 I.320 «q1 
?0^00-|»590^0I 
lo^oo i.eso^oi 

JO^OO 1»9 qq 4 o 1 — 

ro^oo 9.m0400 

>0^00 l^OOO+Ol 
70400 8»770400 
J0400 3r.290^0l- 

r040I 9.97040I 

ro40G^^ •010401 
iO^OO 1.580401 
> 04^0 l“2^r4 504^ o ^— 
t0400 1*360401 
^0400-1^^60-^01 
>0400 l*16Q4gl 
^04 eo 8^6 90 40 G ~ 
30*00 9.300400 

t&400-1^08O40j 

S0400 9.990400 
rtJ400-l,^16O4Gl 


lOr 3000 w 

19-02 6.626-02 
h6-02 9 . 7 65*02 
>8-02 1.S29-02 
^7-0t 2.r^29*01~ 
tS-Dl 3.336-01 
M-Gl 9^707*01 - 
>9-01 6.930-01 
^3-01^-7.969-01 
12400 9.772-01 
!8400 - 1.089400 
2400 1.2124Q0 
►9*00 1.127400 
i34Q0 1.089400 
9 400 -8.639-03 
>9400 6.591-01 
;9-01 2.986-01 
9-01 3.213-02 
10 0.000 
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2^00i 3000* 

•000 0*000 

• 000- OiQGO ' ' 

•000 0*000 

• 000 OiOOO - 

•000 0 * 00 q 

6 000 - OiOOO — — 

•000 0*000 

iOOO 0*000 -- 

•000 0*000 

liOOO OiOOO — 

•000 0*000 

fiOOO OiOOO - — - 

1*000 0*000 

UOOO OiOOO - — - 

1*000 0*000 

1^000 OiOOO 

1*000 0*000 

iiOOO 0*000 — 

1*000 0*00o 

IsOOO- “-OiOOO 

1*000 0*000 

^i 000 ChiOCO- — 

1*000 1*061*03 

fiOOO— — lr33a-83 

1*000 1.8I»9*03 

UOGO 2c350*03 - 

f. 636-04 1*954-03 

^097-03 4*012-03 

1*629-03 5*014-03 

r • 120-03 6*702-0 3 - - - 

!*984-03 7*751-03 

^ . 2 1 7-03 9.705-03 

>•788-03 1*161-02 

r.842-03^~lr*34G-G2 

1*022-02 1*622-02 

U3i7-O2-i*604-B2- — 

I. 679-02 2*246-02 

1*977-02 2*^236-02 - 

?. 344-02 2.441-02 

! .74 0-02 2. 7 31-0^— 

M04-02 3*005-02 

J . 1 B 8 -0 2 ~2 • 9 1 5 -0 2 - 

5.443-02 3.247-02 

5*699-02 3.589-02 
1.092-02 4*098-02 
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300G* 

.ZZH-aZ n* 201 - 0 ? 
. 972-02 9 *^ 75-02 - 

SH7-02 $*03S-02 

.028-^2- 5rf 19-0? 

087-02 9*995-02 

911-02 9*729-02 

223-02 3.S21-02 

751-02- 2 »039-02 

<|87-03 1.859-03 


900. 3000. 

eS9^03 8*928^03 

893403 8,2694G3 

55 5^0 3-5^ 39-»0 3 

199403 9.373403 

G50^03 ^559-403 

931403 2.906403 

69 ?4-03“- 

292403 1.713403 

2«4^2~lrr283^&3 

8314q2 9.0104Q2 
608 ^ 02 — 8 -» 3 - 354 D 2 ’^" — 
399402 9.110402 

805402-5*169^02 

192402 1.378402 

36^0 1 - 8 . 9 - 1 3 -40 1 

>105400 8.3384Q0 

287^08 -8*981--01 

DOG 0.000 

300 — O»GO0- 

300 0.000 

3 0 0 0*0 08 

300 O.QOO 

300 - — >- 0*000 

300 0.000 

)00 ___0^POO 

)00 0.000 

>00 0*000 

300 0.000 

>00 0.000 

100 0.000 

IQO 0.000 


AMPLE PROBLEM 
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3000 i- - “ “ 

»000 0*000 

»00q^ o»ooo — — - 

»000 0*000 

rOOO o«ooa 

.000 0*000 

iOOO — 0*000 - 

.000 0*000 

iOOO - OsOOO “ - - 

.000 0*000 

i^3«r^02 6r23«r’^02 

*23<I402 6.23«I402 

• 23«*'»02ft*23%»o 2 - 

• 23<l«02 6.234+02 

i lTT+^02" t •344'+03 

•350+03 2*239+03 

r317^03- 2^1 2^03 

• U9+03 1*580 + 03 

f 321^03 2 Tl 18^03 

.470+03 2.408+03 

i706»03 2^800+03 

*140+03 3.545+03 

rl^2+-o3 5* 2 44-+^ 

•213+03 8*924+03 

V?30+03H5 rT53-+03 

*579+03 8*863+03 

t16’1+^03' 6»6 24^+03 — — 

.245+03 5*226+03 

»252+Q3-2^3rlr+cr3 

i*571 + 02 8.472+02 

Ilf- 12 1 +Q 2 — 4-r564‘'*'02— 

1*116 + 02 2.998 + 02 

rw64^2+02' ^+367-^02 — ” - 

t.094+02 3*230+02 

^7'57 ♦ 02 — inr 1‘22^+G 3 — — 

1*530+03 2.595+03 

U419+03-4rll8-+a3 

>.263+03 3.942+03 
U318 +032^1 88+03 
r. 939+02 1.275+Q3 

> « 0 7 7+ 0 2 3 « i 2 3+02 

J. 606+01 5.372+01 
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C 101 | 6 | ABSORPTION COEFFICieNTS - 

UAve NO. TEMPERATURES (0E6 Kl 

U/CH» 3D0. 600* 12oOi 1800* 280O*”- 3000* SOOO. 

2000. 6.077-OR 2.«»*l5-02 8.392-02 8.912-02 7.876-02 5.896-02 2.793-02 

2025. -8.000-03-1.000-02- 5.000-02 1.188-01- 8»6<l3-02 -6.259-02 -2i789-02 - 

2050. 9.801-02 2.645-01 2.209-01 1.407-01 9.026-02 6.102-02 2.589-02 

2075 . 1*000-01 5*000-01 2.400-01 1.427-ol 8.294-02-5.419-022^466-02 

2100. 1.871*00 8.408-01 2.621-01 1.197-01 6.617-02 4.528-02 2.542-02 

2125. 2.198*00 5.000-01 1.587-01 7. 665-q2 5.226-02 4.22S-622.850-02 

2 I 5 Q. 1 . 099*00 2.866-01 1.112-01 0.114-02 6.624-02 5.617-02 3.569-02 

2175.-3.261*00 1.103*00 3.372-0 11.768-01 1.185-01 8.757-02-4.457-02 

2200. 1.340*00 1.001*00 4.204-01 2.286-qI 1.476-01 1.057-01 5.005-02 
2225. 1.270-01 3.963-01 3.251-01 2. 093-qI 1.432-01-1.045-01-4.924-02 - 
2250 . 2.128-03 6.365-02 1.567-01 1.388-qI 1.Q86-01 8.478-02 4.299-02 

2275 . 3.451-06 3.156-03 4.209-02 6.256-02-5*^162-02 5.465-02-3.220-02 

2300. 1.623-10 2.708-05 4.832-03 1.636-o2 2.347-02 2.544-02 t. 936-02 

2325 . 0 rOOO 9 r4 23-09-^*202-04-1*656-03-4*563-0 3-6*939-03-9*364-0 5 

2350. 0.000 0.000 1.725-08 7.456-q6 1.153-04 4.958-04 2.958-03 


NOTE THAT FINE STRUCTURE PARAMETERS (1/d) 
WERE NOT INPUT FOR CO. THIS IS THE ONLY 
CONSTITUENT FOR WHICH THIS CAN OCCUR WITH 
THE PRESENT CODING. 


SAMPLE PROBLEM 2 
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THETA PHI 

degrees 

45.00 180.000 


NUHERlCAt PLUHE SHjPE F*CT0R= " iO09T0 ptUHE RADIANCE = “r;9223'»OO HArTS/SO=CfP=STER 

NUMERICAL TOTAL SHAPE FACTORr .00970 OVERALL RADIAnCEr 1.9223«00 wATTS/SO-CH-STER 

EXACT TOTAL SHAPE FACTOR r iOOOGj TOTAL FLUX = 5i8S57-02 MATTS/SQ-CM - 


SHAPE 

FACTOR 

9.696-q3 


FLUX-UATTS/SQ-CM 
DELTA SUM 

5.86>02 5.86-02 


RAO-WATTS/SQ-cM-SR 
LOS AVG 
I.922«00 1.922*00 


AVG ComB OpT path 
TRANS cN“*TM«STP» 
.5j 4 i.eoi^oi 


LBLOCK S2ER0 SHIN 


294. 


CALCULATION TIME 4 SECONDS 




<^AMP1 P PRflRI PM P 
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— XC/X* XO/X* - -AVE-FiNE- STRUCTURE 

parameters 




^097- 

- 9.772-02^- 

-4^^36CHG2 

1.000 

•999 

9,131400 

3,391400 

• 231 

.099 

6.996^03- 

-9*,23-03 

,386 

• 057 

9.703-02 

1.299-02 

.999 

1.000 

7,189400 

2,780400 

»22r 

.092- 

7.293-03 

9^598-03“ 

.9 37 

- .07$ 

5.0e$-02 

lr919-02 

• 996 

.999 


2.793400 


*16 6 - * 02R 6 *7 9 l»o3 - R rin 7»03 

— *506 *109 5*790-52 — i-*535-02 

.995 .999 7.006*00 2.796*00 

- *l6l *023 6*760-03^ 9.969-03 


— *539 .119 3.965-02 1 .362*52 

.999 ,999 6.998*00 2.795*00 

— .161 .023 6.653-03 9*356-03 

— .555 .132 5*306-02 — 1-.539-G2 

.987 .995 6.789*00 2.767*00 

— *202 .036 - 6.678-03- -9,916-03 


LONG FORM OF 
SPECTRAL OUTPUT 


— .68'9- -*23^ ‘^T*968*02 — 2.39 8—02 

.967 .965 6.199*00 2.599*00 

— .220 — ^0^2 6*®53-03— 9.591-03 -- 

—*590 .230 6.303-02 — 1*870-02 - 

.920 .839 5.392*00 2.329*00 

— . 138 — *« 8 6*739-03 — 9*563-03 ^ 


— *793 *30«f 7*565-02 - 2 . 270-02 

.828 .571 9.511*00 2*009*00 

—.129 .015 6.729-03 -9.691-03 


— *788 — -.375 8,3* l-o2— 2,532-02 

.669 .303 3.209*00 1.957*00 

.137 -iOl’ 6.796-03- 9.828-03 

- ,897 ,998 9.261-02 2.805-02 

.995 .150 2.152*00 9.952-01 

- *381 — ^ io’l 6,926-03 — 9.952-03 — 

.867- .559 1,006-01 3*089-02 - 

.380 .088 1.608*00 7.575-01 

.288 .073 7.119-03 5.o73«03 
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«AVE have IRRADIAnCE 6 OPTiCAt COWSTITUENT X* X/X* ‘ XC/X* " XO/X* ■ A VE FINE STRUCTURE " 

NUMBER LENGTH W*TTS/ DEPTH IOENtIFICA- PARAMETERS 

1/CM MICRON CM TiON AC AO 

-STER 


2300. <*.3<I8 1.140-03 .0000 1.462*01 

l*001-Ol--H20ilS» -lrirl&-01— -*898 .893 .646 1*103-01 —3.402-02 

1.448*01 C102I6) 4.728*01 .306 .302 .056 j.i88*00 5.614-01 

3.975-02 C10H61 7.0‘'8-02 - *564 .542 *256 7*329-03 5,196-03 

2325. 4.301 8.912-04 .0000 l.lsZ^Ol 

8.453-02 H20i<G) — 9*180-02 *921 *917 *738 1. 213-01 3.772-02 

1.173*01 C102J6) 5.43o*01 *216 .214 .030 6.530-01 3.096-01 

6*021-03 ClOKei 6*461-03 *932 *907 *899 7*507-03 5.309-03 

2350. 4.255 9.333-04 .0002 8.317*00 

7.173- 02 HaOltgl 7*638-02 -*939 *935 *815 1*335-01 4*180-02 

8.246*00 C102I6I 5.512*01 .150 .149 .015 3.125-01 1.455-01 

2*975-05 CtOtIGJ 2*975-05 1*000 1*000 1*000 7*720-03 5*453-03 - 

2375. 4.211 1.2o8-o3 .0349 3.356*00 

5*801-02 H2011B1 ”6^068-02 *956 *952 *886 3*470-01 4^627-02 

3.298*00 C102«G» 3.910*01 .084 .084 .005 6.982-o2 3.128-02 

2400* 4.167 7.948-04 .s900 5*276-01 

5.237-02 H201(Gt 5.430-02 .964 .960 .918 1.606-01 5.099-02 . 

4*752-01 C102181 Iji866*00 *098 *097 *007 1*166-02 5*741-03 

2425. 4.124 9.797-05 .9575 4.342-02 

4*342-02 «20l IS) 9-*4S8-02 »^69 *951- 1.740-01 — St57 8-02 

2450. 4.082 8.496-05 .9637 3.696-02 

3.696-02 ~ H2eHS» 3*766-0? *982 .^76 —.970 1*924-01 -6.193-02 _ 

2475. 4.040 7.326-05 .9691 3.142-02 

3.142-02 H20H6) 3*183-02 *’87 *982 *982 2.100-01 — 6.821-02 - 

2500. 4.000 6.740-05 .9715 2.896-02 

2.896-02 H201IS) 2.926-02 .990 — .-984 ,987 — 2.253-01 — 7,433-02 - 

2525. 3.960 6.647-05 ,9719 2.853-02 

2.853-02 - H201IG) 2,878-02 *992 .986 .989 2.446-01 8,153-02 

2550. 3.922 6.644-05 ,9720 2.837-02 

2*837-02 H20H 6) ^2,858-02 *’93 - *987 .991 2.655-01 8,940-02- - 

2575. 3.883 6.779-05 .9716 2.877-Q2 

2*877-02 — H2eif S5 *992 *98 6 *990 2*535-01 — e,6L5-02 

2600. 3.846 6.931-05 .9711 2.931-02 

2,931-02 H201I SI - 2.956-02 *’91 *985 ,989 -2,432-01 -6*336-02 - - 

2625. 3.810 7.782-05 .9678 3.268-02 

3*268-02 H201I6) - 3*3o4-02 *989 - *983 *986 2*356-01 8*116-02 

2650. 3.774 8.203-05 .9662 3.443-02 

- - 3*443-02 -H201te) 3*987-02 *987 *981 *983 -2*220-01 — 7*751-02 

2675. 3.738 9.916-05 .9591 4.174-02 

4.174- 02 H201IG) 4.25o-02 *’62 *976 *974 2.120-01 7.476-02 

2700. 3.704 1.184-04 .9508 5.044-02 

5*044-02 H201(6) 5*172-02 *975 *’69 *959 1*999-01 7*147-02 

2725. 3.670 1.484-04 .9381 6.386-02 

6*386-02 H201I6)- 6»62o-02 *’65 »9S9 *930 1.881-01 6.822-02 

2750. 3.636 1.695-04 .9298 7.284-02 

7,284-02 H201I6) -7,612-0? *’57 *951 .907 -1.803-01 6.590-02 

2775. 3.604 1,878-04 .9220 8.120-02 

8.120-02 H201IG) 8.561-02 .949 .943 .880 1.710-01 6.315-02 
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WAVE WAVE tPPtOtAHCt G OPTICAL - CONSTITUCNT X* XXX* - XC/X* - XO/X* AVC FINE STPUCTORC - 

NUMBER LENGTH MATTS/ DEPTH lOENTlFICA- PARAMETERS 

1/CM MICRON CM TIOM - - AC AD • 

-STER 

2800. 3.571 2.279-0*1 .90*18 1*000-01 

— 1^000-01 M20Hfl» ^ls073-01 *932— *926 *623 1*617-01 6*052-02 

2825. 3.5*10 2.*t9S-0*l .8963 1.09*1-01 

1*09*1-01 H20HG1 1*187-01 *922 *916 *788 l.SAS-Ol 5*829-02 - 

2850. 3.509 2.758-0** .8866 1.203-01 

1*203-01 H20tt6» 1*286-01 iTjg - *929 *893 2*036-01 7*907-02 - - 

2875. 3.*«78 3.188-0** .8682 l.**13-0l 

- 1**»13-01— H20U8I ^1*535-01 *920 *913 *790 1*933-01 7*601-02 

2900. 3.***»8 3. ****1-09 .8577 1.535-01 

1*535-01- H20H6J ls753»or *87^ *869 *639 1*3**9-01 5.232-02 - 

2925. 3. **19 3.666-0** .8**77 1.652-01 

- - - — 1*652-01 - H201f6> 1*921-01 *860 “ *852 *595 1*275-01 ***^99**-02 - 

2950. 3.390 **.071-0** .8311 1.851-01 

— — — —1*851-01 — M201tO» 2*20«-Ol *8 38 *830 ^539 - — 1*221-01 — ** *817-02 

2975. 3.361 <*.369-0<* .8178 2.011-01 

- - 2*011-01 M20KOJ 2i*f6t=oi i«»7 TiSOB ***92 1* 160-0 t--**»6l6-02 

3000. 3.333 <*.855-0** .7984 2.251-01 

— ” 2*251-01 — H26iiSr 2*64/— 01 *791 *781 *439 1*115-01 — 4*463-02 

3025. 3.306 5.180-04 .7849 2.422-01 

- —2.422-01 — H20H6I 3.352-01 *768 *758 *400 1*067-0 3“4, 294-02 

3050. 3.279 6. 056-04 .7505 2.870-01 

2*870-01- H20H 87 3*965-01 *924 *713 *333 1.027-01 -4*372-02- 

3075. 3.252 5.992-04 .7531 2.835-01 

2.835-01 -H20H6I 3.956-ei *736 *706 .325 9.822-02 4,014-02 

3100. 3.226 7.095-04 .7077 3.457-01 

; ■ 3 .452— O'! — S . 26 1 “O’l ■’656 ,645 .256 9 -.3 93— 02 — 3.873—02 ' ' ■ " . . ~ 

5.117-04 C102(6I 5.117-04 1.000 1.000 .971 2.e84«00 1.693«00 

3i25, 3.200 7.867-04 .6746 3.936-03 — - - - 

3.921-01 H201IG) 6.461-01 .607 .597 .211 8.924-02 3.715-02 

1.516-09 - C102(6> 1,516-03 - 1.000 — 1,000 .998 4 ,634*00 - 2,756*0G- - - - 

3150. 3.175 8.502-04 .6520 4.277-01 

— — 4 .249—01 — — H20H6I 7^,414-— 81" *593 ,563* ,365 8,61‘4— 02 3,616—02 — 

2.825-03 C102(6) 2.825-03 1.000 1.000 .995 4.494*00 2.708*00 

3175. 3«lS0 9.019-04 .6288 4.640-01 - — - - 

4.618-01 H201I6) 8.665-01 .533 .524 .158 8.191-02 3.468-02 

2.187-03 C102I61 2,187-03 — 1.000 -1.000 .98q -4,571*00 -2,761*00 

3200. 3.125 9.661-04 .6026 5.065-01 

5.027-01 -^M203l«l l,029*oo »9«9 ^480 ,132 7,717-02 — 3,303-02 

3.748-03 C102(6I 3.748-03 1.000 1.000 .991 4.044*00 2.942*00 

3225* 3.101 1,105-03 .5340 6.274-01 - - - - 

6.215-01 H20tl6t 1.193*00 .521 .512 .156 1.057-01 4.690-02 

5.839-03 €102(61 5.839-03- 1.000 - 1.000 .999 5.280*00 3.225*00 - — ^ 

3250. 3.077 1.081-03 .5488 6.000*01 70 

— — 5.924-01“-4l201(6} 1.488*00 .398 i392- “ .088 6*755-02 — 2*974-02 

7.665-03 C102(GI 7.665-03 1.000 1.000 .998 6*199*00 3.781*00 

3275, 3.053 1.162-03 .5115 6.704-01 - 

6.557-01 H20K6) 1.820*00 .360 .355 .073 6.571-02 2,925-02 

1^474-02 C102(S» - 1,474-02 1.000-1.000 .999 7.170*00 4.444*00 


SAMPLE PROBLEM 2 


1014-9 




(*- — xc/x* “xc/x* xvr riNr sTHucTURE 

parameters 

- — - AC AO 


56 iSSI ?D6<t 6i1*t-0?- 2i87T-02 

DO 1.000 1.000 9.oS8«00 S.675«00 

61 .355 .076 8. 754-02 4.068-02 

DO -liOOO l.Ojjo" li254-»01 7i922«00 - 

23 *319- *062 8*074-02 3*782-02 

00 .999 l.Ooo 1.420«01 9.131«00 

14 .310 .059 8.481-02 4.014-02 

00 ^999 liOOO 1*523«01 9*969-*00 

85 .2 8 1 i Ol9 T* 73 4—02 — 3*685—0 2 

00 .999 1.000 1. 447*01 9.640*00 

75 .272 .047 7.908-02 3.798-02 

00 *996- -1*006 8*420*00 - 5*666*00 

47 *245 i038 6*855-02 3*315-02 — 

99 .990 .999 4.l60*oo 2.888*00 


37 .235 .036 6.549-02 3.180-02 

94 - *977 *994— 2*280*00- 1*619*00 

'32 *230- - *034 6*417-02 3.H6-02 

'80 .956 .977 1.420*oo 1.020*00 


134 .232 .o35 6.796-02 3.313-02 

i58 *932 *946 1*003*00 7*241-01 

r37 . 235 *036-^ 6 *7 4^1~G 2—5*297-02^ 

134 .908 .904 7.668.Q1 5.525-01 


»43 .240 .037 6.484-02 3.143-02 

,28 *902 *891 7.507-01 5.376-01- 

,44 .242 - *U37 5*7 54-02 - 2* 762-02 

>65 .939 .956 1.505*oo 1.072*00 

227 .225 .q32 5.215-02 2.476-02 

>87 .966 .986 2. 651*oo 2.032*00 

245 i24 3 iOjS $*849-02—2*875-02- 

»94 .976 .993 4.884*oq 3.534*00 

243 .241 .038 6.446-02 3.204-02 

99O .968 .989 4.291*00 3.145*00 
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W*Vf WAVE IRRAOIANCE 6 OPTICAt CONSTITWENT X* X/X* XC/X* - XD/X* - -jtWE FiNr- STRiJCTURC— 

NUMBER LENGTH WATTS/ OEpTH IDENTlFICA- PARAMETERS 

I/CM MICRON CM TIOn - - - AC AO ' - 

-STER 

3700« 2.703 1.693>03 .0945 2.360->00 

— - - - 1.149 *00 -- H20UGI 4^«T3*00 »?35 »2^3 *036 7*021-02 — ^3»4e«-02 

I. 211*00 C102(6| 1.240*00 .977 .951 .973 2.929*00 2.171*00 

3725* 2.685 1.599-03 .0988 2.315*00 - . 

1.084*00 h201IGI 4.4o3*00 .246 .244 .039 6.950-02 3.392-o2 

1.231*00 C102(6L 1.282*00 .960 .932 .950 2. 1 35*00 -1 .601.00 - - 

3750. 2.667 1.551-03 .1386 1.976*oo 

-- - 9*4 08-01 — H20K6) ly470*00 *271“ *268 *045 6*716-02 — 3*22tr02 

1.036*00 C102(G) 1.164*00 .889 .861 .816 8.330-01 6.392-01 

3775, 2.649 1.009-03 .5014 - 6.903-01 - 

6.044-01 H201(6I 1.573*00 .384 .377 .092 6.523-02 3.307-02 

8*590-02 C102< 6) 8*919-02 - *963 *933 *955 1*510-01 1*180-01 

3800. 2.632 9.452-04 .5143 6.650-01 

fr*650-07 — M201l€^l l*e69-*00 *-3S6 *350^ *079 6*511-02 — ^3*289-02 

3825. 2.614 9.688-04 .4767 7.409-01 

7*409-01- H201IG1 2r3o8*00- *321" *316 *065 6.410-02 3*206-02 

3850. 2.597 8.965-04 .4682 7.589-01 

7.589-01 — H201(6» 3.08o*00 *246 *244 *037 4*877-02 2.287-02 - 

3875. 2.581 1.061-03 .3923 9.358-01 

- “ ^ 9-.358— 0 1 — H20H G ) 3'.789*oo *247 *244 *04 0 6*009—02 — 1* 062—02 - 

^ 3900. 2.564 1.088-03 .4042 9.058-01 

9.058-01 — H201C61 — ^*7G^^^oo *245 *242 .040 5*752-02 — 2*997-02 

3925. 2.548 1.161-03 .3835 9.583-01 

9.583-01 H20K6I S.Vqo^OO *246 *243 .042 6.106-02 - 3.272-02 

3950. 2.532 1.214-03 .3745 9.822-01 

9.822-01 — H281IOI -3.932-*tl0- — *250— —.247 .043 6.366-02—5.471-02 

3975. 2.516 1.105-03 .4633 7.693-01 

7.693-01 H201t6) 2.744*00 *280 .276 .053 5.670-02 — 3.065-02 

4000. 2.500 1.035-03 .5079 6.774-01 

6.774-01 -H20HG1 2.096*00 .323 *318 .070 5.880-02 — 3*206-02 - 


SPECTRAL OUTPUT TIME 2.00000 SECONDS 


SAMPLE PROBLEM 2 


ie-Ho'aiy 



EACH SYWBOt ptPntStHJS t RAOTATlOK STEFOf S 192231 — WATTS/S9*CHr 

TH^TA scale goes To 9Q.000 DEGREES 

<Q 0E6) 


>> ■ 


y> 


1270 OEBI <• 


^ ^ -|-9q - 0£6| •- - 


* ^—SING LE LINE>OF-SIGH T 


# 




♦ 

U80 OEG> 


RTR 1014-9 



FINISH — ■ - 

NORMAL EXIT. CPU TINEt 24712 TOTAL SU^Si 57141 (MILLISECONDS > 


aXQT PLOT 


BEGINNING OF PLOT SECTION 


1 PLOT PROGRAM OUTPUT 

END OF PLOT SECTION. PLOT TIME 8.00000 SECONDS 

1 ■” ' - -■ - - .. 


>• “ 

^ NORMAL EXIT. CPU TIME: 4715 TOTAL SOPS: 9780 (MILLISECONDS I 
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ACDATA 

ACDATA 

BLOCKM 

BLOCKM 

BLOKIN 

BLOKIN 

FFPREP 

FFPREP 

FLOUT 

FLOUT 

FLOW 

FLOW 

FLOWAX 

FLOWAX 

FLOW IN 

FLOWIN 

FL0W3O 

FL0W3D 

MAIN 

MAIN 

NOGO 

NOGO 

PLANCK 

PLANCK 

PLUMPT 

PLUMPT 

RAD 

RAD 

SKIP 

SKIP 

SLG 

SL6MLT 


YNGH20 

SLIMIT 

SLIMIT 

SPCOUT 

SPCOUT 

VIEW 

VIEW 

YAPPRX 

YAPPRX 

YF 

YF 
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SUBROUTINE NAME: ACDATA 

DES C RIPT I Qtl 

ACDATA reads the band model parameter input (absorption coefficients, F, 
and line densities, 1/d). The data are scaled if specified by input and averaged 
over the spectral integration interval (DENU) if the spectral interval of the 
data (IDNU(I)) is smaller than the intergration interval (IDNU(I) < DENU). The 
band model data to be used in the problem is printed in either a short form 
(initial and final values in each table) or long form (the complete table). 

CALLING SEQUENCE 


CALL ACDATA 


COMMON BLOCKS AND EXTERNAL REFERENCES 

COMMON/ABCOEF/ 
COMMON/RADSUB/ 
EXP 
SQRT 

METHOD OF SOLUTION 


COMMON/CONTROL/ 

COMMON/CONID/ 

COMMON/FLOSCL/ 

COMMON/GOMTRY/ 

COMMON/BLOCKO/ 


Band model parameters are read in the format specified for cards 17 
through 19 (Appendix Al), and they are placed in the table COEF (vl,K) where 
the J index represents temperature and the K index represents wavenumber. The 
wavenumber interval for each constituent is specified by the initial and final 
values (NUL(I) and NUU(I)) and the spacing between entries is IDNU(I). The 
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subroutine assigns the index (K) for the initial values of absorption coefficient 
(ICOEF(I)) and the line density (IFIN(I)) to enable later location of the param- 
eters while maintaining a minimum size sequential table. 

If scale factors are used (Card 2, ISCALA=1) the band model parameters are 
multiplied by the nonzero scale factors input on Card 15. 

If the input spectral integration increment (DENU) Is larger than the input 
interval in the band model parameter table (IDNU (I)), the parameters are aver- 
aged as indicated in Eqs. 19 and 20 (Section 3.3). 

Note that there must be an integral number of IDNU(I) intervals in DENU, 
Selection of the format for printed output is made based on NABS (Card 2). 
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SUBROUTINE NAME: BLOCKT-I 

DESCRIPTION 

BLOCKM tests a Hne-of-sight vector for intersection with specified blocking 
(or shading) surfaces. If intersections occur, the upper limit of the vector 
(corresponding to the upper limit of integration) is adjusted to correspond to the 
closest surface intersection. 


CALLING SEQUENCE 


CALL BLOCKM 


COMMON BLOCKS AND EXTERNAL REFERENCES 


COMMON/CONTRL/ 

COMMON/CONID/ 

COMMON/ FLOSCL/ 

COMMON/GOMTRY/ 

COMMON/BLOCKO/ 

COMMON/BLOKOP/ 

METHOD OF SOLUTION 


COMMON/FLOWS/ 

COMMON/SLIMTC/ 

N060 

SQRT 

COS 


The upper limit of each line-of-sight ^LOS) is the value SMIN which is 
initially assigned equal to the input value SMAX. This initial value may be 
reduced in subroutine SLIMIT and BLOCKM. 

The subroutine operates on surfaces in the order of input starting with 
surface sequence number I BLOCK and proceeding to NBLOCK. Surfaces with sequence 
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index below IBLOCK are either the point-of-interest (1) or the plumes 
(2 to 1+NEUG). The dimensionless vector length (DELTA) is computed for 
each surface using Eq. 39. The minimum value of DELTA is maintained in 
DCOMP. After all surfaces have been tested, the value of SMIN is redefined as 
SMIN = DC0MP*SMIN 
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SUBROUTINE NAME: BLOKIN 

DESCRIPTION 

BLOKIN reads all surface geometry input and prepares the necessary 
transformation vectors and matrices for later use in other subroutines. 

The surface input includes the point-of-1nterest, plume origins and 
blocking surfaces. The sequence of the surfaces is defined as follows: 

1 Point-of-interest 

2 to 1+NEN6 Plumes 

IBLOCK to NBLOCK Blocking surfaces 

The parameter NENG is defined from the input for IFLOW (Card 2), and 
IBL0CK=2+NENG. The upper index, NBLOCK, is set by the program when surface 
input is terminated. 

CALLING SEQUENCE 

CALL BLOKIN 


COMMON BLOCKS AND EXTERNAL REFERENCES 


COMMON/CONTROL/ 

COMMON/CON ID/ 

COMMON/FLOSCL/ 

COMMON/GOMTRY/ 

COMMON/BLOCKO/ 

COMMON/BLOKOP/ 

COMMON/FLOWS/ 


COMMON/SLIMTC/ 

COMMON/CHI POF/ 

ATAN 

SIN 

COS 

TAN 
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METHOD OF SOLUTION 

The transformation matrices (T(J,K,I)) and vector components (R1(I), R2(I), 
R3(I)) are computed for each surface as described in Section 4.1 (Eqs. 30 and 
31). Values of conical parameters for the approximate plume boundaries are com- 
puted and surface coefficiences A,B, C (Eq. 38) and CC (Eq. 39) are defined for 
the plumes and blocking surfaces. 


A93 



f=? e: -r ^ o j— 1 IPv4 o- 


RTR 014-9 


SUBROUTINE NAME; FFPREP 
DESCRIPTION 


FFPREP prepares a binary file of flowfield gas properties on Unit 8 
using either card input or binary input on Unit 11. The input can be scaled 
if desired to convert dimensional units or change flowfield scale. A limited 
amount of output is provided to assist in diagnostics if the program encounters 
an error. 

CALLING SEQUENCE 

CALL FFPREP 

COMMON BLOCKS AND EXTERNAL REFERENCES 

COMMON/CONTRL/ 

COMMON/CON ID/ 

COMMON/FLOSCL/ 

COMMON/GOMTRY/ 

METHOD OF SOLUTION 

Gas property input is read, scaled if necessary and output to Unit 8. 


COMMON/ BLOCKO/ 
COMMON/ BLOKUP/ 
COMMON/FLOWS/ 
COMMON/SLIMTC/ 
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SUBROUTINE NAME; FLOUT 
DESCRIPTION 

FLOUT reads the files on problem description from Unit 10 and prepares 
printed output of the problem description. Following the output, all dimensioned 
parameters are converted to metric units and radius., 

CALLING SEQUENCE 


CALL FLOUT 

COMMON BLOCKS AND EXTERNAL REFERENCES 

COMMON/CONTRL/ 

COMMON/CON ID/ 

COMMON/ FLOSCL/ 

COMMON/COMTRY/ 

COMMON/BLOCKO/ 

COMMON/BLOKOP/ 

METHOD OR SOLUTION 

Data on the problem description is read from Unit 10 and printed. If the 
case is being run from a Unit 10 generated on a previous run, tests are made 
to. assure compatibility and diagnostics are generated if inconsistencies are 
noted. 


COKMON/FLOWS/ 

COMMON/SLIMTC/ 

COMMON/CHIPOF/ 

SCLOCK 

SKIP 

COS 
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SUBROUTINE NAME: FLOW 


DESCRIPTION 


FLOW controls sequencing of the set of overlay subroutines involved 
in generating line-of-sight property data on Unit 10. 

CALLING SEQUENCE 


CALL FLOW 

COMMON BLOCKS AND EXTERNAL REFERENCES 

COMMON/CONTRL/ 

COMMON/CON ID/ 

COMMON/FLOSCL/ 

COMMON/GOMTRY/ 

COMMON/BLOCKO/ 

COMMON/ BLOKOP/ 

COMMON/ FLOWS/ 

METHOD OF SOLUTION 

Calling of Subroutines FLOWIN, FLOUT, FL0W3D and FLOWAX are controlled 
by the input parameter IFLOW (Card 2). 


COMMON/SLIMTC/ 

COMMON/CUIPOF/ 

FLOWIN 

FLOUT 

FL0W3D 

FLOWAX 


AO 
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SUBROUTINE NAME; FLOWAX 
DESCRIPTION 

FLOWAX uses axisymnetrlc flowfield gas property data Input on Unit 8 
to prepare line-of-sight gas property output on Unit 10. 

CALLING SEQUENCE 


CALL FLOWAX 

COMMON BLOCKS AND EXTERNAL REFERENCES 


COMMON/CONTRL/ 

COMMON/CON ID/ 

COMMON/FLOSCL/ 

COMMON/GOMTRY/ 

COMMON/BLOCKO/ 

COMMON/BLOKOP/ 

COMMON/FLOWS/ 


COMMON/SLIMTC/ 

PLUMPT 

CPUTIM 

BLOCKM 

SLIMIT 

SIN 

COS 


METHOD OF SOLUTION 


Gas property data input from Unit 8 describes properties in an axisymmetric 
format with arbitrary radial positions (R) specified at selected axial positions 
(Z). Linear interpolations are made in the flowfield gas property data (as 
described in Section 4.3) to provide gas property data at each position on a 
line-of-sight defined by the input problem geometry. As each line-of-sight Is 
processed, the integration limits in s are redefined If necessary by subroutines 
SLIMIT and BLOCKM. 


A97 



F=9 e: rs/1 -r e: o I— I i rvj o. 


RTR 014-9 


After each line-of-sight is processed, it is checked for negative property 
values, temperatures outside the usual range of the band model parameter table, 
and zero temperatures. If all temperatures are zero or if the adjusted inte- 
gration limits indicate the line-of-sight does not intersect the plume flowfield, 
the line-of-sight sequence number, LSET, is set negative to indicate the line-of- 
sight missed the gas. 
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SUBROUTINE NAME; FLOWIN 
DESCRIPTION 

FLOWIN reads or controls other subroutines (FFPREP and BLOKIN) which 
read flowfield and geometry input. The data, along with that read in MAIN, 
are output as the initial records on Unit 10. 

CALLING SEQUENCE 


CALL FLOWIN 

COMMON BLOCKS AND EXT ERNAL REFERENCES 

COMMON/CONTRL/ 

COMMON/CONI 0/ 

COMMON/FLOSCL/ 

COMMON/GOMTRY/ 

COMMON/BLOCKO/ 

COMMON/BLOKOP/ 

METHOD OF SOLUTION 

Problem description and integration limits and intervals are read in this 
subroutine. If the flowfield properties are input on cards. Subroutine FFPREP 
is called to prepare the property data on Unit 8. Subroutine BLOKIN is called 
to read all surface geometry input. 


COMMON/FLOWS/ 

COMMON/SLIMTC/ 

COMMON/CHIPOF/ 

SCLOCK 

FFPREP 

BLOKIN 
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SUBROUTINE NAME; FL0W3D 
DESCRIPTION 

FL0W3D uses 3-dimensional flowfield gas property data input on Unit 8 
to prepare line-of-sight gas property output on Unit 10. 

CALLING SEQUENCE 



CALL FL0W3D 

COMMON BLOCKS AND EXTERNAL 

REFERENCES 

COMMON/CONTRL/ 

PLUMPT 

COMMON/CON ID/ 

CPUTIM 

COMMON/FLOSCL/ 

BLOCKM 

COMMON/GOMTRY 

SLIMIT 

COMMON/BLOCKO/ 

SIM 

COKiMON/BLOKOP/ 

COS 

COMMON/FLOWS/ 

ATAN2 

COMMON/SLIMTC 

SQRT 

METHOD OF SOLUTION 



Gas property data input from Unit 8 describe properties in a cylindrical 

coordinate system with arbitrary radial positions (R) specified at selected 

axial (Z) and angular (n) positions. Linear interpolations are made in the flow- 

field gas property data (as described in Section 4.3) to provide gas property 

data at each position on the lines-of-sight defined by the input problem geometry. 

As each line-of-sight is processed, the integration limits in s are redefined if 

necessary by Subroutines SLIMIT and Pf OCKM. 

TO 
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After each line-of-sight is processed, it is checked for negative 
property values, temperatures outside the usual range of the band model 
parameter table, and zero temperatures. If all temperatures are zero or 
if the adjusted integration limits indicate the line-of-sight does not 
intersect the plume, the line-of-sight sequence number, LSET, is set 
negative to indicate the line-of-sight missed the gas. 


AlOl 
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SUBROUTINE NAME: MAIN 
DESCRIPTION 

MAIN reads the problem title, control parameters, and integration limits, 
but the primary function is to control overlay for the upper level program 
segments. 

CALLING SEQUENCE 


NONE 


COMMON BLOCKS AND EXTERNAL REFERENCES 


COMMON/CONTRL/ 

EXIT 

COMMON/CON ID/ 

FLOW 

COMMON/ FLOSCL/ 

RAD 

COMMON/COMTRY/ 

VIEW 

COMMON/BLOCKO/ 



METHOD OF SOLUTION 

Reads the first 6 input cards and checks to see if the program limits 
for number of constituents or spectral intervals will be exceeded. Calls are 
made to the three segments FLOW, RAD, and VIEW, then the next case is read. 
When the next case title is blank, the system routine EXIT is called. 
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SUBROUTINE NAME: N060 

DESCRIPTION 

NOGO is an integer function which determines if an intersection between 
a line-of-sight and a blocking surface Is within the portion of the surface 
described by the y limits defined in Fig. 4.3. 

CALLING SEQUENCE 

FUNCTION = NOGO (X,Y,A,B) 

where X,Y are two coordinates of the point of intersection in the surface 
coordinate system, and A and B are the y limits, GMIN(I) and GMAX(I). 

COMMON BLOCKS AND EXTERNAL REFERENCES 

COMMON None ATAN2 

METHOD OF SOLUTION 

The angle of the intersection in a plane normal to the Z-axis is 
determined in degrees, and it is converted to a positive angle if it was 
in the 3rd or 4th quadrant. The value is then checked against the y-liraits 
for the surface to determine if the intersection is within the limits 
(NOGO=1) or outside the limits (N0G0=-1). 


A 103 



F=? ^ rs/i ”T ^ o hH I rvj o. 


RTR 014-9 


SUBROUTINE NAME: PLANCK 

DESCRIPTION 

PLANCK provides integrated values of the Planck function using a 
tabulated function. 

CALLING SEQUENCE 

CALL PLANCK (PNUl, PNU2, T, PLA) 
where PNUl=PNU(l)=Lower wavenumber limit (cm“M 
PNU2=PNU(2)=Upper wavenumber limit (cm“M 
T = Temperature (K) 

PLA = Integrated value of the Planck function (watts/cm^-sr) 

COMMON BLOCKS AND EXTERNAL REFERENCES 

COMMON None EXP 

METHOD OF SOLUTION 


The value of the spectral integral of the function 
PL(I) = IT N^d^^/oT** 

is determined for the lower (1=1) and upper (1=2) limits of wavenumber. In this 
function N^ is the Planck function, m is the wavenumber and a is the 
Stefan-Boltzman constant. The integrated value of the Planck function is then 
PLA (PL(2)-PL(D) 

For the range 0.05 <_ PL(I)/T 1, a table of values (Ref. 10) for the 

integrals are linearly interpolated, while outside this range, integrals are 
ba 

based on limiting approximations of the Planck function. For PNU(I)/T < 0.05 
N° - 1 3.7405 X 10-^^ 0)^ 

^ TT 1 /1-5Q70 ,.,/T * 
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and the integral is 

TT N°d/aT'* = 0.15284(w/T)*. 

o ^ ^ 

For PNU{I)/T > 10, 

N° = i. 3.7405 X 10-^^ 0)^ 

~ ir exp (1.43879 w/T) 
and the integral is 

ir f = 1.0-0.4585 exp(-1.4388 w/T)[((w/T 

+ 2.0851)o)/T + 2.8984)w/T + 2.0145] 
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SUBROUTINE NAME; PLUMPT 


DESCRIPTION 

PLUMPT determines the location in the plume coordinate system of a 
point location specified in the central coordinate system. If there are 
multiple axisynmetric plumes, the plume system chosen is the one in which 
the point is closest to the plume axis. 

CALLING SEQUENCE 

CALL PLUMPT (X,Y,Z.XW,YW,ZW,RW,N1,N2,NP) 
where X,Y, Z = coordinates of the point in the central coordinate system, 

XW, YW, ZW = coordinates returned by the subroutine in the plume system, 

RW = distance of the point from the plume Z-axis, 

N1,N2 = lower and upper limits of the surface sequence index, 

corresponding to the plumes {usually 2, 1+NEN6), and 
NP = surface sequence index of the plume the point is located in. 


COMMON BLOCKS AND EXTERNAL REFERENCES 


COMMON/CONTROL/ 
COMMON/CON ID/ 
COMMON/ FLOSCL/ 
COMMON/COMTRY/ 
COMMON/BLOCKO/ 


COMMON/BLOKOP/ 

COMMON/FLOWS/ 

COMMON/SLIMTC/ 

COMMON/CHIPOF/ 

SQRT 


AlO" 
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METHOD OF SOLUTION 

Initialize RW=10^® and NP=0., then enter a loop for M=N1,N2 and compute 
the location in the plume XN,YN,ZM and radius RN for each plume. If the ZN 
value is outside the range of Z currently stored for the flowfield» proceed 
to the next N, but if is not, test the radial position. If RM <_ RW assign 
this point and plume as the solution 

RW=RN, 

XW=XN, 

YW=YN, 

ZW=ZN, 

and NP=N, 

and continue the loop through the plumes in search of another solution with 
RN < RW. When the loop is complete, the solution for the plume closest to 
the point is returned. 
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SUBROUTINE NAME: RAD 

DESCRIPTION 


RAD controls sequencing of the set of overlay subroutines involved in 
radiation calculations. 

CALLING SEQUENCE 


CALL RAD 

COMMON BLOCKS AND EXTERNAL REFERENCES 

COMMON/CONTRL/ 

COMMON/CON ID/ 

COMMON/ FLOSCL/ 

COMMON/ CONTRY/ 

COMMOM/BLOCKO/ 

METHOD OF SOLUTION 

Suborutine ACDATA is called unconditionally. Since all Input cards 
may have been read at this point, control is returned to MAIN if an error 
is detected sc that the next case may be attempted. If IERR0R=0, Subroutine 
SL6 is called for the radiation prediction, and if spectral output is called 
for (NSPEK > 0), Subroutine SPCOUT is called. Note that there are two decks 
with the subroutine name SLG, so the subroutine used will depend on the deck 
name used in the mapping process (refer to Section Al.l). 


COMMON/ABCOEF/ 
COMMON/ RADSUB/ 
ACDATA 
SLG 

SPCOUT 


AlOP 
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SUBROUTINE NAME; SKIP 
DESCRIPTION 


SKIP renews the page for output and returns the line counter to zero. 


CALLING SEQUENCE 


CALL SKIP 

COMMON BLOCKS AND EXTERNAL REFERENCES 


COMMON/CONTRL/ COMMON/GOMTRY/ 
COMMON/CON ID/ COMMON/ BLOCKO/ 
COMMON/FLOSCL/ 


METHOD OF SOLUTION 


Not applicable. 
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SUBROUTINE NAME: SLG (DECK SLGMLT) 

DESCRIPTION 

SLG (SLGMLT) reads line-of-sight property data from Unit 10* predicts the 
transmissivity using an exponential line strength distribution (Eqs. 10 and 
12) with the Curtis-Godson approximation (Eq. 15) for inhomogeneous gas 
effects. Radiation is predicted using Eq. 14. See also SLG(YNGH20). 

CALLING SEQUENCE 

CALL SLG 

COMMOM BLOCKS AND EXTERNAL REFERENCES 


COMMON/CONTRL/ 

CPUTIM 

COMMON/CONID/ 

PLANCK 

COMMON/FLOSCL/ 

SQRT 

COMMON/GOMTRY/ 

EXP 

COMMON/BLOCKO/ 

ALOG 

COMMON/ABCOEF/ 

SIN 

COMMON/RADSUB/ 



METHOD OF SOLUTION 


Integration of Eq. 14 is performed using the numerical form 
(})£ tor J 

q/A = AN®(u),i)[x(a),i-l)-T(u,i)] sinS cose ^ A6 


where t(o), 0)=1 and AN°(co,i) is used in preference to N°(co,i)Aw to more 
accurately represent the Planck function when large wavenumber intervals 
are used. In addition. Ax is used as equivalent to (dx/ds)ds. Actual 
values of the line-of-sight increments indicated by i may be determined 


Ar 
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by either the input As value (DS) or a multiple of As (DSS) set by the 
temperature step criteria TDIFF (described in Section 3.4). 

Values of the transmissivity are determined using Eqs. 10 and 12 with 
Curtis-Godson approximations of the band model parameters indicated by Eqs. 15. 
The pressure (p), mole fractions (c), and temperature (T) are obtained from 
the line-of-sight property values read from Unit 10. The line half-widths are 
computed using Eqs. 8 and 9, and the absorption coefficient (F) and line density 
(1/d) are determined from tables prepared by Subroutine ACDATA. The tables are 
prepared to correspond to the exact values of the wavenumbers to be used in the 
integration to represent the center of the intervals. Linear interpolation is 
performed between temperature entries in the tables to obtain values at. the gas 
temperature. 


Alll 
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SUBROUTINE NAME: SLG (DECK YNGH20) 

DESCRIPTION 

SLG(YNGH20) reads 1 ine-of-sight property data from Unit 10, predicts the 
transmissivity and radiance using an exponentially- tailed-inverse line strength 
distribution with the Intuitive Derivative Method for inhomogeniety (Eq. 17). 
This method is for water vapor only and does not include Doppler broadening. 

See also SLG(SLGMLT). 

CALLING SEQUENCE 


COMMON/CONTRL/ 

CPUTIM 

COMMON/CON ID/ 

YAPPRX 

COMMON/FLOSCL/ 

SQRT 

COMMON/GOMTRY/ 

EXP 

COMMON/BLOCKO/ 

SIU 

COMMON/ABCOEF/ 

ALOG 

COMMON/RADSUB/ 



METHOD OF SOLUTION 


The method of solution uses Eq. 14 in the font 
Gf (|)f (i)f smin s . / 

q/A = E Z Z Z N®(w,s)[ Z sin0 coscj) A<|) A6 

szero s-As 

where the terms in the inner [] brackets are defined in Eq. 17. Property values 
for the integration are based on the input As increment (DS) which may be lengthened 
to a multiple of As (DSS) by the temperature step size parameter TDIFF (described in 
Section 3.4). These incremental lengths are usually too large for accurate integra- 
tion of Eq. 17 in regions of moderate to strong absorption, so the inner bracketed 
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term is integrated using a smaller step size determined by the input parameter 
IRAD and the strength of absorption represented by the initial value of the 
parameter DFO (=F c p As). 

The pressure (p), mole fraction (c), and temperature (T) used in Eq. 17 
are obtained from line-of-sight property data on Unit 10. The line half-width 
(Yc) is computed from Eq. 8, and the absorption coefficient (F) and line density 
(1/d) are determined from tables for water vapor prepared by Subroutine ACDATA. 
The tables are prepared to correspond to the exact values of the wavenumbers 
representing the centers of the integration intervals. Linear interpolation is 
performed between entries in the tables to obtain values at the gas temperatures. 
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SUBROUTINE NAME: SLIMIT 


DESCRIPTION 


SLIMIT adjusts the integration limits on each line-of-sight to correspond 
to the portion of the line within the conical surfaces approximating the plume 
boundaries. 

CALLING SEQUENCE 

CALL SLIMIT 

COMMON BLOCKS AND EXTERNAL RE FERENCES 

COMMON/CONTRL/ 

COMMON/CONID/ 

COMMON/FLOSCL/ 

COMMON/GOMTRY/ 

COMMON/BLOCKO/ 

METHOD OF SOLUTION 

Solutions for the intersection of the line-of-sight with the limiting 

cone for each plume are obtained from Eq. 39. If no solutions exist, both 

limits (SZERO and SMIN) are set at zero. If solutions exist for a plume, 

the solutions are assigned as the upper and lower limits based on tests on: 

(1) the direction of the line relative to the plume axis, (2) whether the 

origin of the line is inside or outisde of the cone, (3) the position of the 

line-of-sight origin relative to the cone apex, and (4) the relative values of 

the solutions. After all plumes are processed, solutions are compared to select 

the smallest for SZERO and the largest for SMIN. 

Ail4 


COMMON/BLOKOP/ 

COMMON/FLOWS/ 

COMMON/SLIMTC/ 

COMMON/CHIPOF/ 

SQRT 
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In the process of obtaining solutions, the range of the cone is limited 
to the Z-range of the plume by 3 limits {Fig. 4.3-surface Type 5) which are 
set in the flowfield subroutine FLOWAX or FL0W3D. 


Alls 



F? ^ IN/1 1“ ^ O HH I rs4 O. 


RTR 014-9 


SUBROUTINE NAME; SPCOUT 
DESCRIPTION 


SPCOUT provides printed output for spectral radiation data. The spectral 
radiation is the only parameter output for the integrated problem. All other 
parameters are for the final line-of-sight only . The subroutine also prepares a 
a binary file on Unit 4 for use by the plot program PLOTIT. 

CALLING SEQUENCE 

CALL SPCOUT 


COMMON BLOCKS AND EXTERNAL REF E RENCES 

COMMON/ABCOEF/ 
COMMON/RADSUB/ 
CPUTIM 
ALOG 
SQRT 

METHOD OF SOLUTION 


COMMON/CONTRL/ 
COMMON/CON ID/ 
COMMON/FLOSCL/ 
COMMON/GOMTRY/ 
COMMON/BLOCKO/ 


Spectral parameters produced in Subroutine SL6 are used for spectral 
output and to compute desired spectral parameters. The equations used are 
those for the Curtis-Godson approximation with an exponential line strength 
distribution which were used in Deck SLGMLT, so the parameters printed in the 
"long form" output (Card 2 Col. 27 = 2 or 4) are not correct for use with Deck 
YNGH20 or any other band model. Because of this, the code In SPCOUT will dis- 
able the long form output if the parameter IRAD used by Deck YNGH20 is greater 

A^ 


than zero. 
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SUBROUTINE NAME: VIEW 

DESCRIPTION 

VIEW prepares a polar plot on the printer depicting line-of-sight radiance 
as a function of elevation angle 6 and azimuth angle This is a visualization 
tool to assist in judging the correctness of the point and plume orientations and 
give an approximate indication of the areas of intense radiation. 

CALLING SEQUENCE 

CALL VIEW 

COMMON BLOCKS AND EXTERNAL REFERENCES 

COMMON/CONTRL/ 

COMMON/CON ID/ 

COMMON/FLOSCL/ 

COMMON/GOMTRY/ 

METHOD OF SOLUTION 

A table IFIELD(51,83) is used to represent 51 lines with 83 symbols per 
line. Symbols are assigned based on the results on each line-of-sight read 
from Unit 10. An averaging method is used to approximate a result when several 
lines-of-sight must be represented by a single symbol. But no interpolation is 
performed, so widely spaced lines-of-sight result in widely spaced symbols. 

With the maximum radiance on any line-of-sight as RADMAX (watts/cm^-sr) 
the symbols used have the following meanings with the line-of-sight radiance 
noted as N. 


COMMON/BLOCKO/ 

SIN 

COS 
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Symbol Index 


Symbol Meaning 


1 

2 

3 

4 

5 

6 

7 

8 
9 
9 

blank 

B 

0 

* 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


0.1 RADMAX 

0.2 


0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 

0-9 

Missed the gas 
Blocked 

0 < N < 0.1 RADMAX 


N < 


0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.0 


Polar grid marks 


RADMAX 


8 



^ K/l “T E O h-H I rvj O- 


RTR 014-9 


SUBROUTINE NAME: YAPPRX 
DESCRIPTION 

YAPPRX evaluates an approximation for y(x,r,q) described in Eq. 17i. 
CALLING SEQUENCE 

CALL YAPPRX(X,R.Q,Y) 
where X = TrXe(ti>.s) (defined by Eqs. 17b and 17d) 

R = r(w,s) (defined by Eqs. 17c and 17f) 

Q ® q(o),s) (defined by Eqs. 17e and 17g) 

Y - y(x,r,q) (defined by Eq. 17i) 

COMMON BLOCKS Af>ID EXTERNAL REFERENCES 

COMMON None SQRT 

YF 

METHOD OF SOLUTION 

The method of approximating y(x,r,q) is that described by Dr. Stephen 
Young in Ref. 7. 
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FUNCTION NAME; YF 
DESCRIPTION 

YF is the value of y(x,r,l) used in the approximation (Subroutine 
YAPPRX) of y{x,r,q) 

CALLING SEQUENCE 

FUNCTION = YF(X,R) 
where X = u Xg(w,s) (defined by Eqs, 17b and 17d) 

R = r(ui,s) (defined by Eqs. 17c and 17f) 

COMMON BLOCKS AND EXTERNAL REFERENCES 

COMMON None SQRT 

METHOD OF ANALYSIS 

The equation for y(x,r,l) (Ref. 17) is 

YF(X.R) = y(x,r,l) = C2r(l+x)+(l+r=^)(l+x^)^^]/C(l+x")^*(r+(l+x")^)^] 
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A5 PLOT PROGRAM 


The plot program is the single module PLOTIT along with all the 
external references for plotting subroutines supplied by the M5FC Library, 
The only input required is the binary file on Unit 4 prepared as output in 
Subroutine SPCOUT (described in Appendix A2, Table A2.2). Description of 
the plot parameters (Card 16) is read in Subroutine ACDATA and passed to 
PLOTIT by the binary file. 

Output consists of the plots, diagnostic messages described in Section 
A2.2, and messages identifying the beginning and the end of the program to 
verify execution. 



